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Abstract 

Various pion and photon production mechanisms in high-energy nuclear col- 
lisions at RHIC and LHC are discussed. Comparison with RHIC data is done 
whenever possible. The prospect of using electromagnetic probes to charac- 
terize quark-gluon plasma formation is assessed. 
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1. INTRODUCTION 



The production of photons in heavy ion collisions is rather complex and, in the standard approach, one 
roughly distinguishes three types of mechanisms: 

— 1) the photon is produced in the hard interaction of two partons in the incoming nuclei similarly 

to the well-known Chromodynamics (QCD) processes (QCD Compton, annihilation, bremsstrahlung) in 
nucleon-nucleon collisions. The rate is calculable in perturbative QCD and falls off at large transverse 
momentum, p^,, as a power law. Photons are also produced as decay products of hadrons, such as of 
7r°, 77, • • •, which are emitted in hard QCD processes and at large the decay photon spectrum is also 
power behaved; 

— 2) in the collision of two nuclei the density of secondary partons is so high that the quarks and 
gluons rescatter and eventually thermalize to form a bubble of hot quark-gluon plasma: at LHC the initial 
temperature of the plasma is expected to be of the order of 1 GeV. It is assumed that the plasma evolves 
hydrodynamically. Photons are emitted in the collisions of quarks and gluons with an energy spectrum 
which is exponentially damped but which should extend up to several GeV; 

— 3) the QGP bubble expands and cools until a temperature of 150 to 200 MeV is reached and a 
hadronic phase appears. As they collide the hot hadronic resonances (tt^ , p, uo) emit photons until the 
freeze-out temperature is reached. The typical energy of such photons ranges from several hundred MeV 
to several GeV. 

In this standard picture it is expected that the thermal production mechanisms will produce an 
excess of photons with an energy of a few GeV. In the following we critically discuss each stage of the 
production processes. One of the main problems is the fact that the theoretical predictions suffer from 
very large uncertainties at each step of the above scenario. 

The report is organized as follows: 
chap. 2: nomenclature of photons according to their production mechanisms; 
chap. 3: results of RHIC; 

chap. 4: experimental aspects of photon detection at LHC; 

chap. 5: inclusive photon and (background) vr*^ production from perturbative QCD; 
chap. 6: inclusive thermal photon and 7r° production; 

chap. 7: comparison of thermal and non-thermal photon and -k^ production mechanisms; 

chap. 8. preUminary studies on photon-jet and photon-particle correlations; 

chap. 9: theoretical considerations on non equilibrium effects; 

chap. 10: theoretical considerations on lattice calculations of dilepton rates. 

Useful information on nuclear cross sections on the one hand, and on luminosities, acceptances, etc., on 
the other hand is collected in two appendices. 

Chapters 3 and 4 contain the experimental part of the report. The PHENIX results on vr^ and 7 
production are first reviewed before turning to a discussion of the experimental capabihties of ALICE, 
CMS and ATLAS for detecting pions and photons. 

In Chaps. 5 to 8 quantitative studies of inclusive photon production are presented. Both "signal", 
i.e. direct photons, and "background" i.e. photons from decay of resonances, are discussed. These 
studies are carried out using presently available tools and models and, whenever possible, uncertainties 
in the predictions are given. Proton-proton, proton-nucleus and nucleus-nucleus collisions are treated in 
parallel. Predictions are made for the ratio 7dircct/7aii or 7dircct/7r'^ which determine if the extraction of 
a direct photon signal is feasible. The production of low mass lepton pairs at large transverse momentum 
is presented as a channel complementary to real photon production: it is based on similar dynamics of 
production but suffers from different backgrounds. 



Whenever possible we use two alternative models. On the one hand, the "standard" approach 
based on next-to-leading-order (NLO) QCD calculations to describe the hard processes together with a 
hydrodynamic model to describe the thermal evolution of the fireball produced in heavy ion collisions, 
and, on the other hand, the Dual Parton Model (DPM) which combines soft and hard leading-order 
(LO) dynamics and has been extremely successful in describing hadron-hadron scattering as well as 
fixed target nucleus-nucleus scattering. Surprisingly, predictions of the two models in nucleus-nucleus 
collisions turn out to be very similar despite the fact that the treatment of final state nuclear effects are 
quite different. 

Exploratory studies on photon-jet, photon-hadron and photon-photon correlations are presented 
and the comparison between proton-proton and nucleus-nucleus scattering is made. Only the LO ap- 
proximation is used and, for the latter two cases, only the signal is considered. 

We conclude from these studies that: 

- thermal photon production manifests itself by an enhancement in the inclusive photon spectrum at 
values below 10 to 15 GeV/c at the LHC. The shape of the transverse momentum spectrum may also be 
indicative of the production mechanisms; 

- the 7direct/7r° ratio in nucleus-nucleus collisions should be large enough to allow for the extraction 
of a direct photon signal, however the uncertainties on the predictions are large mainly due to the poor 
knowledge of the model parameters; 

- the lepton pair channel at large momentum transfer looks promising but further detailed studies are 
necessary to determine if the large background from uncorrelated pairs can be reliably subtracted; 

- correlation studies show characteristic changes of shapes in A-i-A collisions compared to pp collisions 
but, here again, further studies are necessary concerning the background. 

Many of the results presented are new, in particular, the comparison between two alternative mod- 
els as well as the extensive discussion on the uncertainties at the LHC (role of the initial conditions of 
thermahzation, chemical equihbrium vs. non equilibrium, • • •), the studies of the lepton pair spectrum as 
well as that of the correlations involving a photon. 

Chapters 9 and 10 are of a different nature and address more fundamental issues. Some of the 
basic hypotheses upon which the thermal production studies are based are critically analysed. 

In Chap. 9, the relevance of the finite life-time of the thermal system is adressed and it is shown that 

recent claims of a very large production rate of photons due to this finite life-time are not tenable since 
they are based on a defective modeling of the system. This is an original piece of work. Furthermore, 
ways to deal with this problem in a realistic way are sketched. 

In Chap. 10, the improved perturbative methods upon which the calculation of thermal photon rates 
are based are compared with the non-perturbative lattice based method in a simple example, namely the 
production of a static lepton pair. The two approaches seem to indicate a large discrepancy both in the 
magnitude of the rate as well as in the functional dependence on the lepton pair mass. The error bars in 
both predictions are very large however and it is premature to conclude if a real contradiction between 
the two results exists. This is a very interesting problem and, clearly, further studies are called for. 

This report is far from being the definitive work on photon, dilepton and pion production in heavy 
ion colhsions at the LHC. The results are sufficiently interesting however to motivate more detailed 
theoretical and phenomenological studies on these topics. 

2. PHOTONS AND PHOTONS 

We define the "inclusive" photon spectrum in the usual sense: it is the unbiased photon spectrum ob- 
served in a collision between two hadrons or a hadron and a nucleus or two nuclei. This spectrum is 
built-up from a "cocktail" of many components: 



- "prompt" photons which are produced, early in the coUision, in hard QCD processes. They 
are directly emerging from a hard process or produced by bremsstrahlung in a hard QCD process. The 
associated spectrum is power behaved and dominates at large transverse momentum; 

- "thermal" photons which are emitted in the collisions of quarks and gluons in the QGP phase 
or in scattering of hadronic resonances in hot matter; their spectrum is exponentially damped at large 
enough energy; 

- "decay" photons are decay products of hadronic resonances (essentially vr'^ and rj). These res- 
onances can be either produced in hard QCD processes (and the corresponding decay photon spectrum 
will be power behaved) or at the end of the thermal evolution of the system; 

- "direct" photons are the sum of "prompt" and "thermal" photons. They can be obtained exper- 
imentally by subtracting from the inclusive spectrum the contribution from the "decay" photons which 
constitute a reducible background. 

In heavy ion collisions, the aim is the extraction of the "thermal" signal: it can only be done 
by subtracting from the "direct" photon spectrum the contribution of the "prompt" photons (which are 
an irreducible background to direct thermal photons) calculated from theory. Therefore it is of utmost 
interest to correctly estimate the latter component. A prerequisite is a precise control of the photon 
production rate in proton-proton and proton-nucleus collisions. 

One also defines in the context of thermal production "soft" and "hard" photons. The "soft" 
photons have an energy much less that the temperature of the medium while the "hard" ones have an 
energy of the order of the temperature or larger. This terminology is somewhat different from that used 
in the context of perturbative QCD. Only hard thermal photons are of interest for phenomenological 
studied since soft ones are overwhelmed by the background. 

3. NEUTRAL PION AND PHOTON RESULTS FROM RHIC 

G. David 

The first three years of RHIC experiments brought spectacular results at an impressive pace. After 
producing the first Au-i-Au collisions at ^/s = 130 GeV lune 12, 2000, the accelerator delivered signif- 
icant integrated luminosity (Table Q and by QM'Ol (lanuary 2001) many exciting analyses were com- 
pleted and presented. Maybe the most intriguing observation reported was the large suppression of high 
Pj, neutral pions and charged hadrons in central collisions with respect to peripheral or pp collisions 1 1 1 ' 
scaled with the number of nucleon-nucleon collisions. In Run-2 RHIC operated at ^/s = 200 GeV^ 
producing both Au-i-Au |2l|3l and polarized pp collisions 0]. The suppression of high hadrons in 
central Au-i-Au collisions has been confirmed and the measurement extended to 10 GeV/c, while the 
pp data provided neutral pion spectra up to 13 GeV/c. Therefore, the nuclear modification factor could 
be established with measured in the same experiment. However, it remained an open issue whether 
the suppression is an initial state or final state effect. Proving its versatility in Run-3 RHIC delivered 
D-i-Au collisions |5l|6l (and once again polarized pp) which were analyzed extremely fast and the results 
were published less than three months after data taking. These results essentially ruled out initial state 
effects as cause of the high Pj, suppression observed in Au-i-Au collisions at RHIC energies. 

Meanwhile, few and only preliminary results on inclusive and direct photon production became 
public. This is understandable since the photon measurement is much more difficult than the (correlation- 
type) vr*^ measurement, and also because claiming any excess photons over the abundant background from 
hadron decays assumes that spectra of the contributing hadrons themselves (vr*^, rj, ...etc.) are known to 
high precision. First published photon results from RHIC are expected by the end of 2003, initially 
addressing the high Pj, region where photon identification is least problematic. 

'pp spectra were interpolated from measurements at lower and higher ^/s 
Full design energy for Au+Au 



One of the strengths of the RHIC program is a certain redundance within and overlap between 
experiments. In particular, photon and measurements - while mostly done in PHENIX with the 
electromagnetic calorimeters Q - are also possible both in STAR and PHENIX via photon conversion 
which serves not only as a cross-check but helps to extend the spectra to very low p^,. In addition, 
even within PHENIX there are two different electromagnetic calorimeters using different technologies 
and analyzed separately. The fact that one can make independent measurements of the same observable 
within the same experiment greatly increases the confidence in the ultimate results. 



Run/date 


Species 




Int. luminosity 


Submitted publications 


Run-1 










June-Sep 2000 


Au+Au 


130 GeV 




PRL 88, 022301 (2002) 


Run-2 










Sep-Nov 2001 


Au+Au 


200 GeV 


24/ib-i 


PRL 91, 072301 (2003) 
NPA 715 (2003) 683c 
NPA 715 (2003) 691c 


Dec 2001 - Jan 2002 


pp 


200 GeV 


0.15pb-i 


hep-ex/0304038 v2 (^PRL) 


Run-3 










Nov 2002 - Mar 2003 


D+Au 


200 GeV 


2.7nb~i 


PRL 91 072303 (2003) 


Apr 2003 - Jun 2003 


pp 


200 GeV 


0.35pb^i 





Table 1: Overview of RHIC runs as of June 2003. Integrated luminosity is given for PHENIX (which is typically the highest 
for the four RHIC experiments). In the last column we only list the final (submitted) publications on vr" and published QM'02 
preliminary results on photons. 



3.1 7r° spectra at RHIC 

One of the first and still most intriguing results from RHIC was the observation in Run-1 (^/s = 
130 GeV) that in central Au-i-Au collisions the yield of high p^, vr'^-s was strongly suppressed with 
respect to expectations from pp results at comparable energy^ scaled by the calculated number of binary 
nucleon-nucleon collisions (Fig. Q, although the same collision scaling described the peripheral data 
adequately [1 J. Despite the large errors the effect was significant (> 2.5a), but low integrated luminosity 
and an only partially instrumented calorimeter prevented PHENIX from exploring it past Pj, = 4 GeV/c. 
In Run-2 the combination of higher cm. energy {^/s = 200 GeV), much higher statistics, and a fully 
instrumented detector'* made it possible to extend the minimum bias vr*^ Pj, spectra up to 12 GeV/c, and 
the semi-inclusive spectra up to 6-10 GeV/c, depending on the centrality class (Fig.|2li. 

The deviation from simple scaling with the number of nucleon-nucleon collisions Ncou (or, more 
rigorously, with the nuclear overlap function Tab) is usually characterized by the nuclear modification 
factor Raa defined as 

o . ^ {l/N^^X)d^Ni\/dp,dy 

the expectation being that starting at some relatively low Pj. (1.5-2.0 GeV/c) which marks the approx- 
imate transition from "soft" to "hard" physics and where jets become dominant Raa reaches unity be- 
cause higher p^ particles are produced in incoherent, large momentum transfer parton-parton scatterings 

''since there are no pp data at y/s — 130 GeV, the reference spectrum was an interpolation of SPS and Tevatron neutral and 
charged pion results; note that the uncertainty on the reference spectrum is comparable to the total error of the data on Fig.Q 

''in Run-1 only 3 of the 8 PHENIX electromagnetic calorimeter sectors were instrumented and read out; in Run-2 all 8 
sectors were operational and included in the analysis 




Fig. 1: PHENIX results from 130 GeV Au+Au collisions 
(Run-1). The yields per event at mid-rapidity for {left) 
charged hadrons and {right) neutral pions are shown as a 
function of for 60-80% {lower) and 0-10% {upper) 
event samples, with the tt" results from the PbSc and PbGl 
analyses plotted separately. The error bars indicate the sta- 
tistical enws on the yield; the surrounding brackets indicate 
the systematic errors. Shown for reference are the yields per 
collision in N-l-N collisions, of charged hadrons and neutral 
pions respectively, each scaled up by (Ncoii) for the class. 
The bands indicate both the uncertainty in the N-l-N refer- 
ence and in the determination of (Ncoii)- 



Fig. 2: PHENIX results from 200 GeV Au-l-Au collisions 
(Run-2). Invariant yields of tt" at mid-rapidity are plotted as 
a function of for minimum bias and 9 different centrality 
selections (0-10% is the most central, 80-92% is the most 
peripheral). The yields are scaled for clarity. 



whose (small) probability in turn is proportional to the number of nucleon-nucleon collisions Ncoii- At 
lower ^/s (SPS-energies) Raa even rises above unity (Cronin-effect) due to multiple scattering of partons 
before the large process initiates the jet. In stark contrast to that expectation of Ncou scaling or some 
enhancement Raa > 1> PHENIX found already at ^/s = 130 GeV (Run-1) that for vr^ in central Au-i-Au 
collisions Raa never reaches unity. Instead, after reaching its maximum at Pj, ~ 1.5 GeV/c it decreases 
for higher transverse momenta (Fig.|3ll to about 0.3. In other words, there is a factor of ~3 suppression 
already around = 3 GeV/c. On the left panel of Fig. |5] is plotted for vr^ and charged hadrons 
((/i"*" + /i~)/2 in the most central Au-i-Au collisions with Ncou calculated from the Glauber-model. The 
suppression at higher p^, is even more dramatic when compared to the enhancement observed in Pb-i-Pb 
at ^/s = 17.3 GeV and a + q at = 31 GeV, also shown on the plot. On the right panel of Fig.|3]the 
ratio of central to peripheral data - both normalized with Ncoii - is shown. The behavior is very similar 
to Raa, but the coverage is somewhat less since the peripheral spectrum has a smaller range. It 
should be pointed out that Raa and the central/peripheral ratio (often referred to as Rc p) are dominated 
by different systematic errors, so the combined "message" of the two results (significant suppression) is 
even stronger than suggested by the error-bars on any one of them. Also, uncertainty on cJpp^' cancels to 
first order in the central to peripheral ratio (right panel). 

The observed suppression was a very strong indication that an extremely hot and dense medium 
has been created in central Au-i-Au collisions which by some mechanism depleted the number of high 
Pj, jets. However, it was unclear whether the very creation of the jets was suppressed {e.g. due to initial 
state parton saturation) or the jet production itself was the same as in pp, scaled by Ncoii, but the jets lost 
a significant fraction of their energy while traversing the medium (as predicted by different models of jet 
quenching). Also, in case of energy loss the nature of the medium causing the loss (partonic or hadronic) 
was unclear. Finally, various theoretical scenarios predicted similar suppression at 3-4 GeV/c (the up- 
per limit of the Run-1 results), but predicted different p^ dependence of the suppression at transverse 
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Fig. 3: PHENIX results from 130 GeV Au+Au collisions (Run-1). Left panel: the ratio Raa for charged hadrons ((/i"*" + 
h~)/2) and neutral pions in central Au+Au collisions. The error bars indicate the statistical errors, the surrounding bands 
(shaded for 7r"'s, brackets for {K^ + h^)/2) indicate the combined statistical and systematic errors on the ratio, including the 
uncertainty in thepp data and the uncertainty in (Ncoii)- Also shown for reference are Raa for a + a (y/s = 31 GeV) and for 
central Pb+Pb collisions {^/s = 17.3 GeV) measured at the CERN-SPS. Right panel: Ratio of central to peripheral p^, spectra 
(both normalized with the calculated Ncott) for charged hadrons and n". 



momenta beyond that range. 

In Run-2 RHIC delivered Au+Au and pp collisions at ^/s = 200 GeV, both at sufficiently high 
integrated luminosity such that the p^. spectra could be extended to 10-13 GeV/c. The pp results are 
described elsewhere in this Report and in |4|, but it should be emphasized once again that measuring 
the reference pp spectrum in the very same experiment {i.e. with similar systematics) greatly reduces 
the systematic error on Raa proper. The nuclear modification factor using vr'^-s in the most peripheral 
and most central Au-i-Au collisions is shown on Fig. 0] where this time Raa was calculated using the 
PHENIX pp measurement. Perfect scaling with Ncoii would mean Raa = 1- Although averaging below 
one, the peripheral data are certainly consistent with Ncou scaling within errors^. However, in central 
collisions the suppression is unambiguous in the entire p^ range^. Raa reaches its highest value (small- 
est suppression, ~2.5) around 2 GeV/c, then decreases, and is constant within errors above 4 GeV/c, 
giving a suppression factor of 4-5. While this result disfavored those quenching models that predict apj, 
dependence of the suppression^, it did not differentiate between initial state and final state effects^, nor 
could it distinguish between partonic and hadronic energy loss scenarios. 

Therefore, the major part of Run-3 at RHIC has been dedicated to D-i-Au collisions at ^/s = 
200 GeV, under the assumption that in D-i-Au collisions the gold nucleus remains "cold", and any effects 
due to a hot, dense medium in the final state would be absent. On the other hand if an initial state effect 
in the Au nucleus is responsible for the observed suppression, such suppression should manifest itself in 
D-i-Au collisions, too. Note, that the choice of deuterium (D-i-Au) instead of proton (jjAu) was motivated 
by purely technical reasons^, and it has been demonstrated that there is little if any difference in the 

^Errors are dominated by the fully correlated normalization error shown as a grey band at the left side of the plot. 
*Also, it should be pointed out that in the region of overlap there is a very good agreement between the Raa values 
from Run- 1 and Run-2 

For instance models including the LPM-effect predict that Raa oc y/p^ asymptotically 

^Although the presence of a hot, dense medium that causes energy loss was further supported by the observation that 
back-to-back jets virtually disappeared in central Au+Au collisions 1 3 1 

'Easier to collide because the magnetic rigidity of the two beams is closer to each other than for pAu 
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Fig. 4: PHENIX results from = 200 GeV. Nuclear mod- 
ification factor Raa {p^ ) foi" 

central (closed circles) 
and peripheral (open circles) Au-l-Au collisions. The error 
bars include all point-to-point experimental (pp, Au-l-Au) 
errors. The shaded bands represent the fractional uncertain- 
ties in (TauAu) and in the tt" yields normalization added 
in quadrature, which can move all the points up or down 
together (in the central case the shaded band shown is the 
fractional error for the first point). 



Fig. 5: Top: Nuclear modification factor 7?dA for 
(/i'*'+/i~)/2in minimum bias D-l- Au compared to Raa 
in the 10% most central Au-l-Au collisions. Inner bands 
show systematic errors which can vary with p^, and 
outer bands include also the normalization uncertainty. 
Bottom: Comparison of Rua for {K^ +h~)/2 and the 
average of the -k^ measurements in D-l-Au. The bar at 
the left indicates the systematic uncertainty in common 
for the charged and vr" measurements. 



physics of +DAu and pAu collisions at these energies 1 5 1 - an interesting observation in its own right. 
The nuclear modification factor for D+Au is shown on Fig.|5lfor charged particles and neutral pions. The 
suppression observed in central Au+Au collisions is clearly absent here^". This result indicates that the 
suppression in central Au+Au collisions is not an initial state effect, nor does it arise from modification 
of parton distribution functions in nuclei. Further analyses including a detailed study of the centrality 
dependence of R^a may shed further light on the mechanism of suppression. 

There is a substantial effort in PHENIX to extract the t] yield from Run-2 Au+Au data and to 
establish the asymptotic rj/ir^ ratio, but no results have been presented yet. STAR published spectra 
at ^ = 130 GeV using the Kg vr+vr" channel |8|, but only at < 1.5 GeV/c. Both r] and 
have TT^ decay channels and thus "feed down" into the spectrum but due to the (two- and three-body) 
decay kinematics this contamination (as well as contributions from higher mass mesons) is negligible if 
compared to current experimental errors ||2|. 

Using the published vr'^ data at 130 GeV and the preliminary 200 GeV 7r° data a first study of 
XT scaling has been presented at the Fall 2002 DNP meeting |9|. At both energies suppression of high 
p^ t:^ with respect to point-like scaling from pp collisions has been observed. If the effect is due to 
shadowing of the structure functions rather than a final state interaction with the hot medium, the yields 
at a given xt and centrality should exhibit the same suppression and the scaling exponent n{xT, 
should remain unchanged from pp to Au+Au collisions. The study found that, within systematic errors, 
XT scaling with n{xT, = 6.3 ± 0.6 applies for both peripheral and central collisions in the range 
0.025 <XT < 0.06. 

3.2 Inclusive and direct (non-hadronic) photons 

While 7r° spectra have already been published from all three RHIC runs, so far only preliminary results 
were shown on photons. Though counterintuitive at first, measuring photons - even the inclusive photon 
spectrum - is much more difficult than measuring tt'^-s. Neutral pions are usually reconstructed from a 
correlation'^ and - except for very high p^, and/or very low multiplicities - they are measured only sta- 
tistically, calculating invariant mass from all photon candidates in an event rather than trying to identify 
both decay photons from a particular vr^. Moreover, at low multiplicities {pp, D+Au, peripheral Au+Au) 
an accurate vr'^ spectrum can be extracted without any photon identification at all - just by calculating 
the invariant mass from pairs of all clusters in an event. Even in high multiplicity events good photon 
identification and effective rejection of other particles is not crucial'^. Also, since the true mass of vr'^ is 
known, the measurement is "self-calibrating" in the sense that any shift in the observed centroid and any 
widening of the mass peak corresponds to a calculable shift in the (apparent) transverse momentum and 
can be corrected for. Finally, many of the pions created outside the colUsion vertex'^ do not reconstruct 
with the proper invariant mass due to the mismatch between the true and the apparent opening angle of 
the two photons'^. Therefore, many of those "background" vr'^-s don't contribute to the raw yields, and 
don't have to be corrected for. 

The situation with the inclusive photon spectrum is entirely different. Obviously particle identifi- 

'"Xhese results have been published less than three months after the data were taken and the analysis included only a subset 
of all available data. Therefore, within the quoted errors the tt" data can not differentiate between RdA = 1 or some small 
Cronin-type enhancement. A new analysis using the entire dataset is underway, it is expected to have smaller errors and it will 
investigate the centrality dependence of RdA as well 

"invariant mass nij-y of two photons 

'^What's really important is to know the efficiency of the photon identification (if any) and the smearing of the photon energy 
due to overlaps with other particles - which may push m~f^, out of a reasonable invariant mass window causing a loss in the raw 
tt" count. However, even large number of hadrons mistakenly identified as photons (contamination) rarely cause any problems 

''in nuclear interactions with the detector material, i.e. as "classic" background, or in decays of long-lived hadrons, i.e. as 
irreducible "physics" background. 

'''Unless the direction of the photon is measured, e.g. with a preshower detector or by an e'^e" conversion, one has to 
calculate the "apparent" opening angle of all photon pairs under the assumption that they came from the vertex. 



cation and a very precise knowledge of its efficiency is essential. Usually there is no straightforward way 
to check the energy (and p^,) scale, although - depending on the slope - just 1% error on the energy can 
give up to 10% error on the cross-section. Unless the direction of the photon is measured, there is little 
if any rejection/identification of photons not coming from the collision vertex. 

While inclusive photon spectra (7inc) are important in their own right, the exciting new physics is 
hidden in the "direct" or "excess" photon spectrum (7dir), the difference of inclusive photons and photons 
from electromagnetic decays of final state hadrons (7dccay)- Since 7dccay is simulated using fits to the 
measured Tr^,rj, ... spectra, reliable hadron spectra are a prerequisite to the direct photon measurement 
and many (although not all) of the errors on hadron spectra will propagate into the direct photon spectrum. 

Finally, a good direct photon measurement will have to reveal not only the magnitude of photon 
excess over hadronic sources 7inc — 7decay but also make the decomposition of at least three overlapping 
spectra possible (hard scattering, plasma phase, final state radiation). Therefore, the errors both on shape 
and absolute normalization should be very small^^. Once this ambitious goal is reached, one can work 
backwards (starting with the highest available p^) in unfolding the contributions of the different photon 
production processes'^. Once a high quality 7dir = 7inc — 7dccay "direct" (excess) photon spectrum 
is available, starting at the highest transverse momenta one can try to fix the pQCD scale and improve 
upon fragmentation functions in pp as well as establish the effects of the medium in Au+Au. Next, 
this pQCD contribution has to be subtracted from 7dir- The result at a few GeV/c is expected to be 
dominated by radiation from the early plasma. (Cross-checks with D-i-Au collisions will decrease the 
uncertainties). Then once again one can go lower in p^, (i.e. later in evolution time) to look for radiation 
from the plasma - hadronic gas phase transition, and so on. Whereas the spectral shapes of different 
contributors are different, deconvolution of these spectra obviously requires unprecedented accuracy on 
the 7dir excess photon spectrum (which implies a comparable accuracy of the vr*^ spectra themselves). 

Preliminary results from PHENIX at ^/s = 130 GeV (presented at QM'02) are shown on Fig. 13.31 
for peripheral collisions (left panel) and central collisions (right panel). The plots show the ratio 7inc/ 7deca 
of the inclusive photon spectrum and the photon spectrum expected from final state hadron decays. A fit 
to the measured vr'^ spectrum is used to calculate 7dccay with a fast Monte-Carlo program; m^, scaling 
is assumed for r] (higher mass mesons are not included). The band indicates the systematic errors, com- 
pletely dominated by the uncertainties of the fit itself. Note that in Run-1 the vr" measurement extended 
only until p^ = 4 GeV/c with large en^ors (Fig.Q. Therefore, the fit was not very well constrained, both 
the absolute normalization and the shape can vary significantly. Photon identification is based solely on 
shower shape and timing - charged particle veto is not applied. Fluctuations of the points indicate that 
the la errors are not overestimated. Within errors the results are consistent with no photon excess over 
known hadronic sources, but obviously the sensitivity of the measurement is low. 

Preliminary results from PHENIX at ^/s = 200 GeV (presented at QM'02) are shown on Fig. 13.31 
for peripheral collisions (left panel) and central collisions (right panel). Different from the 130 GeV 
data, here the double ratio {jmeas / T^meas) / ilsim/T^sim.) plotted, where the subscript meas refers to the 
measured inclusive photon and vr" spectra, vr^j^ is a fit to the measured vr'^ spectrum which is then used 
in a Monte-Carlo to generate the photon spectrum ^sim expected from hadronic decays. In the absence of 
any non-hadronic sources this double ratio would be exactly one, which is clearly the case in peripheral 
collisions (left panel). One should appreciate the large, non-statistical fluctuations of the points even at 
low - indicating that the errors are once again not overestimated. Data do not extend as far in as 
for vr'^-s, to avoid large uncertainties of the vr*^ fit in the tail and beyond the measured range. The double 
ratio for central events (right panel) is still consistent with no excess within errors, but obviously exhibits 
a suggestive trend. 

'^Less than 10%, preferably -^5-6%. Results from WA98 provide a useful context: their final publication quotes a 
dependent 5.7-8.9% systematic error on the photon excess. To squeeze the errors below 5% is almost impossible in a real-life 
experiment. 

'^Although this will never be a purely experimental process: it will always rely to some extent on theoretical calculations 
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Fig. 6: Ratio of inclusive photon spectrum to the photon spectrum expected from hadron decays in = 130 GeV Au+Au col- 
lisions, in peripheral events (left panel) and central events (right panel). The shaded bands are the systematic errors dominated 
by the uncertainties of the (measured) vr" spectrum 



2 
1.8 
1.6 r 
1.4 r 
1.2 r 

1 '-■ 
0.8 r 
0.6 r 
0.4r 
0.2E- 

0, 



PHENIX preliminary 



PbGI, peripheral (60-80%) 
\|s^ = 200 GeV Au+Au 



5 6 
Pt (GeV) 




0.6 
0.4r 
0.2E- 




PbGI, central (0-10%) 
200 GeV Au+Au 



5 6 
Pt (GeV) 



Fig. 7: Measured over expected 7/7r° ratio in Au+Au collisions at y^=200 GeV in peripheral (left panel) and central (right 
panel) collisions. "Expected" means that all photons are assumed to come from hadron decays, and the double ratio eliminates 
certain systematic errors. A ratio of unity means that all photons come from hadron decays. The shaded boxes represent the 
estimated la systematic errors on the data points. 



The STAR collaboration also presented preliminary results on excess photons at QM' 02 fill. In 
this measurement photons converted to e+e^ pairs in the silicon vertex tracker and the inner field cage 
of the TPC were analyzed. The experiment found that above = 1.65 GeV/c the contribution from 
TT^ to the inclusive photon spectrum decreases in the most central events (0-11% centrality), indicating 
an increase in contribution from other photon sources, possibly other electromagnetic decays or direct 
photons. However, there is a significant uncertainty in the normalization of the tt^ spectra. 

3.3 Outlook 

Obviously, the carefully worded statements from both experiments reflect the technical difficulties and 
do not exclude a possibly substantial yield of direct photons. Results with reduced systematic errors can 
be expected within a year and at larger p^, will provide the first tests of pQCD calculations. However, as 
emphasized in other sections of this Report, measuring the photon excess over the yield expected from 
hadron decays is only the first step: much of the interesting physics lies is hidden in the composition of 
this excess. In order to disentangle the contributions from hard scattering, the QGP, the mixed and the 
hadronic phases one has to measure this excess both in a very wide range and with unprecedented 
accuracy. This is and continues to be the biggest challenge to the experimentalists. 



Therefore, it is instructive to briefly review the dominant sources of systematic errors on 7inc/7dccay 
in the current analyses. At RHIC the errors on vr*^ are currently around 14-17%, distributed about equally 
between yield extraction, particle identification and effects of the energy scale, and this determines the 
accuracy of the calculated 7decay- Contributions from other mesons are currently not measured (although 
analyses of t] as well as Kg above 1.5 GeV/c by charged pions are under way). 

For the inclusive photon analysis photon efficiency has about the same error as in the vr*^ measure- 
ment, but hadron contamination is added. The energy scale uncertainties have bigger weight (because of 
the lack of "self-calibration", the possibility to cross-check the energy scale using the invariant mass), 
and a major contributor to the errors is the instrumental background. Note that three of these contribu- 
tors (particle identification, contamination and background) become worse as one tries to move to lower 
transverse momenta (thermal region). 

On the other hand it is exactly at low p^, where the complementary measurements via conversion 
electrons in both PHENIX and STAR offer advantages: since they provide directional information, much 
of the instrumental background, and, in general, photons not from the vertex can be eliminated. Some 
upgrades of current RHIC detectors point in this direction, too. Also, new types of analyses are devel- 
oped to reduce systematic errors^^. While trying to defeat all sources of systematic errors is probably 
a futile excercise one can make independent measurements of the same quantity, within the same ex- 
periment, with very different systematics, thus increasing the level of confidence in the results. Photon 
measurements are notoriously difficult: the author thinks future LHC experiments are well advised to 
include such "redundancies" from the very beginning. 

4. DETECTOR STUDIES, RESOLUTION 

In this section one briefly describes the main features of ALICE, CMS and ATLAS concerning photon 
and hadron detection. The acceptances for each detector are summarized in appendix 2. 

4.1 Photon detection at ALICE 

H. Delagrange, D. d'Enterria, T. Peitzmann, Yu. Kharlov 

The ALICE detector 1 12 1 aims to study the physics of strongly interacting matter at extreme energy 
densities, where the formation of a new phase of matter, the quark-gluon plasma, is expected. Among 
all the probes of the quark matter, photons could be used to study the dynamics of strong interactions 
in hadronic collisions. Owing to their small electromagnetic coupling, produced photons do not interact 
with the surrounding matter and thus probe the properties of the matter at the time of their production 

4.11 Photon spectrometer PHOS in ALICE 

In ALICE, the photon spectrometer PHOS is designed |14| to detect, identify and measure with high 
resolution the 4-momenta of photons. Photon studies in heavy ion collisions require from the detec- 
tor a high discrimination power between photons and any other kind of particles, charged and neutral 
hadrons or electrons. The best possible resolutions in energy and position, as well as a reasonably large 
acceptance provide a high neutral meson identification and, thus, the high precision determination of the 
background for the direct photon spectrum. 

The final design of the photon spectrometer PHOS consists of five identical modules positioned at 
the bottom of the ALICE detector (Fig|8ll. The PHOS modules are positioned at the distance 4.6 m from 
the beam interaction point and installed at the azimuth angles ±40.7°, ±20.3°, and 0°. Each module 
consists of the electromagnetic calorimeter detector (EMC) and a charged particle veto detector (CPV). 

"One possibility is to compare results from analyses where the origin of systematic errors is very different, like comparing 
Raa and Rcp in the 7r° analysis. 



Fig. 8: The PHOS detector: PHOS inside the ALICE solenoid magnet (left), and one PHOS module (right). 



Each EMC module is constructed as a matrix of 64 x 56 cells of lead tungstate (PbW04) scintillator 
crystals. Each crystal, elementary unit of the calorimeter, is an 18 cm long parallelepiped providing 
20 units of radiation length (Xq = 0.89 cm). It is shaped with a squared cross-section of 22x22 mm^, 
to be compared to the Moliere radius of lead tungstanate, tm = 20 mm. The scintillation light, in the 
visible near UV- wavelength range, is read out by a 5 x 5 mm^ avalanche photo-diode (APD) integrated 
with a low-noise pre-amplifier. The calorimeter is operated at low temperature, —25 °C, stabilized to 
±0.3° C. This operation mode on one hand enhances the scintillation light output by a three fold factor 
and provides the required high and constant energy resolution even for the less energetic photons and, 
on the other hand, keeps the noise of the APD low enough to provide a high signal to noise ratio. The 
electronic chain associated to each crystal of the PHOS spectrometer delivers two energy signals, one 
with low- and one with high gain, proportional to the energy deposited in the crystal and a timing signal 
that measures the time of the particle impact with respect to a trigger time reference. 

The CPV consists of multiwire proportional chambers with cathode pad read-out. Each calorime- 
ter module is covered by a CPV module with an active area of 144.6x144.6 cm^, 1.3 cm deep and 
filled with a gas mixture 70% Ar and 30% CO2. The total thickness of the CPV module is 5.1 cm. Low- 
mass construction materials are used for the CPV construction to minimize the material budget, radiation 
length and detector mass. The anode wires of the proportional chamber are strung along the beam direc- 
tion with a wire pitch of 5.65 mm and placed at 5 mm above the cathode. The cathode is segmented in 
64 (along beam direction)x 128 (across beam direction) rectangular pads of 2.26x 1.13 cm^, elongated 
along the anode wires. The electron shower induced by the passage of a charged particle is collected on 
the cathode and an induced charge is retrieved from each pad of the CPV. 

The PHOS acceptance in pseudorapidity is defined by \ri\ < 0.13. Each of five modules covers 
17.8° in azimuth angle. 

4.12 Intrinsic performance of PHOS 

The two parameters that describe the response of the EMC spectrometer and play the most important 
role for photon identification are the resolutions in energy and position. The energy resolution depends 
on the spectrometer ability to collect most of the electromagnetic shower energy, convert it into visible 
light and transmit it to the APD, as well as on the APD photo-efficiency and photon electron gain factor. 
The position resolution depends on the segmentation of the spectrometer and energy resolution. 



To determine experimentally these features an electron beam of energy ranging from 0.6 to 4.5 GeV 
irradiated the central module of an array of 3 x 3 EMC modules. A Gaussian function was adjusted to 
the distribution of total energy, E, collected in the array. The resulting resolution, a/E, was compared 
(Fig. 121) to the one obtained by the simulation performed in exactly the same conditions as the experi- 




Fig. 9: Energy resolution of a 3 x 3 EMC array measured in response to mono energetic electrons (•) or calculated with 
simulations of mono energetic photons (A). The continuous line represents the result of the fit of Eq. to the data and 
extrapolated to 10 GeV. The dashed line represents the values quoted in the TDR 1141 . 

ment. The following parametrization was adjusted to the experimental resolution dependence on electron 
energy: 

where E is in units of GeV, a represents the contribution of the electronic noise, b the stochastic term, 
and c the constant term. 

The impact position on PHOS is reconstructed by calculating the position of the center of gravity 
of the reconstructed cluster. This position is further corrected for the incidence direction of the impinging 
photon. The previously discussed test beam measurements were extended to verify the influence of the 
photon incidence on the position resolution by tilting the array of EMC modules by 0, 3, 6 and 9°. Fig.fTOl 
shows the r.m.s. of the Gaussian fit of this distribution versus the photon energies from 1 to 50 GeV for 
several incidence angles. 

The CPV detector is sensitive to any particle which initiates an ionization process in the CPV 
gas volume. Therefore it will detect charged particles with almost any momentum. The only parameter 
which defines the response of the CPV, is the position resolution of the charged track passing through 
the detector. The effective spatial resolution of CPV was measured during beam-tests as ax = 0.138 mm 
(across the wires) and o"^ = 0.154 mm (along the wires). Fig.^Jillustrates the coordinate resolution of 
the CPV. 

4.13 Particle identification 

Particle identification in PHOS is based on three methods: 
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Fig. 10: EMC position resolution versus photon energy for 
the incidence angles a = 0, 3, 6 and 9° as well as for all 
possible incidence angles of photons emitted from the inter- 
action point. 
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Fig. 11: Difference between the exact coordinate of the charged particle impact on the CPV, and the reconstructed coordinate. 
The plot for the x-axis (across the anode wires) is fitted by a single Gaussian, and for the z-axis (along the anode wires) is fitted 
by a sum of two Gaussians. 



• Time-of-flight (TOF) of the particles from the interaction point to EMC which allows to discrim- 
inate Ught particles (photons and high-energy hadrons) from slow heavy particles (low-energy 
nucleons); 

• Charge particle rejection which is based on matching of the reconstructed points in CPV and EMC; 

• Shower shape analysis which is based on the knowledge of the shower shape produced by different 
particles in the calorimeter. 

The performance of the time-of-flight depends on the time resolution of EMC electronics. The 
time-of-flight of the photons from the interaction point to PHOS is 15.3 ns. Figure^Jshows the spectrum 
of photons compared to those of neutrons and antineutrons identified as photons by TOF with two TOF 
resolutions, 1 and 2 ns, in the most central Pb-Pb collisions, versus their reconstructed energy in EMC. 
The final time resolution for the EMC electronics is not selected yet, so this figure provides the guideline 
for the design of the TOF system. 

Charged particles can be rejected in PHOS if a reconstructed point in the EMC matches a CPV 
reconstructed point. Figure ^] shows the average deviation between the CPV and EMC reconstructed 
points for charged pions versus their reconstructed energy, in the ALICE magnetic field 0.5 T. 

Shower shapes in EMC can be characterized by several parameters. A set of 7 shower parameters 
has been chosen to identify photons in EMC and discriminate them from hadrons and vr'^-mesons: 



I 1 ns I I 2 ns I 




Fig. 12: Spectrum of photons, compared to those of neutrons and antineutrons identified as photons by the TOF criterion with 
two TOF resolutions, 1 (left) and 2 ns (right) in Pb-Pb HIJING events. 

• lateral dispersion which is a mean squared deviation of the fired cells from the shower center; 

• shower main axes, Aq and Ai which are eigen values of the shower tensor in the {x,z) plane; 

• shower sphericity which is a relative difference between Aq and Ai. 

• the core energy which is a shower energy within a radius of 3 cm around the shower center; 

• the largest energy fraction in a single cell; 

• shower cell multiplicity. 

These 7 parameters can be statistically correlated, and a set of 7 statistically independent parameters is 
found by diagonalizing the covariance matrix of the shower shape in this 7-dimensional space. 

Figure ^5 shows the purity of identified photons by all three identification criteria, as well as the 
contamination by antineutrons. Purity is a fraction of the reconstructed particles identified as photons, 
which are really photons, from the total number of reconstructed particles identified as photons. Con- 
tamination in defined as a fraction of antineutrons which are identified as photons, from the total number 
of reconstructed particles and identified as photons particles. The purity and contamination are shown 
for three definitions of photon quality: low, medium and high purity photons. 

Shower shape analysis allows to distinguish photons and vr^-mesons at high px, where both parti- 
cles produce a single shower in EMC. Figure shows the true identification probability of a photon at 
high pt, a misidentification probability of a photon as a vr*^ and the ratio of the misidentification prob- 
ability to the true identification probability, for low, medium and high purity photons. This plot clearly 
illustrates that photons can be distinguished from the vr'^-mesons in a wide dynamical range. 
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Fig. 13: RMS of the Gaussian fit of the EMC-CPV distance along the a:-axis (i.e. across magnetic field) and a-axis (along the 
beam) for the charged pions produced with a uniform px distribution versus the deposited energy. 
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Fig. 14: Photon purity and antineutron contamination to the identified photon spectrum in HIJING Pb-Pb collisions, for three 
definitions of photon quaUty: low (A), medium (■) and high (•) purity photons. 
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Fig. 15: Photon true identification probabifity, photon misidentification probability as tt" and the ratio of the misidentification 
probability to the true one, for three values of photon purities. 



4.2 Photon detection at CMS 

O.L. Kodolova, J.H. Liu, LP. Lokhtin, A.N. Nikitenko, LN. Vardanyan, P. Yepes 
4.21 CMS detector 

The Compact Muon Solenoid (CMS) is a general purpose detector designed primarily to search for 
the Higgs boson in proton-proton collisions at the LHC |15|. The detector is optimized for accurate 
measurements of the characteristics of high-energy leptons and photons, as well as hadronic jets in a 
large acceptance, providing unique capabilities for "hard probes" in both pp and AA collisions L16J . 
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Fig. 16: The longitudinal view of the CMS detector. 
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Detailed description of the detector elements can be found in the Technical Design Reports 1*1711181 
[^123. The longitudinal view of the CMS detector is presented in Fig.[^ The central element of CMS is 
the magnet, a 13 m long solenoid with an internal radius f« 3 m, which will provide a strong 4 T uniform 
magnetic field. The Att detector consists of a 6 m long and 1.3 m radius central tracker, electromagnetic 
(ECAL) and hadronic (HCAL) calorimeters inside the magnet and muon stations outside. The tracker 
and muon chambers cover the pseudorapidity region |r/| < 2.4, while the ECAL and HCAL calorimeters 
reach |?7| = 3. A pair of quartz-fiber very forward (HE) calorimeters, located ±11 m from the interaction 
point, cover the region 3 < |7/| < 5 and complement the energy measurement. The tracker is composed of 
pixel layers and silicon strip counters. The CMS muon stations consist of drift tubes in the barrel region 
(MB), cathode strip chambers in the endcap regions (ME), and resistive plate chambers in both MB 



and ME dedicated to triggering. The electromagnetic calorimeter is made of almost 76000 scintillating 
PbW04 crystals and the hadronic calorimeter consists of scintillator sandwiched between brass absorber 
plates. The main characteristics of the calorimeters are presented in Table|2l 



Table 2: Energy resolution, a/E, and granularity of the CMS calorimeters in the barrel (HB, EB), endcap (HE, 
EE) and very forward (HF) regions. The energy resolution is shown for the total energy of electrons and photons 
(ECAL) and transverse energy of hadronic jets (HCAL, HF). 



Rapidity coverage 


< 7? < 1.5 


1.5 <| 7? |< 3.0 


3.0 <| r] |< 5.0 


Subdetector 


HCAL (HB) 


ECAL (EB) 


HCAL (HE) 


ECAL (EE) 


HF 


a/E = a/VE^b 

a 

b 


1.16 
0.05 


0.027 
0.0055 


0.91 
0.05 


0.057 
0.0055 


0.77 
0.05 


granularity 

Ar] X A(p 


0.087 X 0.087 


0.0174 X 0.0174 


0.087 X 0.087 


0.0174 X 0.0174 
to 0.05 X 0.05 


0.175 X 0.175 



4.22 Photon triggering and identification 

Photon identification, measurement and triggering in Pb— Pb collisions with the maximum estimated par- 
ticle density, dN^/dy{y = 0) = 8000, have been studied {WTTl with a full GEANT-based simulation 
of the CMS calorimetry and a parameterization of HIJING |22| data for the background. 

The CMS electron/photon trigger algorithm developed for pp collisions, see Fig. is suitable 
for triggering on energetic photons produced in the heavy ion collisions. Programmable thresholds on 
the cluster variables used in the algorithm have to be tuned to make it efficient. Estimates of the photon 
trigger rates have been made with two Algorithm Vetoes, see Fig.ini the Hadronic Veto {H/E < 0.4) 
and Neighbour Et Veto (^^ Neighbours Et < 25 GeV). The rate of the single photon trigger is less 
than 1 (10) Hz for a 50 (20) GeV threshold. With a such a threshold, the trigger efficiency is close to 
100% for the 7+jet events useful for the off-line analysis 1211 . 

Apart from the trigger selections, single photon identification based on calorimeter isolation or on 
the use of calorimeter cells above a certain Et threshold (labeled "cell Et cut") was considered. The 
photon energy may be measured in a 5 x 5 crystals cell (the trigger cell size) centred on the crystal with 
the highest response. Such a cell contains about 97% of the photon energy. Identification may be based 
on a cut in the transverse energy, Ef°^, deposited in an area of 3 x 3 or 5 x 5 such cells, not including 
the central one. Disttibutions of Ef''\b x 5) and Ef°\2, x 3) are shown in Fig. ITsT a) and (b) for the 
energy from a Pb— Pb event deposited in the ECAL only and in the total ECAL+HCAL system. Only 
about 6% of the transverse energy in the isolation area is measured by the hadron calorimeter, reflecting 
the softness of the charged particle spectra. 

The "cell Et cut" criterion is another method of photon identification which has been studied. It 
requires no energy above a given threshold deposited in every cell of the area around the central cell 
containing the photon. The transverse energy distribution in the cell is shown for Pb— Pb events in 
Fig. fTSl c). A threshold of Et = 6.5 GeV has been chosen for this distribution. This has been applied in 
an area of 7 x 7 cells, not including the central 3x3 trigger matrix since the trigger criteria must still 
be optimized and applied separately. "Cell Et cut" gives us a rejection factor of ~ 2.7 for vr'^ decays 
compared to a 14% reduction of the single photon rate. 
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Fig. 17: The CMS electron/photon trigger algorithm. 



The photon energy resolution is degraded due to the large "pile up" contribution in heavy ion 
collisions. In a 5 x 5 crystal matrix, we have extra energy deposited, with ~ 1 GeV, as seen in Fig. fTSl d). 
so that for a 120 GeV photon, the 0.64% test beam resolution [23} will be degraded to 0.80% due to "pile 
up". This photon resolution is still much better than the jet energy resolution. This can be improved by 
using 3x3 crystals matrix for energy measurement. 

4.23 Photon reconstruction efficiency and resolutions 

The capability of the CMS ECAL to reconstruct photons in heavy ion collisions was investigated in 
several Et intervals using a standard electromagnetic cluster reconstruction algorithm implemented in 
CMS object oriented reconstruction package ORCA (version 6). The algorithm looks for crystals with 
energies above a certain threshold and creates a cluster in a 5 x 5 crystal matrix. The full GEANT-based 
simulation of the CMS calorimetry (CMSIM_125 package) responses on single photons and HIJEMG 
central Pb— Pb event as a background were used in the analysis for the barrel and endcaps. 

Photon reconstruction efficiency as a function of photon transverse energy in Pb— Pb and pp events 
is shown in Fig.^J The estimated photon reconstruction efficiency for Pb— Pb collisions appears to be 
high enough, > 80%, starting from Et ~ 15 GeV. The efficiency dependence on the pseudorapidity is 
rather weak. The work on improvement of the photon reconstruction algorithm to increase reconstruction 
efficiency in high multiplicity environment is in progress. The photon spatial resolutions, o"^ ~ o"^ ~ 
0.005, are practically the same for pp and Pb— Pb collisions and do not depend significantly on the 
transverse energy. Such spatial resolution is better than rj — if size of an ECAL cell, 0.0174 x 0.0174. 
However, the influence of Pb— Pb background on energy resolution is more significant, as shown in 
Figure |20l At Et = 10 GeV, the transverse energy resolution degrades strongly in Pb— Pb events 
relative to pp, from 2% to 10%. The difference decreases with increasing Et and becomes insignificant 
at Et ~ 100 GeV. 
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Fig. 18: Distribution of the transverse energy of Pb— Pb event deposited in trigger cells of 5 x 5 (a) and 3 x 3 (b) 
crystals, not including the central one. (c) The transverse energy deposited in one trigger cell. The solid histograms 
in (a)-(c) give the energy in the electromagnetic calorimeter while the dashed histograms give the energy in the 
ECAL+HCAL system, (d) A fit of the transverse energy deposited in the electromagnetic part of one trigger cell. 
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Fig. 19: Photon reconstruction efficiency in Pb— Pb (soUd) and pp (dashed) events. 
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Fig. 20: Photon transverse energy resolution in Pb— Pb (triangles) and pp (squares) events. 



4.24 Jet reconstruction 



A detailed description of the jet reconstruction procedure in heavy ion collisions using the sliding 
window-type jet finding algorithm which subtracts the large background from the underlying event and 
a full GEANT-based simulation of the CMS calorimetry can be found in the chapter on jets. The ef- 
ficiencies and background contamination levels are shown in Table |3j along with the transverse energy 
resolution for several values of jet transverse energy in central Pb— Pb collisions assuming dN^/dy{y = 
0) = 8000. The jet energy resolution is defined as a(£;^^'=°/£;|f°)/(£^f?™/£;|f'') where is the 

reconstructed transverse energy and £'|f'" is the transverse energy of all generated particles inside the 
given cone radius R. Starting at Et = 100 GeV, jets can be reconstructed with 100% efficiency and 
purity. The purity is defined as number of events with true QCD jets divided by the number of events 
with all reconstructed jets. 

Table 3: The jet purity, noise (contamination levels, false jets / generated jets) and transverse energy resolution 
(|r/| < 0.3, R = 0.6) in central Pb-Pb collisions with dN^/dy{y = 0) = 8000. 



Et min (GeV) 


Purity 


Noise 


a{ET)/ET{%) 


75 


0.88 ± 0.03 


0.083 ± 0.009 


17.8 


100 


0.97 ± 0.03 


0.01 1 ± 0.003 


18.4 


125 


0.99 ± 0.03 


0.004 ± 0.002 


16.8 


200 


0.99 ± 0.03 


0.001 ± 0.001 


12.7 



Although the jet transverse energy resolution is degraded by a factor ~ 2 in high multiplicity 
central Pb— Pb collisions compared to pp, the average measured jet energy in Pb— Pb collisions is the 
same as in pp. Thus pp interactions can be used as a baseline for heavy ion jet physics. 

It is important to note that the jet angular resolution, a^p = 0.045 and cr^ = 0.05 at E';^* = 100 
GeV, is still better than the r] — ip size of an HCAL tower 0.087 x 0.087. Thus the angular position of a 
hard jet can be reconstructed in heavy ion collisions at CMS with high enough accuracy for analysis of 
jet production as a function of azimuthal angle and pseudorapidity. 



4.3 Photon studies in the ATLAS detector 

H.Takai, S. Tapprogge 

The ATLAS detector is designed to study high pT physics in full luminosity proton-proton col- 
lisions at the LHC. Most of the detector subsystems will be available for heavy ion collisions as well. 
We are interested in the detection of photons from heavy ion collisions. The highly segmented electro- 
magnetic calorimeter will be heavily used for this studies. We report on early assessment of the detector 
capabilities in the heavy ion environment. 

4. 31 The ATLAS Detector 

The ATLAS detector is designed to study proton-proton collisions at the LHC design energy of 14 TeV 
in the center of mass. The physics pursued by the collaboration is vast and includes: Higgs boson search, 
searches for SUSY, and other scenarios beyond the Standard Model, as well as precision measurements 
of processes within (and possibly beyond) the Standard Model. To achieve these goals at full machine 
luminosity of 10^^ cm~^s~^, ATLAS will have a precise tracking system (Inner Detector) for charged 
particle measurements, an as hermetic as possible calorimeter system, which has an extremely fine grain 
segmentation, and a stand-alone muon system. An overview of the detector is shown in Figure l2Tl 




Fig. 21: The overall layout of the ATLAS detector 



The Inner Detector is composed of (1) a finely segmented silicon pixel detector, (2) silicon strip 
detectors (Semiconductor Tracker (SCT)) and (3) the Transition Radiation Tracker (TRT). The segmen- 
tation is optimized for proton-proton collisions at design machine luminosity. The inner detector system 
is designed to cover a pseudo-rapidity of | rj |< 2.5 and is located inside a 2 T solenoid magnet. 

The calorimeter system in the ATLAS detector that surrounds the solenoid magnet is divided into 
electromagnetic and hadronic sections and covers pseudo-rapidity | r/ |< 4.9. The EM calorimeter is 



an accordion liquid argon device and is finely segmented longitudinally and transversely for | rj |< 3.1. 
The first longitudinal segmentation has a granularity of 0.003 x 0.1 {Arj x A0) in the barrel and slightly 
coarser in the endcaps. The second longitudinal segmentation is composed of Arj x Acp = 0.025 x 0.025 
cells and the last segment Ar] x A(p = 0.05 x 0.05 cells. In addition a finely segmented (0.025 x 0.1) pre- 
sampler system is present in front of the electromagnetic (EM) calorimeter. The overall energy resolution 
of the EM calorimeter as determined in test beam measurements is 10%/\/E © 0.5%. The calorimeter 
also has good pointing resolution, 60mrad/ \fE for photons and timing resolution better than 200 ps for 
showers of energy larger than 20 GeV. 

The hadronic calorimeter is also segmented longitudinally and transversely. Except for the end- 
caps and the forward calorimeters, the technology utilized for the calorimeter is a lead-scintillator tile 
structure with a granularity of At] x A</) = 0.1 x 0.1. In the endcaps the hadronic calorimeter is imple- 
mented in liquid argon technology for radiation hardness with the same granularity as the barrel hadronic 
calorimeter. The energy resolution for the hadronic calorimeters is 50%/\/^ ® 2% for pions. The very 
forward region, up to = |4.9| is covered by the Forward Calorimeter implemented as an axial drift 
liquid argon calorimeter. The overall performance of the calorimeter system is described in l24l . 

The muon spectrometer in ATLAS is located behind the calorimeters, thus shielded from most 
hadronic showers, and has a coverage of | t] |< 2.7. The spectrometer is implemented using several 
technologies for tracking devices and a toroidal magnet system, which provides a field of 4 T strength 
to have an independent momentum measurement outside the calorimeter volume. Most of the volume is 
covered by MDTs (Monitored Drift Tubes). In the forward region, where the rate is high, the Cathode 
Strip Chamber technology is chosen. The stand-alone muon spectrometer momentum resolution is of the 
order of 2% for muons with -p^ in the range 10 - 100 GeV. 

The trigger and data acquisition system of ATLAS is a multi-level system, which has to reduce the 
beam crossing rate of 40 MHz to an output rate to mass storage of O(IOO) Hz. The first stage (LVLl) 
is a hardware based trigger, which makes use of coarse granularity calorimeter data and dedicated muon 
trigger chambers only. It has to reduce the output rate to about 75 kHz, within a maximum latency 
of 2.5 ;us. The High-Level Trigger (HLT) is composed of two stages, the second level trigger (LVL2) 
and the event filter (EE), where further reduction of the rate is achieved using algorithms implemented 
in software, making use of the full granularity and all sub-detectors. Eor LVL2, the Region-of-Interest 
(Rol) concept is used to reduce the amount of event data needed to only a few per cent. Eor heavy ion 
physics where an interaction rate of 8 kHz is expected for full luminosity Pb-Pb collisions, we expect to 
be able to record data with minimal trigger requirements, e.g. centrality trigger. 

The performance results mentioned have been obtained using a detailed full simulation of the 
ATLAS detector response with GEANT and have been validated by an extensive program of testbeam 
measurements of all components. 

432 ATLAS and Photon physics 

Photons in ATLAS are detected in the EM calorimeter. Early in the detector design it was decided to 
include the capabilities of detecting the Higgs boson through its decay into two photons. Therefore 
the calorimeter is highly optimized for high pT photon detection and good rejection of tt^'s. The EM 
calorimeter sampling is divided in fine strips of Arj x Acf) = 0.003 x 0.1 to aid in the vr*^ rejection, by 
identifying overlapping photon showers. Extensive simulations of the performance in proton-proton en- 
vironment have been carried out. Eigure|22l shows the results of this studies for photons with transverse 
energies up to 100 GeV. Heavy ion events are different in character when compared to high luminos- 
ity proton-proton runs. Although the multiplicity is higher there is no event pile-up. The underlying 
background to signals of interest, for a typical central collision HIJING event, comes from the soft par- 
ticles produced in the collision. The majority of the soft charged particle tracks {pT < 200MeV/c) are 
curled up in the solenoid magnetic field. At the face of the calorimeter the density of charged particles 
is approximately dN/dy ~ 2000. The neutral particles are in their vast majority low pT vr'^s and they 
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Fig. 22: Photon to 7r° rejection as function of 77. 
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Fig. 23: Jet rejection after photon selection cuts as function of jet Et for low and high luminosity proton-proton conditions. 



will tend to deposit most of their energy in the first compartment, if not in the calorimeter and solenoid 
cryostat walls (~ IXq) However, the second electromagnetic compartment could be relatively quiet and 
used for photon studies. 

The major challenge in the study of single photons with the ATLAS detector will be the mis- 
identification of vr^'s (or jets with a leading vr'^) as a photon. Studies performed for high luminosity 
proton-proton runs indicate good performance for single photon identification. This is shown in Figs.l22l 
and |23] Because the performance is due to the fine transverse segmentation coupled to a longitudinal 
segmentation, we do expect similar performance for heavy ions. Detailed simulation studies are under 
way. 

5. NON THERMAL PRODUCTION MECHANISMS 

P. Aurenche, F. Bopp, H. Delagrange, S. Jeon, P. Levai, J. Ranft, I. Sarcevic, M. Tokarev, M. Werlen 



For the discovery of the quark-gluon plasma using single photon production the p^, range of in- 
terest is roughly 1 GeV/c < p^, < 10 GeV/c: indeed it is expected that thermal photon production in 
heavy ion collisions will increase the production rate somewhere in this domain. It is therefore important 
to also understand the non thermal production mechanism in the same energy range. We shall consider 
both 7 and vr'^ production since the latter is a background to the former. To calculate the relevant rates 
the usual theoretical tools at our disposal are the next-to-leading logarithm (NLO) QCD calculations. 
However, the considered Pj, values are very small compared to the center-of-mass energy of the colli- 
sion and one is not far from the (small x) kinematical boundary where perturbation theory may not be 
reliable. We therefore supplement the QCD predictions with those from a model which includes soft 
physics dynamics as well as semi -hard physics: such a model (Dual Parton Model as implemented in 
Phojet/Dpmjet) has been successfully confronted with data over a wide energy domain. We review 
each approach in turn, before presenting phenomenological results. 



5.1 Theory : perturbative QCD approach at next-to-leading order (NLO) 

P. Aurenche, H. Delagrange, S. Jeon, P. Levai, I. Sarcevic, M. Werlen 



5.11 Proton-proton collisions 

The production cross section of a particle at large transverse momentum in perturbative QCD is well 
known and has the usual factorisable form: 

da / . , dz „ , , ,^ „ , , ,s „ , , ^ ■. da 



dp.dy 



/dz do^ 
dxadxb^F^^Jxa, M)F^^^{xb, M)Dc/c{z, M^)-^——{fi, M,M^), (3) 



where the functions F are the parton densities and q^"^^" is the hard cross section between partons a 
and b, in the hadrons A and B respectively, to produce parton c. The function D^^^ is the fragmentation 
function of parton c into particle C. In the case a photon is produced, an extra term has to be considered 
where the fragmentation function reduces to a 5(1 — 2) function: in this case the photon participates 
directly to the hard collision (c = 7) in contrast to the bremsstrahlung process where the photon is 
produced in the fragmentation of a quark or gluon c}^ The calculations have been carried out up to 
next-to-leading order [26, 27 1 in QCD (all functions F, D, a"^ " are known in the NLO approximation) 
but there remains a residual ambiguity related to the choice of the unphysical renormalization scale 
factorization scale M and fragmentation scale Mp . 

'^Both directly produced photons and bremsstrahlung photons are "prompt" in the sense of Chap. l2l and they contribute to 
the "direct" photon spectrum. 



For photon production extensive phenomenological studies have been carried out in proton/(an- 
ti-)proton scattering for ^/s from 20 GeV to 1.8 TeV and the situation was found to be rather con- 
fused I28ll29ll30t . Concerning results on fully inclusive photon production in pp and pp, the theory is 
in satisfactory agreement with all data from (fixed target) 20 GeV to (ISR) 63 GeV with the same set of 
parameters (all scales set around or slightly smaller). There is one exception: the E706 data |l28| 
(at 31.6 and 38.8 GeV on Beryllium, but corrected by the experimental group to be compared to proton- 
proton scattering) which are at least a factor 2 to 3 above the other data and which do not have the same 
p^ dependence. For vr'^ production NLO theory and data are in agreement as far as the shape of the 
spectrum is concerned but the data are systematically above theoretical predictions 1 30 1 (with large fluc- 
tuations in the normalization of experiments when compared to theory): one possible explanation is the 
poor knowledge of fragmentation functions in the dominant large z region which is hardly constrained 
by e^e~ data from which the fragmentation functions are mostly derived. Indeed a large variation in 
the theoretical rates is observed when using different fragmentation functions (BKK |31|, KKP |32|, 
Kretzer L33.I ). The recent RHIC data at 200 GeV |4| will certainly help clarify the phenomenology. 

In the present work, one will need the NLO calculations in a brand new kinematical regime: 
-y/i = 5.5 and 14 TeV and p^, > 3 GeV, corresponding to very small x values. "Large Pj," NLO 
calculations have never been tested before, in hadronic collisions, in this small x kinematical regime. For 
transverse momentum of 3 GeV, the typical x = 2p,^/y/s values are of the order of 10^^ and much less 
if forward/backward rapidities are explored. One may question the reliability of straightforward NLO 
calculations in this domain. This is where "recoil" resummation H^ll is important but further studies are 
needed since the present results are rather dependent on non-perturbative parameters and no convenient 
phenomenological calculations are available at present. As mentioned above the NLO predictions also 
suffer from the usual (factorization, fragmentation, renormalization) scale uncertainties. In the following 
we follow the usual (albeit arbitrary) practice to choose a common scale and let it vary between p^ /2 
and 2p^}'^ The uncertainties on structure/fragmentation functions will be probed by using different sets. 

A specific feature of photon production at very high energy is related to the fact that the brems- 
strahlung component becomes large and dominant at small x. However this component is not really 
under control: in particular the G ^ 7 fragmentation channel is very important but it is not constrained 
by previous data l35ll . This point was not relevant for lower energies because of larger x values or 
because it was eliminated by isolation cuts in the collider data. However, at LHC, for small p^ this will 
introduce a large uncertainty on prompt photon production. 

Among other possible uncertainties in the predictions, one should mention those related to the 
"intrinsic" . As will be seen below, the vr'^ spectrum obtained at RHIC turns out to be in very good 
agreement with NLO predictions using standard scales (of the order of p^) thus alleviating the need of 
introducing "k^. effects". 

To summarize: the main uncertainties in the predictions are, as usual, related to the choice of scale 
values and the choice of fragmentation functions, the latter being important even in the case of prompt 
photon production because of the importance of the bremsstrahlung mechanism at small p^ and large 
energy. As for the uncertainties associated to the structure functions they turn out to be relatively small 
(a few % only). These points will be illustrated quantitatively in the phenomenological sections. At 
a more fundamental level, we must admit that the NLO machinery is applied in a kinematical region 
where it may not be justified but we have no possibility to gauge the associated uncertainty until recoil 
resummation is understood, a problem which needs an urgent solution. 

"Different choices were also tried such as varying M — /i and independently in the range specified above but the 
phenomenology for ^/s = 200 GeV and above is not affected. 

^"For an alternative approach see the recent papers of the Budapest group 1361 . 



5. 12 Nuclear effects in pA collisions 

One explains, in Appendix I, how to relate a pA hard cross section to the corresponding proton-nucleon 
cross section [Eqs. (I135I) . (Il36t 1. The incoming proton undergoes multiple scattering on the nucleons, 
constituents of the nucleus. The number of collisions Ncoii depends on the value of the impact parameter 
b or equivalently on the "centrality class", a high centrality being obtained in collisions at small impact 
parameter. The Glauber model used to describe the multiple collisions is based on the eikonal approxi- 
mation (independent scattering) and it is assumed that the parton distributions of the nucleons, confined 
in the nucleus, are the same as those of the free nucleon. It is known however that nuclear effects mod- 
ify the partons distributions. The nuclear structure functions are measured in Deep Inelastic Scattering 
(DIS) of leptons on nuclei [.371. At small values of x, for x < 0.07, the nuclear structure function is 
found to be less than nucleon structure function scaled by A, exhibiting the so-called nuclear shadow- 
ing. As X grows, the nuclear structure function gets bigger than the nucleon structure function. This 
is known as anti-shadowing. The kinematic region of interest at LHC energies is the region of nuclear 
shadowing |38 1. 

To calculate the rate of a hard process in proton-nucleus collisions one therefore uses Eq. ^ to 
describe the proton nucleon collision, where one of the partonic distribution F^/Aix, M), say, is a nuclear 
structure function properly normalized. Due to nuclear shadowing we expect at LHC, a suppression of 
photon and pion production in proton-nucleus collisions compared to the nucleon-nucleon case. The 
modification of the parton distribution is written in a factorized form as 

F,/a(x, M) = 5,/a(x, M) F,,n{x, M) (4) 

where Fa/N{x, M) is the parton distribution function in a nucleon and Sa/A{x, M) is the parton shad- 
owing function. We assume here Fg^jAix, M) to be normalized to one nucleon in the nucleus. Recent 
parametrizations of the shadowing function of Eskola, Kolhinen and Salgado (EKS98), which are M ^ 
dependent, distinguish between quarks and gluons l39l l40l and are shown to be in very good agreement 
with the NMC data on dependence of Ff'^/Fg-' 1411 . Another parametrization of nuclear parton 
distributions has been given in |42|. A detailed comparison between these various sets is given in Eskola 
etal. L38J . 

In the infinite momentum frame where the nucleus is moving very fast, shadowing is caused by 
high parton density effects at small x. The small x partons have a large longitudinal wavelength and can 
spatially overlap and recombine. These recombination effects reduce the nuclear parton number densities 
and hence the nuclear cross sections. Working in this frame enables one to treat nuclear shadowing and 
parton saturation in nucleons on the same footing due to the identical physical mechanism involved 
in both. Anti-shadowing is due to longitudinal momentum conservation (momentum sum rule) in this 
frame. 

Even though there has been considerable amount of theoretical work done on nuclear shadowing 
and impressive progress made in understanding the physical principles of nuclear shadowing |41|, we 
are far from having a precise and quantitative description of nuclear shadowing. The scale dependence 
of the nuclear structure functions is even less understood due to the limited range of Q ^ covered in fixed 
target experiments. Also, shadowing of gluons is not well understood due to the fact that they cannot be 
directly measured in DIS experiments. The working assumption is that high parton density effects are 
negligible and DGLAP evolution equations are valid in which case the gluon distribution function can be 
obtained from the scaling violation of the F2 structure functions. This assumption, however, will break 
down at small values of x due to high parton density effects l43l and one will need to measure the gluon 
distribution function differently. 

Another nuclear effect that may be considered is that of the Fermi momentum in the nucleus. In 
some approaches this contributes an extra nuclear "k^" which compounds with the intrinsic "k^" in the 
nucleon, thus leading to an appreciable increase of the cross sections |36l. As there is no need of a 
nucleon "k^" to describe proton-proton collisions, we likewise neglect the small nuclear "k^" one may 



introduce for proton-nucleus collisions. 



5.13 Nucleus-nucleus collisions 

The effect of multiple collisions is treated according to the Glauber model [see Appendix I, in particular 
Eqs. (I135t . (I136I) 1. Besides nuclear shadowing discussed in the previous section, an important effect in 
A+A collisions is the medium induced parton energy loss effect |44l. Fast partons produced in parton- 
parton collision propagate through the hot and dense medium and through scatterings lose part of their 
energy l45l and then fragment into hadrons with a reduced energy. While a dynamical study of the parton 
propagation in a hot and dense medium created in a realistic heavy-ion collision and the modification of 
the hadronization is most desirable, there is a phenomenological model f46 1 that takes this effect into 
account. Given the inelastic scattering mean-free-path, Xa, the parton type a scatters n times within a 
distance AL before it escapes the system. The modified fragmentation function, zD^/^i^^ ^f) given 
in terms of the photon fragmentation function zD^i^{z, Mp) by l46l 



1 ^ r " 

D^/,{z, AL,Mp) = —Y, Pain) z'^D^^/aK, M^) +J2 ^iD^./gi^, M,) (5) 



where = z/{l — (XlILo ^i)/^T)' ~ zE^/ej and Pa{n) is the probability that a parton of flavor a 
traveling a distance AL in the nuclear medium will scatter n times. It is given by 

Pain) = (^^/^")" e-AVA.^ (6) 
n! 

and Cff = Yln=o Pain)- The first term in Eq. (|5Jl corresponds to the fragmentation of the leading parton 
a with reduced energy E'^ — Yli^=o ^^^1" n gluon emissions and the second term comes from the j-th 
emitted gluon having energy ei, where ei is the energy loss of parton a after j scattering. One should 
keep in mind, however, that Landau-Pomeranchuk-Migdal (LPM) effect in QCD has been derived for 
static scatterers 1451 . which may not be suitable approximation in case on hot QGP. One can study the 
effect of parton energy loss on prompt photon and neutral pion production at the LHC by considering 
the following cases of parton energy loss |47 48|: 1) constant parton energy loss per parton scattering, 
e"^ = const, 2) Landau-Pomeranchuk-Migdal energy-dependent energy loss, e ~ \/E?i and 3) Bethe- 
Heitler energy-dependent energy loss, e ° ~ E!^. Alternative parametrizations of parton energy loss 
in a hot medium will be considered in Chap.l8lf49^. 

The energy loss mechanism affects the production of pions and the bremsstrahlung production 
of photons, but not the production of photons directly emitted from the hard scattering process. One 
therefore expects a stronger reduction of the vr^ rate than the 7 rate when going from the nucleon-nucleon 
collisions to A-i-A collisions and the ratio 7prompt / should be enhanced in A-i-A collisions compared 
to pp collisions. 



5.2 Theory: the combined Pomeron and (LO) perturbative QCD approach 

F.W.Bopp, J. Ranft 

Since, as mentioned previously, the relevant kinematical domain of interest is at rather small p^ 
we turn now to the discussion of a model which includes a "soft" physics component based on Regge 
theory as well as a "hard" component based on lowest-order perturbative QCD. This model reproduces 
soft and semi-hard data on production of hadrons from fixed target energies until 2 TeV. The production 
of "prompt" photons is not implemented, while that of "decay" photons is possible since the radiative 
decays of hadrons are included in the model. 

^'it is shown that recently observed suppression of tt" production in Au+Au collisions at RHIC, which is found to increase 
with increasing from 3 GeV to 8 GeV, is compatible with parton energy loss, when = 0.06K f48l (see sec.l3.4>. 



5.21 Proton-proton collisions, the Monte Carlo Event Generator PHOJET 

Hadronic collisions at high energies involve the production of particles with low transverse momenta, 
the so-called soft multiparticle production. The theoretical tools available at present are not sufficient to 
understand this feature from first QCD principles and phenomenological models are typically applied in 
addition to pertubative QCD. The Dual Parton Model (DPM) |50| is such a model and its fundamental 
ideas are presently the basis of many of the Monte Carlo (MC) implementations of soft interactions. 

Phojet-1.12 l5ni52l is a modern DPM and perturbative QCD based event generator describing 
hadron-hadron interactions and also hadronic interactions involving photons. PHOJET replaces the orig- 
inal Dtujet model ll53l . which was the first implementation of this combination of perturbative QCD 
and the DPM. 

The DPM combines predictions of the large N^, Nj expansion of QCD |54| and assumptions of 
duality l55l with Gribov's reggeon field theory |56|. PHOJET, being used for the simulation of elemen- 
tary hadron-hadron, photon-hadron and photon-photon interactions with energies greater than 5 GeV, 
implements the DPM as a two-component model using reggeon theory for soft interactions and (LO) 
perturbative QCD for hard interactions. Each PHOJET collision includes multiple hard and soft pomeron 
exchanges, as well as initial and final state radiation. In PHOJET perturbative QCD interactions are ref- 
ered to as hard pomeron exchange. In addition to the model features as described in detail in l57l . the 
version 1.12 incorporates a model for high-mass diffraction dissociation including multiple jet produc- 
tion and recursive insertions of enhanced pomeron graphs (triple-, loop- and double-pomeron graphs). 

High-mass diffraction dissociation is simulated as pomeron-hadron or pomeron-pomeron scatter- 
ing, including multiple soft and hard interactions |58j. To account for the nature of the pomeron being 
a quasi-particle, the CKMT pomeron structure function ll59l with a hard gluonic component is used. 
These considerations refer to pomeron exchange reactions with small pomeron-momentum transfer, \t\. 
For large \t\ the rapidity gap production (e.g. jet-gap-jet events) is implemented on the basis the color 
evaporation model 160|. 

For hard coUisions PHOJET uses the LO parton structure functions GRV98(LO)||6D- All color 
neutral strings in PHOJET are hadronized according to the Lund model as implemented in Pythia ll62l 
l63l . No parton fragmentation functions are needed separately. 

Phojet has been extensively tested against data in hadron-hadron collisions [57 1. In a number 
of papers the four experimental LEP Collaborations compare many features of hadron production in 7 
-7 collisions to PHOJET, a rather good agreement is usually found. PHOJET has been checked against 
practically all data on transverse momentum distributions in pp and pp collisions from colliders [64|. In 
Fig-EUwe plot this comparison. Please note that the points in this Figure are from the PHOJET Monte 
Carlo while the data are represented by fines, fits to the data points. 

5.22 Collisions involving nuclei, the Monte Carlo Event Generator DPMJET-/// 

The Dpmjet-III code system ll68ll69l . is a MC event generator implementing Gribov-Glauber theory 
for collisions involving nuclei. For all elementary nucleon-nucleon collisions it uses the DPM as imple- 
mented in Phojet. Dpmjet-III is unique in its wide range of applications simulating hadron-hadron, 
hadron-nucleus, nucleus-nucleus, photon-hadron, photon-photon and photon-nucleus interactions from 
a few GeV up to cosmic ray energies. 

Since its first implementations fTOl ItTI l72l DPMJET uses the Monte Carlo realization of the 
Gribov-Glauber multiple scattering formalism according to the algorithms of f73l and allows the cal- 
culation of total, elastic, quasi-elastic and production cross sections for any high-energy nuclear colli- 
sion. This formulation of the Glauber model is somewhat more detailed than the model described in 
Appendix I. In the model 1731 the scattering amplitude is parametrized not only by the inelastic nucleon- 
nucleon cross-section, but it is parametrized by using utot, P = R'e/(0)/iAr/Im/(0)/iAr and the elastic 
slope a. (Ttot and a are taken as fitted by PHOJET while for p a pai^ametrization of experimental data is 
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Fig. 24: Transverse momentum distributions of charged hadrons. The resuhs of Phojet (points) are compared to experimental 
data represented by hnes, fitted to the data points. The data at ^/s = 23 and 53 GeV are from the CERN-ISR |65 1, the data at 
200 GeV are from the UAl Collaboration ^ and the data at 630 and 1800 GeV are from the CDF Collaboration i67j. 

used. However, parameters needed for the collision scaling of the NLO vr'' and 7 cross sections (ainei, 
Ncoii) are in very close agreement in DPMJET with the ones determined in Appendix I. To be consistent 
we use in direct comparisons between DPMJET and NLO results always aimi and Ncou as determined 
by Dpmjet. No collision scaling is used by DPMJET, but of course it is easy to calculate Ncou in DPM- 
JET. Realistic nuclear densities and radii are used in DPMJET for light nuclei and Woods-Saxon densities 
otherwise. 

During the simulation of an inelastic collision the above formalism samples the number of "woun- 
ded" nucleons, the impact parameter of the collision and the interaction configurations of the wounded 
nucleons. Individual hadron-nucleon interactions are then described by PHOJET including multiple hard 
and soft pomeron exchanges, initial and final state radiation as well as diffraction. 

As a new feature, DPMJET-III allows the simulation of enhanced graph cuts in non-diffractive 
inelastic hadron-nucleus and nucleus-nucleus interactions. For example, in an event with two wounded 
nucleons, the first nucleon might take part in a non-diffractive interaction whereas the second one scatters 
diffractively producing only very few secondaries. Such graphs are prediced by the Gribov-Glauber 
theory of nuclear scattering but are usually neglected. Further features of DPMJET-III are a formation 
zone intranuclear cascade f74l and the implementation of certain baryon stopping diagrams II75I . 

The Dpmjet-III code and further information are available from the authors. DPMJET-III and 
earlier versions like DPMJET-II have been extensively tested against data in hadron-nucleus and nucleus- 



nucleus collisions 11681 1691 1701 1721. The code is used for the simulation of cosmic ray showers L76J . 

The transverse momentum distributions according to DPMJET are largely determined by the prop- 
erties of Phojet, which is called for each elementary interaction in DPMJET. Note, all of these elemen- 
tary interactions are treated kinematically correctly. All DPMJET events conserve energy, momentum as 
well as additive quantum numbers like charge, baryon number and strangeness strictly. Because of an 
appropriate treatment of soft physics the transverse momentum distributions start at p^, = 0. 

Applying DPMJET to heavy ion collisions at RHIC the first experience was, that the original 
hadron multiplicities and pseudorapidity distributions were about one third too high compared to the 
data. A new mechanism was needed to reduce N^ii and dNch/dr]\ri=o in situations with a produced 
very dense hadronic system. Such a mechanism, the fusion and percolation of soft chains, had indeed 
been conjectured since the beginning of the Nineties. Introducing percolation and fusion of soft chains 
into Dpmjet I77I . it was indeed possible to get a satisfactory agreement of DPMJET with the A'ch and 
dNch/dr] measured at RHIC. The procedure desribed in [77] is quite time consuming and therefore not 
well suited for high statistics DPMJET runs as needed here. Therefore, we apply here only an effective 
method, reducing A'ch and dNch/dr] in DPMJET like in [77J without spoiling energy-momentum con- 
servation. This method would not distort the transverse momentum distributions, but it will reduce the 
cross sections at all Pj, by approximately the same 30 % as mentioned above. 

Calculating from DPMJET transverse momentum distributions of decay photons (DPMJET or PHO- 
JET do not calculate direct photon production) we use the possibility in the code to declare certain reso- 
nances as stable. In this way we get separately (i) all decay photons (having vr'^ and t] decaying), (ii) all 
decay photons not coming fron vr^ decay (declaring vr'^ as stable) and (iii) all decay photons not coming 
from vr'' or r] decay (declaring vr'^ and r] as stable). 

5.3 The z— scaling model 

M. Tokarev 

During the workshop, a phenomenological model to describe hadron and photon production was 
also discussed: the z— scaling model. This concept was suggested to analyze numerous experimental 
data on high-pj, hadron production in pp,pp and pA collisions l78l . The method of data analysis was 
developed for description of pp, pp and pA interactions with direct photon and jet production in the 
high-py region at the RHIC and LHC energies as well f79ll . The scaling function is expressed via the 
invariant differential cross section Ed ^a/dp^ as follows 



TTS 



i:{z) = --^J-^E—. (7) 
pain dp'^ 

Here s is the collision center-of-mass energy squared, crj„ is the inelastic cross section, J is the corre- 
sponding Jacobian, and p = dN/drj is the particle multiplicity density. The function ip{z) is normalized: 

ij{z)dz = 1. (8) 

The relation allows us to give the physical meaning of the scaling function tlj{z) as a probability density 
to form a particle with the formation length. The variable z can be written in the form 

z = zon-^, where X2) = (1 - (1 - X2)'^^ (9) 

The factor zq is proportional to the transverse energy released in the underlying collision of constituents. 
The expression Q,^^ describes the resolution at which the collision of the constituents can be singled out 
of this process. The Q{xi,X2) represents the relative number of all initial configurations containing the 
constituents which carry fractions xi and X2 of the incoming momenta. The and 62 are the anomalous 
fractal dimensions of the colliding objects (hadrons or nuclei). The momentum fractions xi and X2 are 




a) b) 

Fig. 25: (a) Dependence of the inclusive cross section of vr" meson production on transverse momentum in pp collisions at 
^/s = 30 — 200 GeV. Experimental data are obtained at ISR and RHIC. (b) The corresponding scaling function. 



determined in a way to minimize the resolution Q ^(xi, X2) of the fractal measure z with respect to all 
possible sub-processes which satisfies 4-momentum conservation law. 

Figure |25] shows the data Pj, and z-presentation for vr'' meson production in pp collisions. The 
first one demonstrates the strong energy dependence of the cross section as a function of the transverse 
momentum. The z-presentation reveals the energy independence of ipiz) and, for large enough z, the 
power behavior, Tp{z) ~ z~^, of the scaling function. Based on the obtained results f78ll79l we can con- 
clude that verification of the asymptotic behavior of 'tp{z) for vr*^ meson, direct photon and jet production 
at LHC energies is of interest. Since no predictions are available for A-i-A collisions we do not pursue 
the z— scaling approach further in this report. 



In the next sections of this chapter we compare the model predictions to the vr transverse momen- 
tum distributions obtained at RHIC at y/s • 200 GeV, and then make predictions for vr" and 7 spectra 
at LHC energies. It is found that both DPMJET and NLO QCD agree with the data and, furthermore, 
they give very similar results for the extrapolation to LHC energies. A word of warning is necessary at 
this point concerning the comparison between the two models in nucleus-nucleus collisions . We note 
that the DPMJET predictions represent the full picture of a heavy ion collision while the NLO QCD pre- 
dictions give a partial picture to which the contribution of the quark-gluon plasma should be added (see 
Chap. 1^. In Chap. ^\it is found that thermal contributions, in A+A collisions, are important in the 
low Pj, range only. Comparing the predictions of the two models for A-i-A collisions at LHC is therefore 
relevant only for Pj, > 10 — 15 GeV/c. A detailed discussion of our results in A-i-A collisions below this 
value ofpj, is given in Chap. 



5.4 Phenomenology of vr^ production 

Our standard NLO predictions are obtained using the CTEQ5M fSOl or, equivalendy, CTEQ6M |8ll 
parton distributions and the KKP |32 1 fragmentation functions. In Fig.|26l the experimental vr*^ transverse 
momentum spectrum in pp collisions at 200 GeV [4J is compared to theoretical predictions. One notes the 
remarkable agreement of the data with the NLO QCD estimates as well as with the Dpmjet predictions. 



The Dpmjet model parameters had been tuned to fit charged particle spectra (see Fig. l24Ti and no new 
adjustement was necessary to obtain the vr*^ spectrum. 

The spread in the NLO QCD estimates (grey band in the figure) reflects the uncertainties associ- 
ated to the choice of scales (all scales equal varying from .5 p^. to 2 Pj.): if the uncertainty is roughly 
[+100%, -60%] at = 3 GeV/c it reduces to ± 40% at p^ = 12 GeV/c. The upper limit of the 
NLO predictions is obtained with the smallest scales and the lower limit with the largest scales. Vary- 
ing the scales independently (for example fi = M ^ ) one still obtains results in the grey band of 
the figure with only slight changes in the slope in p^, (for p^ below 6 GeV/c). Using the BKK l3Tl 
fragmentation functions rather than the KKP lowers the theoretical predictions by about 25% while the 
Kretzer |33 1 functions reduce the predictions by about 70% (lowest p^) to 45% (highest p,^). Concerning 
CTEQ6M, the 40 different sets predict cross sections within [—5%, +3%] for < 10 GeV/c and within 
[-6%, +10%] at p^ = 14 GeV/c. Turning now to MRS99 El rather than CTEQ hardly makes any 
changes: decrease of the cross section by about 10% at low p^ and stability for Pj, > 6 GeV/c. 
From this discussion one concludes that the largest uncertainties in the theoretical predictions are re- 
lated to the choice of scales and to the fragmentation functions, a remark which also holds true at LHC 
energies. Concerning the fragmentation functions, it should be recalled that KKP give the best phe- 
nomenology for vr'^ production at lower energies l83ll84l . 
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Fig. 26: PHENIX data on tt" production in pp collisions at 200 GeV, compared to various theoretical predictions. The grey 
band indicates the range the NLO QCD predictions for scales ranging from 0.5 Pj. to 2 p^ and the open points are the DPMJET 
predictions. CTEQ6 1811 structure and KKP [321 fragmentation functions are used. The theoretical predictions are evaluated 
by averaging the invariant cross sections in an interval ±0.35 unit of rapidity around rapidity 0. 

We also show on a linear plot, in Fig. ^\ the data normalized to the theoretical predictions. The 
solid line shows the ratio of a fit of the PHENIX data to the standard NLO predictions while the dashed 
line, with open crosses, shows the data normalised to the NLO estimates: the agreement between theory 
and experiment is remarkable given the fact that no parameters has been adjusted to fit the data. In 
view of the high quality, and large Pj, coverage, of the RHIC present and forthcoming data it would 
be interesting to perform a joint fit of e+e~ data and RHIC data to better constrain the fragmentation 
functions specially in the high z region where z is the fragmentation variable. Another advantage would 
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Fig. 27: PHENIX data on n production in pp collisions at 200 GeV normalized to various theoretical predictions shown in 
Fig. l26l The experimental error bars on the data points (crosses) are not shown. 



be the possibility to obtain a better control on the gluon fragmentation function which appears as a next- 
to-leading effect in e+e^ while it is a leading term in pp collisions: this would help un-correlate the 
gluon fragmentation parameters from the other parameters in the fit. The agreement between data and 
Dpmjet is also excellent even at the lowest Pj, values. The large fluctuations at high Pj, are related to 
the low statistics of the model. 

We turn now to LHC energies. The theory predictions are shown, for pp collisions, at 5.5 TeV 
(Fig. l28l and 14 TeV (Fig. l29l . In each case we display two plots, the first one with log-log scales to 
emphasis the low p^, region and the other one with semi-log scales up to p^, = 100 GeV/c. One notes 
again the good overall agreement between Dpmjet and NLO QCD results in the overlapping region. 
Furthermore the use of NLO calculations appears reasonable down to = 3 GeV/c below which one 
does not dare make perturbative QCD predictions ! In fact, at = 3 GeV/c, one notes a change of slope 
in the Dpmjet predictions, illustrating the effects of soft physics (pomeron exchanges) in the model: the 
slope of the distribution decreases leading to a finite prediction at p^, = 0. 

The scale variations in the NLO QCD predictions are still rather large at low transverse momentum 
(for = 3 GeV/c it is [+60%, -25%] at y/s = 5.5 TeV and [+25%, -20%] at ^ = 14 TeV) but 
reduces to roughly ±10% at = 100 GeV/c for both energies. These NLO results tend to be slightly 
less steep than the Dpmjet prediction and larger even when choosing the largest scales. This is due 
to using different input parameters for the perturbative sector in the two calculations. Keeping with the 
CTEQ family of 40 parametrizations the NLO predictions agree within ±10% at p^ =3 GeV/c and 
within ±2.5% at pj. = 100 GeV/c. As mentioned above using MRS99 and BKK fragmention functions 
would lower the NLO predictions which would then overlap the Dpmjet ones. A power behavior of the 
predictions holds for p^ between 3 GeV/c and 20 GeV/c: one finds for Dpmjet a behavior of type 
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Fig. 28: Theoretical predictions for the 7r° transverse momentum spectrum at 5.5 TeV. The conventions used are as in Fig. 1261 
The invariant cross section is calculated at rapidity. 



> 

^ 10 
£ 

n Q, 
■Ol T3 



7c° - =14 TeV - p + p Collisions 



10 - 



10 



10 E 



10 



10^ = 



10 



~i 1 1 — r 



~1 1 1 1 1 — I — I — 

DPMJET J. Ranftetal. 



NLO pQCD : CTEQ6M, KKP 

M=(i=Mp= 1 X [Jy 

M=(i=Mp : from 0.5 xppto2.0xpp ~| 



_l I I L. 




10 



Pj (GeV/c) 



10 

^ 10^ 

:s 10 
^ 1 

10-^ 

10^ 

10-^ t 

10= ^ 

-( 

10 

io' 

-1- 1- 

10 6~ 

10"^ I 
10" 
10" 



7c - \SNfj = 14 TeV - p + p Collisions 



"T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 

O DPMJET J. Ranftetal. 
NLO pQCD : CTEQ6M, KKP 



M=(i=Mp= 1 X 



M=(i=Mp : from 0.5 xppto2.0xpp 



_l I L. 



20 



_l I L. 



40 



_l I L. 



60 



_l I L. 



_J > > ^ 

80 100 

Pj (GeV/c) 



Fig. 29: Same as Fig.|28|but at = 14 TeV. 
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while for the standard (CTEQ6M, KPP, all scales equal to Pj,) NLO predictions the behavior is slightly 
harder: 
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Fig. 30: The pion transverse momentum spectrum for deuterium-gold collisions at ^/s = 200 GeV. The solid line represents 
the standard NLO predictions using the EKS shadowing parametrization. The theoretical predictions are evaluated by averaging 
the invariant cross sections in an interval ±0.35 unit of rapidity around rapidity 0. 



We consider now deuterium-gold scattering at 200 GeV and proton-lead scattering at ^/s = 
8.8 TeV and we display, in figs. j^HlandjlT] the DPMJET results together with the NLO QCD predic- 
tions, using the shadowing model of |39j||401. At RHIC the agreement between data and the standard 
NLO predictions is good although the latter may not be steep enough and tend to fall above the data for 
> 7 Gev/c. The magnitude of uncertainties of the NLO calculations is the same as for pp collisions. 
The Dpmjet predictions have the right order of magnitude at low p^ but they tend to decrease somewhat 
too fast as p^ increases. Work is in progress to understand this point l85l . 

For LHC in pPb collisions the two models are compatible within the error bars but the slope of 
the cross section is slightly steeper for DPMJET. As before the grey band shows the range of predictions 
for our standard inputs while the solid line is obtained using MRS99 and BKK functions with all scales 
equal to Pj. . The latter result is equivalent to the standard one with all scales equal to 2p^ . 

Turning now to A-i-A collisions, the comparison of the vr*^ spectrum in gold-gold collisions at = 
200 GeV together with the theoretical estimates is shown in Fig. {32] The agreement between theory and 
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Fig. 31: The pion transverse momentum spectrum at ?7 = for proton-lead collisions at ^/s — 8.8 TeV. The grey band 
represents the standard NLO predictions as in previous figures. The solid line is the NLO QCD result using MRS99 1 82 1 and 
BKK L31J parametrizations with all scales equal to . All NLO predictions include EKS L39..40J shadowing. 



experiment is surprisingly good over the whole range. In the DPM approach no jet quenching and 
no shadowing are introduced but a string fusion mechanism is at work which is calibrated to reproduce 
the particle multiplicity. It also reproduces the minimum bias rapidity distributions. In the NLO QCD 
approach the theory is consistent with the data at large p^, provided a jet energy loss mechanism is 
introduced. In the figure we display the result with model 3) of sec. 15.131 i.e. Bethe-Heitler energy 
loss with e 5^ = 0.05E'" for gluons and e = 0.025-E'^ for quarks and with \a = 0.5 for gluons and 1 
for quarks f^T W) (see Eq. (O and the discussion following it). Below = b GeV/c the NLO QCD 
predictions fall below the data: this is the region where thermal production of vr may become relevant 
as will be discussed in Chapters l6l and 171 

The extrapolation to lead-lead collisions at y/s = 5.5 TeV are compatible (see Fig. l33l even 
though the treatment of final state effects in both models is completely different. The DPMJET results 
are slightly steeper and one finds, in the range 3 < Pj, [GeV/c] < 20, 



dp^dy 



for NLO QCD 



,6.24 



, for DPMJET , 



(14) 
(15) 



Atpj, =5 GeV/c one expects, integrated over azrmuthal angle and per unit of rapidity, 160 events/GeV ^/sec. 



5.5 Phenomenology of prompt and decay photon production 

In this section we discuss the photon spectrum and compare the rates of production of "prompt" and 
"decay" photons. While for tt^ production we have two independent models at our disposal, for the 
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Fig. 32: Comparison between PHENIX data on gold-gold collisions at -y/s = 200 GeV and theoretical predictions. The 
theoretical predictions are evaluated by averaging the invariant cross sections in an interval ±0.35 unit of rapidity around 
rapidity 0. 



prompt photon spectrum one can, at the moment, only obtain NLO QCD predictions . For decay photons 
the Dpmjet code exphcitely includes the radiative decay of resonances (vr ^, t] and other hadronic reso- 
nances of the low mass meson and baryons multiplets) so that one will be able to assess the importance 
of various channels to the decay photon spectrum. One can also obtain a prediction for the full decay 
photon spectrum from NLO QCD via standard convolution formulae based on the vr ^ spectrum. Assum- 
ing a power behavior for the vr'^ spectrum, da'^ /dp^, ~ p^", one can derive the following formula for 
the ratio of the decay photon spectrum over the vr*^ spectrum^^ |86J : 

2 

No isolation cuts have been applied when calculating the photon rates. 

The results are displayed in figs. |34]to|^ In all figures, the solid Une indicates the "standard 
NLO QCD" prompt photon predictions based on CTEQ5M or CTEQ6M for the structure functions and 
BFG set II [35 1 for the fragmentation functions of partons into a photon. Using (the less favored but 
still compatible with present data) BFG set I |35| would lower the predictions. The grey band, when 
shown, indicates the range of uncertainties in the NLO calculations: the lowest prediction is obtained 
with BFG set I and all scales set to 2 p^. while the highest one is obtained with BFG set II and all scales 
set to 0.5 p^. For the decay photon spectra, we show as a dashed line the estimates of NLO QCD using 
a generalization of Eq. (fT6t |86|. Concerning the Dpmjet predictions, the open circles indicate the full 
photon decay spectrum, the full circles the spectrum with the tt" decay contribution removed and the full 
squares the spectrum with both vr'^ and rj decays removed. 

At y/s = 200 GeV, one sees that the decay photon spectrum dominates the prompt photon one 
and, furthermore, one notices the excellent agreement between the Dpmjet and the NLO QCD estimates 
which is not surprising since both models agreed very well in their vr*^ spectrum predictions. The range 

^^Slightly more elaborate fits to the tt" cross sections are used to derive the decay photon spectra in the NLO QCD estimates. 
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Fig. 33: Same as Fig.|28|for lead-lead collisions at y/s — 5.5 TeV. 
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Fig. 34: Comparison of "prompt" and "decay" spectra in pp collisions at 200 GeV. CTEQ5 [80] structure functions are used for 
the NLO QCD predictions. For the prompt photon spectrum the BFG 1 35 1 parametrization of fragmentation functions is used 
while the KKP 1321 fragmentation functions are used for the NLO QCD estimates of decay photons. The theoretical predictions 
are evaluated by averaging the invariant cross sections in an interval ±0.35 unit of rapidity around rapidity 0. 
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Fig. 35: Comparison of the theoretical prompt photon spectrum in pp collisions at 200 GeV, shown in the above figure, with the 
preliminary PHENIX data (full dots). The data points are collected from a figure shown in ref. 1871 . The theoretical predictions 
of decay photons are those of Fig. 1341 
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Fig. 36: Comparison of "prompt" and "decay" spectra in pp collisions at ^/s —5.5 TeV. CTEQ6 1 81 1 structure functions are 
used for the NLO QCD predictions. For the prompt photon spectrum the BFG 1 35 1 parametrization of fragmentation functions 
is used while the KKP 1 32 1 fragmentation functions are used for the NLO QCD estimates of decay photons. The cross section 
is evaluated at rapidity. 
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Fig. 37: Same as Fig. l36l at y/s — 14 TeV. 
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Fig. 38: NLO QCD predictions for tlie ratio 7prompt/7r° in proton-proton, proton-nucleus and nucleus-nucleus collisions at 
200 GeV and 5.5 TeV. The Au+Au ratio at = 200 GeV and the Pb+Pb ratio at ^ = 5.5 TeV are modified by thermal 
effects at the lower end of the spectrum and are predicted to be higher than shown (see sec. l7.2t . 



of uncertainties, due to changes of scales and of photon fragmentation functions, in the NLO QCD rate 
of prompt photons is less than it 50% at low and much smaller at high p^. Displayed in Fig.l35lis 
the comparison between the theoretical predictions of prompt photon production and the experimental 
results recently obtained by the PHENIX collaboration 1871 : the agreement is very satisfactory. 

At LHC, Figs. |36l and 1^ the same features are found except that the prompt photon rate is rela- 
tively smaller, at the level of the decay photon spectrum with the vr*^ and rj contributions removed. At low 
Pj, , the uncertainty in the prompt photon spectrum is largely due to the choice of photon fragmentation 
functions: for ^/s = 5.5 TeV or 14 TeV, using BFG set I f35l decreases the predictions by a factor of up 
to 2.5 at ~ 3 GeV/c, where the bremsstrahlung component is important, but by 10% or less for 
above 20 GeV/c: this large variation at low is associated to uncertainties in the gluon fragmentation 
into a photon which is hardly contrained by present data. At LHC the uncertainty associated to the CTEQ 
structure functions is much smaller: ±10% atp^ =3 GeV/c and about ±2.5% at p^ = 100 GeV/c sim- 
ilarly to the results for vr*^ production. As for predictions based on the MRS99 parametrization, they are 
reduced by less than 15% compared to our standard results, with the largest reduction occuring at the 
small values of transverse momentum. 

In more details, the results for pp scattering at LHC can be summarized as follows: 1) at the lower 
end of the spectrum, there is almost 10 times more photons from vr'^ decays than from all other sources 
taken together; 2) decays from rj's are again about 10 ten times larger than those from "other" resonances; 
3) the level of production of prompt photon is similar to that of decay photons from "other" hadronic 
resonances. 

At RHIC as well as at LHC, the ratio 7prompt /Tdecay slowly increases with p^. To see this it is 
customary to display the ratio 7prompt/7r'^: this is done in Fig. 138) where the predictions for RHIC and 
LHC are compared. The ratios are rather small, specially at LHC, and isolation cuts will certainly be 
necessary to increase them. 
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Fig. 39: Comparison of "prompt" and "decay" spectra in proton-lead collisions at 8.8 TeV. 



In Fig.|39lthe results for pPh scattering at LHC are summarized: the DPMJET predictions are in 
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Fig. 40: Comparison of the prompt and the decay photon spectra in Au+Au collisions at ^/s — 200 GeV in NLO QCD. The 
Dpmjet results for decay photons are also shown. The theoretical predictions are evaluated by averaging the invariant cross 
sections in an interval ±0.35 unit of rapidity around rapidity 0. 



good agreement with the NLO QCD estimates based on the MRS99 and BKK parametrizations, although 
they tend to have a slightly steeper slope as already noticed in Fig. for pion production. As expected 
the prompt photon signal is much below the decay photon spectrum. 

The ratio 7prompt /tt'^ is shown in Fig. |38]for dAu collisions at 200 GeV and pPb collisions at 
8.8 TeV. At RHIC very little change is seen compared to pp collisions at the same energy. We also 
looked at the effect of assuming an ^4^^" dependence of the dAu cross sections, with a = 1.08 for vr'^ 
production and a = 1.04 for 7 production. It is seen to have a moderate effect on the estimate of the 
ratio. For pPb scattering at LHC a small energy dependence is observed compared to pp scattering at 
5.5 TeV with a larger decrease of the ratio at the larger p^, values. 

We turn finally to the case of nucleus-nucleus collisions. The results for RHIC are shown in 
Fig- EH! One sees that in Au+Au collisions the ratio 7prompt /7dccay increases with p^ and approaches 1 
at Pj, = 10 GeV/c. The value of this ratio is much higher than for pp scattering at the same energy. In 
the NLO QCD model this is due to the "jet quenching" mechanism which decreases the rate of decay 
photons more than that of prompt photons: indeed only the bremsstrahlung component is affected in the 
latter case. An estimate of the relative increase of prompt photons in A+A collisions is also seen by 
plotting the ratio 7prompt / tt'' as done in Fig. |38j this ratio is approximately 3 times higher in Au+Au 
collisions than in pp collisions. 

The same pattern is seen at 5.5 TeV (Fig. ^3 where prompt and decay photons become comparable 
above p^, = 60 GeV/c, in contrast with the pp case as seen in Fig.|36l Looking at the medium p^ range 
one should be able to extract a prompt photon spectrum from the data for p^, > W GeV/c since the ratio 
7prompt/7r'^ becomes of the order of 5% or larger as shown in Fig. |38] As for RHIC this ratio is two 
or three times larger than in pp collisions. In this figure only the non thermal contribution to heavy ion 
collisions is taken into account. Thermal production of pions and photons occurs at low Pj, values in 
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Fig. 41: Comparison of "prompt" and "decay" spectra in lead-lead collisions at 5.5 TeV. 



heavy ion collisions. At LHC it is important below = 10-15 GeV/c and the ratios turn out to be larger 
than shown as will be discussed in sec. 17.21 



5.6 Small mass lepton pairs at large transverse momentum 

As seen in the previous section the prompt photon rate is much below the rate of decay photons at low 



Pj, values. The dominant backgrounds are the decays tt^ — > 77 and rj — > 77. If one considers instead 
of real photon production, the emission of virtual photons (lepton pairs) part of this background can be 
eliminated. For example, considering the production of an electron pair in the mass range Mg+g- = 
[0.2,0.6] GeV/c^ one gets rid of the vr'' background (the lepton pairs from the Dalitz decays of the vr" 
are below this mass range) and one stays below the p,uj,cl) ~^ e+e~ decays. The r] background is also 
reduced because of the mass constraints on the Dalitz pairs. 

From the theoretical point of view, the production rate of small mass lepton pairs at relatively 
large transverse momentum is very similar to that of real photon. To get a rough estimate of the rate 
of production of a Drell-Yan pair, integrated over a given mass range, we follow the procedure of f88| 
where the virtual photon mass is neglected, compared to the other scales p^ and ^/s, everywhere except 
in the bremsstrahlung component: indeed, in the fragmentation process, the virtual photon mass, rather 
than Aqcd, acts as a cut-off of the final state collinear singularity. Following this procedure, we find 
that the rate of production of an electron pair in a finite mass range is given, with a good accuracy, by the 
simple relation 

- — — ~ Ce+e- a - — — (17) 

where Cg+g- ~ 0.3 for 0.2 GeV/c^ < M ^+^- < 0.6 GeV/c^ valid in the range 2 GeV/c < p^ < 

100 GeV/c, and C^+e 0.2 for 1 GeV/c^ < M^+^- < 3 GeV/c^ in the range 4 GeV/c < p^ < 

100 GeV/c. These numbers are independent of the fragmentation functions used. These estimates re- 
fer only to the "prompt" production mechanism. Using, from Appendix II, the luminosity factors for 
ALICE at ^/s = 5.5 TeV, the rate of production of a prompt lepton pair in the mass range 0.2 GeV/c^ 
< Mg+g- < 0.6 GeV/c^ is estimated to be 6 pairs per second with p^, = 3 GeV/c and 0.1 pair per 
second with Pj, = 8 GeV/c. The corresponding numbers in lead-lead collisions are 0.4 and 0.02 respec- 
tively. We should say at this point, that thermal production of dileptons (see next chapter) will increase 
the direct production of lepton pairs. 

The production mechanism we just discussed suffers from a background with can be estimated 
using Dpmjet: indeed this code contains the Dalitz decays of resonances as well as the semi-leptonic 
decay processes of charmed and heavy flavor resonances. We show in Fig.|42lthe prompt lepton pair 
spectrum for 0.2 Gev/c^ < Mg+g- < 0.6 Gev/c^ (Eq. (fTTb ) together with the background from lepton 
pairs originating from the same resonances (correlated pairs). One sees that, as expected, at high enough 
Pj, the signal becomes comparable to the background. 

Unfortunately, this estimate of the background is incomplete as one expects a huge contribution 
from uncorrelated pairs due to the large combinatoric factor when collecting pairs of leptons: in a typical 
Dpmjet event one expects about 180 e~ and including vr ^ decays but only 17 e~ and excluding 
vr "^'s. This is shown in the figure by the open dots (including vr ^ decays) and full dots (exluding vr ^ 
decays): obviously the situation is not as favorable since the background is now two orders of magnitude 
above the estimated signal! Note however the nice feature that the combinatoric background from vr "^'s 
disappears for Pj. > 2 GeV/c. One way to tame the huge background may be to subtract from the data 
the spectrum of like charge pairs in the same kinematic range or pairs constructed from different events. 
Using low statistics runs of Dpmjet one finds that subtracting the like-charge pair spectrum from the 
unlike-charge pair spectrum one gets numbers which are very close to the correlated pair background 
whenever the transverse momentum of the pair is larger than 1.5 GeV/c l89l . The statistical fluctuations 
are large and clearly further studies are necessary to remove the uncorrelated background. 
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Fig. 42: Lepton pair spectrum at = in Pb+Pb collisions as a function of Pj, for y/s = 5.5 TeV. The solid line indicates the 
prompt (signal) NLO QCD predictions while the squares and the points represent backgrounds of different origins calculated 
using Dpmjet. Thermal production will increase the signal (see sec. l7.3L 



6. PHOTONS FROM THERMAL, EXPANDING FIREBALL 

P. Aurenche, F. Gelis, G. Moore, H. Niemi, R. Rapp, K. Redlich, P.V. Ruuskanen, S.S. Rdsdnen, D.K. Sri- 



vastava 



For the calculation of thermal photons — photons from secondary interactions among particles in 
the expanding fireball of matter produced in a heavy ion collision — the emission rates from hot matter 
and the space-time evolution of the matter are needed. The equilibrium emission rates in both QGP and 
the Hadron Resonance Gas (HRG) are reviewed below. The space-time evolution of expanding matter is 
described in terms of relativistic hydrodynamics with the assumption of longitudinal boost invariance. 

There are uncertainties both in the emission rates and the description of the nuclear fireball. One 
such uncertainty affecting both the emission and the evolution of produced matter arises from the fact 
that the initially produced quark-gluon matter is not expected to be in chemical equilibrium. E.g., in 
the pQCD -I- saturation model |90| described below, most of produced partons are gluons. Assuming 
kinetic equilibrium and describing the deviation of number densities from equilibrium values in terms 
of fugacities, Aj, we typically obtain Xg ~ 0.6...1 and Ag ^ 0.2 in the initially produced parton system. 
The quark fugacity is assumed to be the same at LHC energies as the antiquark fugacity. The effect of 
non-equilibrium values of fugacities on the photon emission rates is studied below. Using rate equations 
for reactions which change the parton numbers |91jj92jj, we will study in detail also the time evolution 
of fugacities. 

The description of the space-time evolution of a fireball produced in a nucleus-nucleus collision 
must begin with characterization of the initial state of produced matter. Information on initial state is ob- 
tained from calculations of particle production from primary interactions in nuclear collisions. There are 
several model calculations, some more detailed than the others, of hadron spectra in heavy ion collisions 



at RHIC energies. We use these calculations or, when necessary, their extensions at the LHC energy 
as the first step in obtaining the initial state for the space-time evolution of produced matter. The most 
important quantities in determining the initial conditions are the total hadron multiplicity, dN/ dy, which 
fixes the the total initial entropy and the production time scale, tq. The initial densities are proportional to 
dN/dy and I/tq. The models for particle production differ also in the transverse dependence of nuclear 
densities. 

After reviewing the emission rates both in the quark-gluon plasma and in the hot hadronic gas, 
we consider in Sec. 16. 51 different models for particle production in order to establish a reasonable range 
of uncertainties in the initial state. The hydrodynamic evolution is dealt with in the following section, 
while the predictions for the rates of thermal pions and photons are discussed in Sees. 16.71 to 16.91 The 
last section is devoted to estimating the yield of small mass lepton pairs in the quark-gluon plasma. 



6.1 Thermal photon emission rates 

The rate of production, per unit time and volume, of a real photon of momentum {E, p) in a system in 
thermal equilibrium is calculated via the formula [t93ll94i : 



where li'^ ^^[E, p) is the retarded photon polarization tensor. The pre-factor 1 / (exp(£'/T) — 1) provides 
the expected exponential damping when E ^ T. 

A loop expansion of 11^ is constructed with effective propagators and vertices. For production of 
photons in the quark-gluon plasma the basic degrees of freedom are quarks and gluons and the corre- 
sponding physics will be discussed in the next two sections, while for production in hot hadronic matter 
the relevant degrees of freedom are mesons (vr, p, ai, ■ ■ ■) and baryons, and this will be the subject of 
Sec. El 



6.2 Photon emission rates from the chemically equilibrated quark-gluon plasma 

P. Aurenche, F. Gelis, G. Moore 

The most complete calculations are those based on the improved perturbation theory of Braaten 
and Pisarski, i.e. on the Hard Thermal Loop (HTL) resummed effective theory l95l l96l . In the HTL 
approach, the quarks and the gluons acquire an effective mass of the order of y/a^T, where ^/a^ is the 
strong coupling, due to their interactions in the plasma. Quark and gluon exchange mechanisms on a 
long distance (equivalently small momentum transfers) are also modified and they become effectively 
screened by mass effects of ©(-^/aJT). Thus, thermal effects tend to regularize the infra-red behavior of 
the theory. Estimates of Eq. ( fT^ . in the HTL effective theory, based on 1-loop, 2-loop and multi-loop 
contributions have been performed: each case corresponds to different physical processes. We review 
them in turn. 



6.21 Real photons 

At one-loop order in the HTL theory, the production of a hard (E » T) real photon is given by the 
diagram in Fig. ^\ One of the fermion leg in the loop may have a soft momentum flowing though it, 
the other one beeing necessarily hard. Following the HTL approach, thermal corrections are resummed 
on the soft line and an effective fermion propagator should be used. When taking the discontinuity of 
the diagram one of the cut fermions is necessarily space-Uke (denoted Q in the figure) in the case of 
real photon production. Cutting through the effective quark propagator of momentum Q exhibits Landau 
damping i.e., the exchange of a virtual quark in the scattering of quarks and gluons in the medium. 
In this way one obtains the processes contributing to photon production, namely Compton scattering 



and qq annihilation into a photon and a gluon. As a result of resummation, the potential singularity at 
momentum transfer in the fermion propagator Q is screened by an effective mass rriq. Taking into 
account this thermal correction to the quark propagator, the imaginary part of the photon polarization 
tensor can be calculated (971 HH. For hard photons, it reads: 



Imn''/(S,p) 



In 



(19) 



Note that the mass rriq is given by m'^ = ttq^ CfT^ with Cf = {Nl - l)/2Nc. The numerical factor 



(E,p) 
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otmaoJ — ^ — 



Fig. 43: Left diagram: one-loop contribution to liard real photon production; right diagrams: examples of diagrams contributing 
to photon production, annihilation and Compton scattering. The symbol • on the quark propagator indicates that the effective 
propagator is used. 

5/9 is the sum of the quark electric charges squared for 2 flavors (u and d); for 3 flavors (u, d and s), this 
factor should be replaced by 6/9. It is clear from the expression above that the quark thermal mass acts 
as the cut-off of a logarithmic collinear singularity. 

For some time this was thought to be the final answer for the photon rates at 0(0;^). It then became 
clear that some formally higher order processes (see Fig. l44b are in fact strongly enhanced by collinear 
singularities. The corresponding processes are bremsstrahlung emission and off-shell annihilation. A 
common property of these two diagrams is that they have an off-shell quark next to the vertex where the 
photon is emitted, and the virtuality of this quark becomes very small if the photon is emitted forward. 




Fig. 44: Left most diagrams: two-loop contributions to hard real photon production; right diagrams: bremsstrahlung and 
off-shell annihilation processes. The symbol • on the gluon propagator indicates that the effective propagator is used. 

Again, the quark thermal mass mq prevents these diagrams from being truly singular. However, contrary 
to the one loop diagrams, the singularity is linear instead of logarithmic, and it brings a factor T ^/m-q. 
Combined with a1 from the vertices, these diagrams turn out to be also of order 0{as) I99l ll00lll01l 
II02I . For 3 colors and 2 light quark flavors, the 0{as) contribution of these two diagrams is exactly: 



^ „-R ,, , ^ , 32 5aas 

Imn ,/CE,p = 

^'' ^ 37r 9 



T3 

E 



(20) 



In this formula, the term inl/E comes from the bremsstrahlung diagram and dominates for soft photons, 
while the term linear in E comes from off-shell annihilation which therefore dominates for very hard 
photons {E > T).2^ 

Given the enhancement in the diagrams of Fig. |S1 one may wonder if higher order diagrams also 
contribute to the same 0{as) order. To discuss the issue in physical terms, it is convenient to define the 
concept of photon formation time. Consider a virtual quark of momentum R = P + Q which splits into 
an on-shell quark of momentum Q and a photon of momentum P (see the processes in Fig. l44l . The 
photon formation time can be identified with the lifetime of the virtual quark, which is itself related to 
its virtuality by the uncertainty principle. A simple calculation gives: 

t-^ ^ 5E = ro - ^r^ + ml^ [ql + m^] , (21) 

where the 3 -momentum of the photon defines the longitudinal axis. The coUinear enhancement in the 
diagrams of Fig. |S1 due to the small virtuality of the quark that emits the photon, can be rephrased by 
saying that it is due to a large photon formation time, of 0{l/asT). If the formation time is large, the 
quark can rescatter in the medium, while emitting the photon, as it is shown in Fig. |45] This can be 
partially taken into account by introducing a coUisional width T ~ asTln{l/as) on the quarks in the 
calculation of the diagrams of Fig.|44]| 103 1. A large sensitivity to this parameter was found at leading 
order, thereby indicating that an infinite series of diagrams must be resummed in order to fully determine 
the 0{as) photon rate. This phenomenon is nothing but a manifestation of the Landau Pomeranchuk 
Migdal (LPM) effect ll()41II()5lfT()6l . 




Fig. 45: A ladder correction to bremsstraiilung. 



A considerable progress was made recently in 1, 107.1 : it was shown that there are infrared cancel- 
lations between diagrams of different topologies, and that these cancellations remove any sensitivity to 
the magnetic scale. Physically, this cancellation can be interpreted as the fact that ultrasoft scatterings 
are not efficient in order to induce the production of a photon. As a consequence, only the ladder family 
of diagrams needs to be resummed in order to obtain the complete leading 0{as) photon rate. The re- 
summation of this series of diagram can then be performed in two steps summarized in Fig. |46l The first 

Fig. 46: Resummation of ladder diagrams. 



one is a Dyson equation for the photon polarization tensor, whose explicit form is I107II1081I109I : 

r+oo 2,2 r j2 

Imn /(i?,p) ^ aN,j^ dqo K(ro) - n,{qo)] Re J ^ ■ f{q^) , (22) 



'For 3 flavors, tiie same formula holds but the numerical pref actor is replaced by a very complicated expression 11021 . 



with ro = E ->rqQ, ng{rQ) = l/(exp(ro/r) + 1) the Fermi-Dirac statistical weight, and where the dimen- 
sionless function f{q±) denotes the resummed vertex between the quark line and the transverse modes 
of the photon (this is represented by the shaded vertex in the above pictures). In the Dyson equation, 
this function is dotted into a bare vertex, which is proportional to Qj^ (the photon has a transverse polar- 
ization and therefore the coupling vanishes for a coUinear emission by a quark). The second equation, 
that determines the value of f{q±), is a Bethe-Salpeter equation that resums all the ladder corrections 
111071 11081 IT09l : 

^f{q±)=2q^+47TasCfTj ^C{l^)[f{q^+l^)-f{q^)], (23) 

where tp is the time defined in Eq. and where the collision kernel has the following expression: 
C{1±) = "idebye/'-L('_L + "^debye) (liSI- Notc that in the Dyson equation, the quark propagators are 
dressed (by resumming self-energy corrections which is akin to introducing the collisional width or 
damping rate) to match the resummation performed for the vertex, so that gauge invariance is preserved. 
It is this dressing on the quark propagators which is responsible for the term —f{q±) under the integral 
in Eq. (l23l . From this integral equation, it is easy to see that each extra rung in the ladder contributes a 
correction of order asTq^rQ/ Em^, in which the drops out. Therefore, all these corrections contribute 
to 0{as) to the photon rate. Note again that the only parameters of the QGP that enter this equation 
are the quark thermal mass mq and the Debye screening mass, mdobyc- This integral equation was 
solved numerically in |108|, and the results are displayed in Fig. In this plot, 'LPM' denotes the 

a =03, 3 colors, 2 flavors, T=l GeV 




E (GeV) 

Fig. 47: 0{aa) contributions to the photon production rate in a QGP. The parameters used in this plot are Us — 0.3, 3 colors, 
2 flavors and T = 1 GeV. 

contribution of all the multiple scattering diagrams, while '2 — > 2' denotes the processes of figure 143) 
The single scattering diagrams (Fig.l44l are also given so that one can appreciate the suppression due to 
the LPM effect. 

6.22 Small mass lepton pairs 

For the reasons given in Sec. l5.6l it is also interesting to consider the production of a small mass leptons 
pair at hard momentum. The calculation parallels that described above with added technical complica- 
tions due to the mass of the virtual photon (lepton pair). The relation between the rate of production of a 



pair of mass M at 4-momentum {E,-p) and the virtual photon polarization tensor is 



dtd^^dEd^p ~ 127r4 M2 eME/T) - 1 m''(^'P) • (24) 

Neglecting the mass of the pair compared to its energy, the imaginary part of the photon polarization ten- 
sor at one-loop is still given by Eqs. ( fl^ but the collinear cut-off (quark thermal mass) in the logarithmic 
term now involves the photon mass IllOllllll . The two-loop 11121 and the multi-loop 11131 calculations 
have been carried out. In the latter case, it is convenient to separate the contributions of the transverse 
and longitudinal polarizations. For the transverse case, one obtains again Eqs. d23l . (l22l while a related 
integral equation has been derived and solved for the case of the longitudinal polarization 11131 . 



6.3 Hard photon production from hot out-of-chemical equilibrium QGP 

F. Gelis, K. Redlich 

As a first step we consider only the basic reactions for photon production due to Compton scatter- 
ing and annihilation of quarks (see Fig. l43b . We assume that off-equilibrium effects on the momentum 
of QGP constituents can be parameterized by a modification of their distributions. A HTL-resummation 
is also supposed to be valid to screen soft and collinear (mass) singularities. 

To account for non-equiUbrium effects in partonic medium we follow the approximations de- 
scribed in detail in I114llll5llll6ll . It amounts simply to replace the thermal particle energy distribution 
n{k q) by their non-equilibrium counterparts, i.e. in general by Wigner distributions n{kQ,X). In 
terms of the cumulant expansion of the initial density operator, this amounts to approximate the initial 
correlations by the second cumulant only, i.e. to neglect all non-gaussian correlations. As a further 
approximation, we ignore for the time being the possible dependence on the center-of-mass coordinate 
X, essentially assuming a homogeneous and isotropic medium^^. Practically, we take the following 
distributions for quarks and gluons: 

^«(^o) - I ^ _ ^^(i^^i)) n,{ko} - I _(i + for ^ Q (25) 

with the Jiittner parameterizations 

^g(g)(l^ol) = e|fcol/T^A ' ^^^^ 

expressed by introducing the fugacity parameter A fTW\, which is assumed to be energy independent. 
Obviously A 7^ 1 in the case of chemical non-equilibrium. 

Neglecting the slow macroscopic scale X compared to the fast microscopic scale, the structure 
of the resulting perturbation theory is very similar to the equilibrium Closed Time Path formalism. For 
instance, the 21 components of the non-equilibrium bare propagators for bosonic and fermionic fields, 

iD2i{P) = 2^e(po)5(P')(l + rigipo)) , iS2i{P) = 27r6{P^)f'e{po){l - n,{po)) , (27) 

have formally the same structure as the corresponding ones in equilibrium, but depend now on the mod- 
ified distribution functions (Eq.|55Jl- However, for a dressed propagator the above is not in general valid, 
as the appearance of some additional out-of-equilibrium terms is to be expected I110llll9lll20l . There- 
fore, in order to calculate a physical rate with this formalism, one needs to account not only for the 
modifications of the particle momentum distributions but also for modifications of their propagators. 

^''This is a valid assumption if there is a scale separation between the macroscopic scale (i.e. the scale of hydrodynamical 
inhomogeneities, encoded in the variable X) and the microscopic scale corresponding to the process under study 11171 . The X 
dependence is reintroduced later by making the fugacities X dependent. 



The rate of real photon emission, already given in Eq. (fTSl in terms of the retarded photon polar- 
ization tensor, can also be written in terms of its component 11 12 in the Closed Time Path formalism: 



E 



dN 



dtd?x(Pp 2(27r)' 



(28) 



As in the equilibrium case the simplest diagrams, in the loop expansion of 11 12, correspond to photon 
production by annihilation and Compton processes 9 + ^^5 + 7, q{q) + g ^ q{q) + 7 , as shown 
in Fig.^l In the off-equilibrium situation one needs first to determine the appropriate generalization of 
the effective quark propagator. In the un-resummed theory it is known that ill-defined terms of the form 
5{Q'^)/Q^ are generated. However, in the resummed approach one can show that quark propagator can 
be approximated by a form similar to the equilibrium one L12l1ll22l : 



-nq{qo) Re 



S„ iq)+i£qo 



(29) 



where S is the quark HTL self-energy, evaluated with non-equilibrium distribution functions. The 
only non-equilibrium effect on this HTL self-energy amounts to a modification of the quark thermal mass 
niq (the pref actor m L ): 



vr- 



;Cf / EdE[ng{E) + nq{E)] = nasCfT 
Jo 



1(^.4 



(30) 



which reduces to the expression given in the previous section in the equilibrium case (A^ 



1). 



The hard photon production rate from a system away from equilibrium can then be calculated 
using well established methods in equilibrium field theory. The final result in the leading order in Us 
is L122J (c/Eq.[^for the equilibrium result): 
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where 7 is Euler's constant, C,' the derivative of Riemann's function and where we have defined: 

i-T 



(32) 



= 1)^0.67, C- 



n=2 



00 

E 

n=2 



ln(n - 1) ~ -0.04 . 



(33) 



The above result shows that, as in the equilibrium medium, the generalized thermal mass provides a 
self consistent cut-off for the logarithmic singularity. The dynamical screening of the mass singularity 
seen in Eq. (13 It (here given for the distributions Eq. d26l ) actually does not depend on the explicit form 
of the non-equilibrium quark and gluon distribution functions. Changing the parameterization of these 
functions enters only through the redefinition of the mass parameter in Eq. dSTT i and the constant C in 



Eq. (I32t . keeping at the same time the functional form of the rate unchanged. For Aq 
(I32I reproduce the well known one-loop results of the previous section. 



An 



1 Eqs. dlD, 



The out-of-equilibrium approach discussed above for the case of Compton and qq annihilation 
processes can also be extended to the more complicated processes of Fig. |44] The main effect of the 
fugacities comes from the fact that the number of emitters is reduced. To take this into account, one just 
need to rewrite the statistical functions in Eqs. ( fTSl and d22t as: 

exp(g/r) - 1 ^^^'-^^ + "'g(Po)] = -nqipo + E)nq{-po) . (34) 

Since the equilibrium quark distribution function satisfies nq{—x) = 1 — nq{x), one sees readily that 
the right hand side of the previous formula gives the combination of distribution functions one expects 
from kinetic theory in each of the kinematical domains: for > (bremsstrahlung of a quark) we 
get nq{pQ + E){\ — nq{po)), for —E < po < (quark-antiquark annihilation with a scattering) we get 
n'q{po+E)nq{\po\), and for Po < —E (bremsstrahlung of an antiquark) we getnq{\po\){l—nq{\pQ+E\)). 
Given its interpretation in terms of kinetic theory, the right hand side of Eq. d34t can be generalized^^ to 
out of chemical equilibrium situations by simply replacing the Fermi distribution by the form of Eq. (126 1 . 
For this to make sense, the out-of-equilibrium distribution at negative arguments must be defined in such 
a way that nq{—x) = 1 — nq{x) holds, which is indeed the case (see Eqs. d25l '). 

The other places where the fugacities enter in the calculation of the multiple scattering diagrams is 

and via the Debye mass ni^^^^^. 



via the quark thermal mass -rni and via the Debye mass ni^,^. We recall here the expression for these 



masses: 



"iq = E (Ag + ^j, (35) 



mL^ye = ^iNcXg + Nf^). (36) 



It is interesting to remark that using fugacities Xg, Xg < 1 will obviously reduce the photon rate, 
because the number of emitters is reduced, however, since the thermal quark mass is also decreased, the 
coUinear enhancement at work in bremsstrahlung and off-shell annihilation processes is increased. The 
two effects partly compensate 1 123il 



6.4 Thermal photon production from hadronic matter 

R. Rapp 

The objective in the following is to assess the emission rate of photons in a hot hadron gas as ex- 
pected to be formed under conditions characteristic for LHC heavy-ion collisions. The rate is calculated 
acording to Eq. (fTSb . 



6.41 vrpai gas and form factors 

Early calculations Wfl of photon radiation from hadronic matter were guided by the analogue to lowest- 
order pQCD processes in a QGP, i.e., Compton and annihilation graphs involving the most abundant 
constituents in a meson gas. These were vrvr p^, irp ttj as well as /? — > ttttj (and to ir^), 
which are based on the vrvr/o-vertex with a coupling constant determined from the free p (u) decay. At 
photon energies in excess of ~ 1 GeV, t-channel pion exchange in the second reaction completely 
dominates the rate. Shortly thereafter 11241 . the ai s-channel pole graph was recognized as a poten- 
tially relevant contribution, due its large coupling to irp states. At sufficiently high energies, s-channel 
resonance graphs are, however, suppressed due to — mj^) powers in the intermediate propagators 
(mji: resonance mass), rendering ai t-channel graphs increasingly important. A systematic evaluation 
of photon production from a irpai gas, based on the Massive Yang-Mills (MYM) Lagrangian approach, 

^^But it would not be correct to use the out-of-equilibrium in the left hand side [nq(po + E) — rig (po)]/(cxp(_E/r) — 1). 
In fact, the equality of Eq. <34t is only valid with equilibrium distributions. 



has been performed in ref. 1 125 1. Here, chiral symmetry provides relations between couplings, leaving 
4 free parameters (including non-minimal coupling terms) allowing for a reasonable phenomenology for 
the free irpai system with electromagnetic decays. Two different sets implying different off-shell behav- 
iors lead to total rates that differ by no more than 50%. The dominant contribution to the rate at photon 
energies beyond ~ 0.5 GeV was confirmed to be due to t-channel (vr and ai) exchanges in the np ttj 
reaction 1 125 1. A convenient parametrization of the pertinent results, which will be used below, has been 
provided in ref. 11261 . 




Fig. 48: Thermal photon production rates from the np tt^ reaction within SU (2) chiral lagrangians including vector mesons 
using the Massive Yang-Mills 11251 (dotted line) and Hidden-Local Symmetry 11271 (full line). The dashed line is from the 
earlier work of ref. 1971 which does not include ai degrees of freedom. 



An important element in applying effective hadronic models at moderate and high momentum 
transfers is the use of vertex form factors to simulate finite hadronic-size effects, which are not accounted 
for in the irpai calcuations of Ref. 1 125 1. For similar reactions their effect has been studied in Ref. [9'T\ 
and found to give a typical net suppression over the bare graphs by a significant factor ~3 at photon 
energies £' ~ 2.5 GeV. We perform a rough estimate of their effect in the present context as follows. 
Consider pion-exchange assuming a standard dipole form factor, 

the average four-momentum transfer in f-channel pion exchange can be approximated according to 

1 1 f 1 1 

7 7^ — 9 , r^9 ! (38) 



(m2-t)2 4£;2 J (m2-t)2 mlAE'^ 







that is 



-tc^2Em^ (39) 

at sufficiently large E > rriT^. Upon multiplying the rate parametrizations of the vr/jai gas 11261 (dotted 
line in Fig.||9ll by F{tf (with A = 1 GeV 1 128 1) we obtain the dashed curve in Fig.|49l The reduction 
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Fig. 49: Effect of hadronic form factors on the photon production rates from a upai gas in the MYM approach of ref. 11251 . 

of the rate in the 2-3 GeV region amounts to a factor of 3-4, quite in line with the exemplary microscopic 
calculation of Ref. |97J. 

We believe that the calculation including formfactors should provide a reasonable baseline for the 
SU{2) case. 

6.42 Additional mesonic sources 

An obvious contribution to photon production not accounted for in the previous section stems from 
radiative decays of heavier mesonic resonances which can not be easily treated in a chiral framework. In 
the many-body approach of Ref s. II128I I1291 the electromagnetic correlator has been evaluated including 
all mesonic resonances up to about 1.4 GeV in mass which exhibit significant coupling to final states 
including rho-mesons (and thus, employing VDM, photons). These are u;(782), /ii(1170), /i(1285), 
7r(1300), 0(1320), i^*(892) and Ki(1270). Pertinent photon rates are displayed in Fig.|50l In particular, 
the tj ^ 7r7 decay generates a large low-energy strength, consistent with the early results of Ref. |97|. 
Note that all hadronic vertices carry (dipole) formfactors with typical cutoff parameters of around 1 GeV, 
as extracted from an optimal fit to measured hadronic and radiative branching ratios within VDM (see 
Ref. fT28i for details). Contributions from still higher-mass mesonic resonances (such as C(;(1420) and 
u;(1650)) have recently been computed 1 13011 and turned out to be negligible. 

Another significant photon source can be expected from additional t-channel exchanges, in par- 
ticular u;(782) exchange in vrp tt^, and pion-exchange in ttK —>■ 7K*(892), 7ri^*(892) —>■ K'j. 
Results for the former, which do not seem to be available in the literature yet, have been obtained re- 
cently I130L indicating an emission strength that is small at low energies but exceeds the in-medium 
many -body results of ref. 11281 beyond E ^ 1.5 — 2 GeV. For the kaon-induced processes, one can use 
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Fig. 50: Various sources of tliermal photon production from a hot meson gas. 

SU{3) symmetry to estimate the relevant couplings from the analogue SU{2) processes, i.e., vrvr —>■ 7/j, 
vrp — > pj, and account for the smaller kaon (K*) densities {nj^j^p^/riT^ ~ 40% at maximal SPS energies 
and beyond). One finds that the strange analogues make up less than 10% in the former case, and about 
20% in the latter. 

6.43 Baryonic sources 

Finally we have to address the role of baryons. Naively, one might not expect significant contributions at 
collider energies where the wef baryon densities are small. However, since the photon is a CP eigenstate, 
the relevant quantity for its production is the sum of baryon and anti-baryon densities, which is not small 
at temperatures close to Tc even in a net bary on-free environment [ 131 1. Moreover, since the correspond- 
ing anti -/baryon content survives the subsequent hadronic evolution in a heavy-ion collision [.132iil33J . 
their effects on photon and dilepton observables remain appreciable also in the later stages. 

Much like in the mesonic case, baryon-resonance decays ought to be rather suppressed at ener- 
gies beyond 1 GeV. No explicit computation of i-channel exchanges has been reported so far. However, 
within the chiral reduction formalism, contributions from nucleons have been assessed in ref. [134|. 
Within a hadronic many-body framework |135 129 128 1, primarily constructed for low-mass dilepton 
applications, the e.m. correlator has been evaluated incorporating a rather extensive set of baryonic pro- 
cesses. These results can be straightforwardly extrapolated to the photon point, and the corresponding 
rate is shown by the solid line in Fig. |^ More specifically, both the mesonic resonance contributions 
as discussed above (short-dashed line), as well as baryonic resonance decays (most notably A(1232), 
A^(1520)) and pion-exchange processes involving nucleons and deltas, are included coherently. The 
baryon-induced emission can be considered as rather reliable for baryon densities up to at least normal 
nuclear matter density, pQ = 0.16 fm^'^, being constrained by photoabsorption spectra on nucleons and 
nuclei [136 1 (at zero temperature). At comparable baryonic densities this approach yields about a factor 
of two more photons than what has been found in the chiral reduction approach 11341 when coupling 
on-shell pions and nucleons (see Ref. 11371 for an update including A(1232) and A^(1520) resonances). 
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Fig. 51: Thermal photon production rate (under conditions resembling LHC energies) for a npai gas 11251 including form- 
factors (long-dashed lines), and within the many-body spectral function approach of ref. 1 129 1 (full line). The total emissivity 
corresponds to the sum of the two rates; uj t-channel and strange-meson processes, not included in the figure, are expected to 
induce an additional increase of about 30-40% at energies above 1.5 GeV. 

Under conditions relevant for LHC energies (as chosen in the figure), processes involving (anti-) 
baryons are the prevailing source of thermal photons at energies E < 1 GeV (similar to what has been 
found for low-mass dileptons [ 131 1). Beyond, the t-channel processes within the light meson gas sector 
take over, dominating the rate at higher energies. 

It is thus suggestive to regard the combination of the vrpai emission rate (including form factor 
effects), supplemented with lo t-channel and kaon-induced proceses, with the hadronic many-body results 
(ai contribution removed) as a semi-realistic approximation to the full hadron gas emissivity. 

6.5 Particle production from primary interactions in nuclear collisions 

K.J. Eskola, J. Ranft, P.V. Ruuskanen 

Results from model calculations of initial particle production are the average numbers, transverse 
energies and momentum spectra of produced gluons and quarks. How to obtain from the momentum 
space densities the space-time densities needed for the hydrodynamic calculation is explained in detail 
in Sec. l6.62l using the pQCD -i- saturation calculation of (mini)jet production |90| as an example. 

The essential quantities in determining the initial state for the hydrodynamic expansion are (i) the 
total multiplicity dN/dy, (ii) the initial time scale tq and (iii) the shape of the transverse distribution. 
Production models differ on all these and our aim in discussing model calculations is to estimate what 
are the present uncertainties in these quantities. We can then vary the initial conditions accordingly 
to explore the range of uncertainty in predicting the thermal photon spectrum in the nucleus-nucleus 
collision at LHC. 

In this subsection, we consider different approaches to primary production: (i) pQCD + final state 
saturation (minijet) |90| and (ii) initial state saturation, the colour glass condensate model |138|, which 
both describe the produced matter as a parton system. We also discuss the string based model Dpmjet, 



of Ranft and collaborators, already presented in Sec. l5.2l (see refs. Il68ll69ll74ll75l ). 



6.51 Perturbative QCD + saturation model 

In the pQCD + saturation model 1901 one assumes that at collider energies the production of final state 
particles is dominated by coUinearly factorized parton (minijet) production above a momentum scale, 
Psat ^ Aqcd> determined from an assumption that saturation reduces parton production below this 
scale and gives only a minor contribution which can be effectively included in the contribution from 
partons above psat- 

The numerical minijet calculation I139II1401 gives the cross sections dai (Ay) / dp^ in the rapidity 
interval Ay for parton species i = g,q and q. For Ay< 1 (around y = 0), the cross section is proportional 
to Ay and the differential cross section can be defined as 

da 1 da{Ay) 1 /" , f . sr^ , , ^2n. . r^2^da'^ 



/ dyi dy2Y,h/A{xi,Q^)fj/Aix2,Q^)-f (40) 

J Ay J ^.f^i at 



dydp\ Ay dp\ Ay 

with nuclear parton distributions /j/^, momentum fractions xi^2> scale Q = Pj, and the partonic cross 
sections da. The basic quantities for the calculation of initial densities in nucleus-nucleus collisions (cf. 
Sec. l6.62l i are the integrated and Pj, -weighted cross sections 

daj 1 p da, (Ay) 1 . . ^ daijpo) 

dPr , , = / dp^— = —ai{po,Ay) = — , (41) 

PO dydp^ Ay Jp^ dp^ Ay dy 

" ri 1 r ^ d^»(Ay) 1 ,^ . . dai{E^){po) 

dPrPr r r = / dp^p^— = —ai{E^){po,Ay) = . (42) 

PO dydp^ Ay Jp^ dp^ Ay dy 

When multiplied with the nuclear overlap function (see Eq. ll02l and appendix 1 for further details), 
ai{po, Ay) gives the average number and ai{Ej,){po, Ay) the average transverse energy of partons i in 
a rapidity interval Ay with transverse momenta above pq. The number of partons is well-defined only 
in the leading order, at higher orders it is not an infrared safe quantity. The transverse energy, however, 
is well-defined, and it has been computed to next-to-leading order in 11411 1142ll . The factor K = 1.6 
from such a computation for the LHC energy^^ is used to estimate the effect of the NLO contributions 
both to the number and the transverse energy of produced partons. In what follows, the number of 
partons is used in estimating the chemical composition of produced partons but the initial energy is 
determined from the produced total transverse energy alone. Nuclear shadowing, as given by the EKS 
parametrization Ii39ll40l l. is included in the calculation of these cross sections. 

To close the calculation, the cut-off momentum pq for central AA collisions is fixed as the satura- 
tion scale, psat, from a geometric condition 



T^.(0)^%^^ = vri?i, (43) 
Ay Psat 



which is solved numerically. We fix Ay = 1. The condition above can be interpreted as giving the 
parton density at which the produced partons begin to overlap leading to fusions which suppress the 
further production of softer partons. The condition can also be interpreted as the beginning of overlap of 
primary interaction volumes that can also be expected to suppress the further yield of partons. 

Once the saturation momentum psat has been obtained from the condition above, the initial produc- 
tion cross sections (Ti(psat, Ay) and ai{Ej,){psa.t, Ay) can be calculated. Although the pQCD-i-saturation 
model is an effective (non-microscopic) approach to describe the bulk of the initial parton production at 
central rapidities of nearly central AA collisions, it provides enough information to fully determine the 
initial state for the hydrodynamic evolution if the initial time is given. With the hydrodynamic evolution 

^^K depends on y/s, scale choice and PDFs. 



11431 it has correctly predicted the multiplicities in central Au+Au collisions at RHIC energies ^/s = 56, 
130 and 200 GeV [ 143 . 11441 . and also given a successful description of the transverse momentum spectra 
of pions, kaons and (anti)protons 11451 . 

Using the notation a{\) = da/dy, the integrated partonic cross sections at ^/s = 5.5 TeV, the 
energy of lead-lead collisions at LHC, are (7g{l) = 135 mb, (7q{l) = 6.46 mb, (7q{l) = 6.14 mb, 
while the Pj,-weighted cross section becomes a{E^){l) = 468 mbGeV. The saturation scale is Psat = 
2.03 GeV. Below in Sec. l6.62l initial energy density and initial values of gluon and (anti)quark fugacities 
are constructed from these cross sections. 

6.52 Initial state parton saturation and color glass condensate 

A much discussed theoretical approach to particle production in heavy-ion collisions at collider energies 
has been based on the assumption that the initial state parton densities saturate and non-linear dynamics 
becomes dominant II138I . This means that the parton (gluon) densities in a nucleus (or in a proton) at 
small X are so high, that gluons interact coherently. From the point of view of the color fields, the high 
density or large occupation numbers of the field quanta can be described as the formation of a "color 
glass condensate" which can be treated in terms of a classical effective field theory 1 146 147J. Quantum 
corrections (to the classical parton distributions of a nucleus) have also been considered II1481 11491 ; see 
also I150II1511 and references therein. 

Applicability of saturation models can be tested, and the saturation scales extracted, in deep inelas- 
tic I A scattering. For the free proton a geometric scaling of the structure function F2 at small values of x 
has been found II152II1531 IT54I. Similar analyses have also been performed for the nuclei II155II156 'I but 
due to the limited amount of nuclear DIS data available, it is not yet fully clear how well the assumption 
of high (saturated) parton density is met in heavy-ion collisions with the nuclei available in nature and at 
the highest collider energy available currently at RHIC or in the near future at the LHC (see also 1 157 1). 

In the effective field theory approach to gluon production in AA collisions, when boost invariance 
is assumed, it becomes possible to choose the gauge in such a way that the problem can be formulated 
as a dimensionally reduced 2-1-1 -dimensional theory. After this the problem is amenable to numerical 
approach II158. I utilizing lattice regularization. For analytical studies of primary gluon production in AA, 
see I1I59I and references therein. 

The earlier numerical calculations with lattice regularization were performed using SU(2) sym- 
metry II 1601 and cylindrical nuclei and the treatment of colour neutrahty was not yet adequate |161 1. 
Recently Krasnitz et al. have formulated the calculation using SU(3) II1621 11631 . more appropriate nu- 
clear geometry for spherical nuclei, and also imposing local colour neutrality in the transverse overlap 
region of collision [ 163 1. Local colour neutrality leads to a rapid decrease of colour field strength outside 
the nucleus and to more a realistic treatment of multiplicity and transverse energy production. However, 
the average transverse momentum, obtained in this model, is too large and has remained a puzzle for 
several years. Quite recently, Lappi has solved this problem, and a factor ~ 2 reduction in the ratio 
E^/N has been discovered [164 J . 

The lattice approach does not, however, give a value for the saturation scale (or the colour source 
density ^u) itself; the overall normahzation must be obtained from elsewhere. For RHIC phenomenology, 
the authors of 111631 suggest two sets of results which now in light of the latest results of II164I lead to 
following, qualitatively different descriptions of the final state: 

• In the case of a smaller scale ^, the total transverse energy ~ g^R^fi^ produced from the 
classical fields approximately equals the measured result, whereas the number of initially produced 
partons (A^ ~ glR\^'^) is only ~half of the multiplicity of hadrons measured in the experiment. 
In this case the only evolution in the final state would be the fragmentation of partons to a couple of 
hadrons on the average. In this picture, which corresponds to the scenario suggested by Kharzeev, 
Levin and Nardi L165J . one would expect the photon and lepton pair emission after the primary 



interactions to be very rare. 
• For a larger saturation scale, the number of partons is close to the measured number of hadrons 
but the initially produced transverse energy is about 2.5 times bigger than the measured one 1 164 1. 
In this case, production must be followed by strong collective expansion, initially dominantly in 
the longitudinal direction, which allows for a transfer of energy into the longitudinal motion. This 
case corresponds to the evolution suggested by pQCD+saturation+hydrodynamics [ 143 1. 

The two cases described above represent the extreme ends of the possibilities for the final state 
and real collisions can be somewhere between these two. In general, the production of partons would be 
followed by some degree of fragmentation, then by thermalization stage and finally by the hydrodynamic 
expansion. How thermalization may occur is studied e.g. in ri66l 11671 fT68i il69i iTTOl 11711 IT72 1. In 
Sec. 16.621 when discussing the range of initial conditions for the hydrodynamic expansion, we take, 
at fixed multiplicity, these alternative scenarios to indicate different choices of the initial time for the 
hydrodynamic expansion. 

In the model of Kharzeev, Levin and Nardi [165] the partonic final state at RHIC energies corre- 
sponds to the first alternative of Krasnitz, Nara and Venugopalan described above. Like in the original 
form 1 138 1, parton saturation is formulated in the initial state in terms of unintegrated gluon distribu- 
tions. The rapidity distributions of produced gluons are computed as 2 ^ 1 gluon collisions and they 
are found to be proportional to the integrated gluon distribution and the number of participant nucleons. 
In this scenario, the produced gluons just fragment into a couple of hadrons and there is no strongly 
interacting medium, no collectivity nor any energy transfer from transverse to longitudinal motion: the 
pseudorapidity distribution of final hadrons is directly that of produced gluons. With negligible amount 
of secondary collisions, photon and lepton pair emission after primary interactions should be negligible. 
(This, of course does not concern the decays of final hadrons which produce plenty of photons and lepton 
pairs.) 

The energy dependence in the model 11651 enters through the dependence of the saturation scale 
on the center-of-mass energy, obtained in I1521ll53lll54l from DIS of the free proton. Also the number 
of participants is y/s dependent. The overall normalization in |165| is, however, left free. If one as- 
sumes that the ratio of the multiplicity of final hadrons to that of the initially produced gluons is energy 
independent, the normalization can be fixed at, say, the RHIC energy ^/s = 130 GeV, and predictions 
can be made for dNAA/drj for final hadrons at other energies. In this way, the model 11651 predicts the 
charged particle multiplicity for central Pb-i-Pb at ^/s = 5.5 TeV to be dNch/dr] ~ 2200. This number 
is somewhat smaller than in the pQCD-i-saturation model 1 90. .143 1 and it is included in the multiplicity 
range considered in Sec. 16.621 

6.53 Particle production in DPMJET dual parton model. 

In the Dpmjet Monte Carlo model all hadrons are produced from soft and hard chains, soft chains 
resulting from soft pomeron exchange and hard chains (jets and minijets) from (LO) perturbative QCD. 
The Dual Parton Model (Dpmjet) is treated in detail in Sec. 15. 21 The assumption is that the initial energy 
density is to be defined by the soft and hard chain end partons before they decay into hadrons and not by 
the final state hadrons as in the original Bjorken formula lll73l . The information on the transverse location 
and the energy-momentum 4-vectors of the chain-end partons are stored for each event. Since Dpmjet 
describes the evolution of the final state in the momentum space, a connection between the rapidity y 
and the longitudinal space-time rapidity i], as well as the production timescale, must be specified. Below, 
when discussing the initial conditions for the hydrodynamical evolution in detail, we compare the energy 
densities from Dpmjet and the pQCD -i- saturation calculation. 

6.6 Hydrodynamic expansion of thermal matter 

H. Niemi, P.V. Ruuskanen, S.S. Rasanen, K. Redlich, D.K. Srivastava 



The dynamic evolution of the expanding matter in local thermal equilibrium is governed by the 
hydrodynamic equations 

which express the energy and momentum conservation as they are transferred by pressure gradients and 
flow from one region to another during the expansion. The energy-momentum tensor 

T^"" = (e + p)u''u'' - pg^"" (45) 

is expressed in terms of e, the energy density, p, the pressure, and u'^, the 4-velocity of flow. At LHC 
energy the net baryon number in the central rapidity region is very small and will be neglected in the 
following. 

At LHC the total rapidity interval is large and, over a few units, the central rapidity region is antici- 
pated to be relatively flat and it should be reasonable to impose boost invariance in the central rapidity re- 
gion. The rapidity of flow is then equal to the space-time rapidity, ryfiow = rj = (1/2) log [(t + — z)], 
not only for the initial conditions but throughout the expansion. We will also restrict the calculations to 
zero-impact-parameter collisions with azimuthal symmetry. The densities and the transverse flow ve- 
locity depend then only on r = ^/t"^ — and r and can be expressed as e = e(r, r) and Vr{T, r). This 
simplifies the equations greatly and numerical solutions are easily obtained once the distributions e(To, r) 
and Vr{TQ,r) are fixed at initial time tq and the Equation of State (EoS) is specified. In all calculations, 
the initial transverse velocity is taken to be zero, Vr{To,r) = 0. 



6.61 The Equation of state 

To solve the hydrodynamic equations an Equation of State (EoS) has to be specified. We study the 
expansion using different EoSs. We investigate the expansion and the photon emission also when quarks 
and gluons are not in chemical equilibrium. The deviation from chemical equilibrium is described using 
a multiplicative fugacity for gluons and (anti)quarks. For the collisions at LHC we assume that the net 
baryon number is zero. 

We mainly use an EoS which describes the high temperature phase as an ideal gas of massless 
quarks and gluons both in kinetic and chemical equilibrium. The low temperature phase is a gas of all 
hadrons and hadron resonances with masses below 1.4 GeV interacting with a repulsive mean field. This, 
as well as a bag constant, is needed for consistent treatment of transition between the two phases as a 
first order phase transition. With A'^^ = 3 the bag constant B and the mean field constant K are chosen 
to be B^/'^ = 243 MeV and K = 450 MeV fm^ giving Tc = 167 MeV for the transition temperature. In 
Ref. G74| this is called the EoS A. 

The expected gluon dominance in parton-based approaches of primary production means that even 
if the produced particles are close to kinetic equilibrium, they will be far from chemical equilibrium. We 
describe this by using multiplicative fugacity factors in the EoS. The equihbrium energy density 

e = eq + tq + eg + B (46) 

with ei{T) = ■iaiT'^, Ug = Wir'^/dO and Og = Og = (21/4)iV/7r790 is then replaced by 

e(r, Xi) = Xqeg + Ag-eq- + A^e^ + B{Tc, X^) . (47) 

Since the net baryon number will be small at central rapidities as indicated also by the result ag ~ ag 
from the minijet calculation at the LHC energy, we will take Ag = Ag. 

A similar expression can be written for the pressure showing explicitly that the relation e*^*^ = 3^*^^ 
remains unchanged for the thermal parts of pressure and energy density. At large temperatures, say 
T > 2Tc, the thermal parts dominate, indicating that the evolution of the flow at high temperatures is 




Fig. 52: Energy density and pressure as function of temperature for 
EoS A and a parametrization of lattice data. 
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unaffected by fugacities becoming different from 1 . However, the relation between the temperature and 
the energy density depends strongly on fugacities: small fugacities indicate a decrease in the number of 
effective degrees of freedom and an increase in the temperature at fixed energy density. This has a large 
effect on the photon emission rates. 

As discussed above, the value of the bag constant B is chosen so that Tc = 167 MeV in chemical 
equilibrium. Obviously this can also be done with fugacities different from 1. In this case B depends 
not only on Tc but also on the fugacities. As there are no counterparts of fugacities in the hadron gas 
we take the phase transition temperature Tc to be a constant, independent of fugacities. This avoids the 
ambiguity of having a phase transition to hadron gas at different temperatures. Once Tc is fixed, the bag 
constant will depend through the Gibbs criterion on fugacities only: 

Pug{Tc) = PQGpiTc, Ai) - B{Xg, Xg) . (48) 

It turns out that the evolution of fugacities is fast enough so that, in the central region of the transverse 
plane, they approach values close to 1 before the temperature reaches Tc. This supports the choice of 
constant Tc. 

This kind of description of phase transition cannot be made with arbitrarily small fugacities. As 
will be discussed in the next subsection, a straightforward application of the minijet results for the parton 
cross sections would lead to very small values of fugacities at the edge of the nuclei in the transverse 
plane. To avoid a situation where the effective number of degrees of freedom in the parton matter 
becomes so small that the phase transition from the hadron gas to parton matter would not occur at all, 
we have to limit the values of fugacities from below. The numerical values of these limits depend on the 
details of the hadron gas phase. If the hadron gas consists of all hadrons and hadron resonances with 
masses below 1.4 GeV, we take limits to be Xg > 0.60 and Xg > 0.17. 

We also consider the expansion using results for the EoS from lattice calculations. As a starting 
point we use a parametrization of lattice results on the pressure and energy density as a function of 
temperature at T > Tc 1.175.1 . We join this parametrization with the EoS of Hadron Resonance Gas 
(HRG) at Tc. We constrain the parametrizations to be consistent with the thermodynamic relations 
s = (e + p)/T = dp/dT. The resulting parametrization, EoSL is shown in Fig. |52lfor the energy 
density and the pressure together with e and p in the EoS A. Since for the energy density the ideal-gas 
Stefan-Boltzmann limit is reached from below on the lattice, the temperature at given energy density is 
higher than in the case of EoS A when the limit is reached from above. The ratio of energy density to 
pressure in plasma is quite similar in both cases when T is not in the vicinity of Tc. 

It turns out that the high temperature region completely dominates the total photon emission from 
the plasma. For this reason the difference in the emission for EoS A and EoS L comes from the difference 



in the temperature when the calculation is done starting from initial conditions with the same energy 
density. As will be seen below, the difference in the spectrum is visible but not very significant. 



6. 62 Initial conditions for hydrodynamics 

Models of particle production are formulated in momentum space. To obtain the initial conditions needed 
to initialize the hydrodynamical calculation, we must correlate the formation time tq and the position 
with the momenta of produced quanta. 

At high collision energy the incoming nuclei are strongly Lorentz contracted and the time it takes 
for them to pass through one another is short. The collision region can then be approximated as a 
plane perpendicular to the beam direction, z = 0, at time t = 0. Particles produced in the collision 
will move with a longitudinal velocity Vz = z/t or rapidity equal to the space-time rapidity, y = rj = 
(1/2) log[(t+z)/(t— z)]. Combining this with the formation time and incorporating the nuclear geometry 
allows us to calculate the produced energy in a volume element dzd^r = Tod-qd^r from the momentum 
spectrum of produced particles. 

Depending on the approach, nuclear geometry enters the calculation somewhat differently. In the 
pQCD + saturation calculation, based on factorization of hard processes, the momentum distribution of 
produced partons per unit transverse area is 

1 dNir) da , 

^ ' - [TA{r)? , (49) 



d^r dydpj. dydp^ 

where the nuclear thickness function r^(r) is given in terms of pA{z,r), the nucleon density in the 
nucleus A, as the integral ^^(r) = f dzpA{z,r). Nuclear geometry enters here explicitly as the factor 
[T^(r)]^ expressing the nucleon-nucleon luminosity per unit transverse area. Nuclear size appears also 
in the calculation of the cross section both through the saturation condition (l43l and nuclear shadowing 
in the structure functions. Since E^, = p^, for massless partons, the energy density of produced partons 
in a volume element with the longitudinal size Todrj = Tody is obtained as the -weighted integral of 
the parton distribution: 

1 f , dN{r) 



'^^'^'^ = ^rj "^P-P-d^- ^^^^ 



Using Eq. (l42b . we obtain from the pQCD -i- saturation model 

, , fr, da \ [TA{r)f 1 ,^ , , , jT^(r)]^ 
^['TQ.r) = ( / dpj,pj,-j--j— ] = — (T(£;j,)(psat, Ay) . (51) 

V-^Psat 



dydp^J To Ay tq 



This formula shows that in the minijet calculation the production is proportional to the number of binary 
colUsions per unit transverse area. 

As mentioned above, in the DPMJET Monte Carlo model the initial energy density is defined in 
terms of the soft and hard chain-end partons before they decay into hadrons. A grid is introduced in 
the transverse plane of the collision and the Bjorken formula is used to calculate the transverse energy 
of chain-end partons separately in each of these transverse space bins. This is done separately in bins 
of different transverse momentum using formation time in the Bjorken formula which depends on the 
transverse momentum, r 1/pt- The formation time is normalized to r = 1 fm/c in the lowest p^, bin 
(0 < pt < 0.8 GeV/c), which contains essentially all the soft chain ends. 

The total energy is obtained by summing the contributions from different transverse energy bins. 
In Fig.|53lwe compare the energy densities from DPMJET and the pQCD -i- saturation calculation. At the 
maximum the DPMJET result is 15 % below that from the minijet calculation and within the uncertainties 
of the models they agree quite well. It should be mentioned that the DPMJET calculation gives also the 
rapidity distribution over the whole rapidity range, whereas the minijet calculation cannot be extended 
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Fig. 53: Transverse dependence of the energy density 
e(r, r) at time r = 1.0 fm/c from the pQCD+saturation 
model compared with results from the PHOJET calculation 
(circles) based on dual parton string model. 
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over more than a few units at central rapidities. Since our hydrodynamical calculation is based on boost 
invariant flow in the longitudinal direction, we cannot utilize the rapidity dependence of the DPMJET 
result. 

In the models based on initial state saturation, production is related to the number of participants 
which for central, zero-impact-parameter collisions is proportional to r^(r). Details in implementing 
the production dynamics may modify the transverse dependence but we will take the binary collision 
density [T^(r)]^ and the participant density T^(r) to represent the range of variation in the shape of 
transverse distribution of initial energy density. 

We next consider the choice of initial time tq for the hydrodynamic evolution. In models based on 
parton production, the formation time scale of primary partons is the inverse of the saturation momen- 
tum. To ^ 1/psat- However, as was discussed earlier, the time needed for the system to thermalize, rth, 
can be longer than the production time of partons 1 169 170 171 1, especially if the primary production is 
followed by fragmentation. Usually it is argued to indicate that thermal and hydrodynamical description 
of expansion cannot be started earlier than at rth- On the other hand, collisions drive the system towards 
equilibrium and they cannot be ignored since they contribute both to the buildup of collective motion and 
the emission of photons and lepton pairs even before full thermalization has been achieved. Applying 
hydrodynamics with initial time tq ~ Tprod < "^th means assuming that it provides a reasonable approxi- 
mation for treating the effects due to collisions which drive the matter to thermal equilibrium. Especially, 
the high transverse momentum photons might otherwise be underestimated since the collision energies 
of partons are largest at the earliest stage. 

In the pQCD -i- saturation model the assumption of early thermalization can be supported with two 
arguments: First, for the kinetic pressure and energy density of massless particles the relation p = e/3 
holds for any isotropic momentum distribution. Second, in the pQCD -i- saturation model the number 
density and the energy density of produced partons is approximately consistent with the thermalization 
assumption. This would indicate that the collision rate can be similar to that in thermal system with the 
same energy density. 

In the pQCD minijet calculation the saturation momentum at ^/s = 5.5 TeV is 2.03 GeV giving 
To = 0.1 fm/c. This is the earliest initial time to be used in exploring the photon emission. 

Results for the multiplicity are very similar from the minijet model and the DPMJET calculations 
as implied by the similarity of initial conditions, Fig.|53] They both give values dNch/dy ~ 3000 for the 
total multiplicity and agree within the uncertainties of the models, like the choice of the saturation scale 
in the minijet calculation or the possible entropy generation after the primary production. 

In the numerical evaluation of classical field equations the saturation scale is a parameter and the 
multiplicity changes with different choices of the scale. 

Calculations based on the initial state saturation formulated in terms of the initial gluon densities 



11651 predict the energy dependence of the multiplicity and once the normalization is fixed from the 
RHIC data, the prediction for the charged multiplicity is around dNch/dy 2200 at LHC. Even though 
the authors indicate that, at least at RHIC, the evolution in final state, other than fragmentation, is not 
important, we exploit their result by considering an initial distribution oc 7a (r) normalized to fit the 
predicted multiplicity from their model. 

Results from classical effective field theory on lattice in the transverse plane can change, depend- 
ing on the choice of the saturation scale, by a large factor (order of 2 at RHIC energy) 11631 for the 
parton multiplicity. The larger multiplicities are comparable with those from the minijet calculation and 
the lower ones lead to a final state similar to the one from the saturation calculation in terms of initial 
gluon distribution of incoming nuclei [165 1. To get an idea how much the multiplicity may affect the pho- 
ton emission we scale the initial conditions also upwards to cover the interval 2000 < clNch/dy < 4000 
or from ~ 3000 to ~ 6000 in the total multiplicity. 

Even though the final state energy density at r = 1 fm/c from the DPMJET model is quite similar 
to other calculations as shown in Fig.|53j it is not clear that the photon emission should be the same. In 
the Dpmjet calculation also strings are present in the final state and the space-time evolution might be 
very different from that of parton gas, predicted from the other calculations. Neither is it clear how well 
the use of the photon emission rates based on pQCD calculation with parton degrees of freedom apply if 
strings play an important role among the degrees of freedom in the final state. 

To summarize: The two most important parameters for the determination of the initial state, the 
multiplicity dNch/dy and the initial timescale tq, have a considerable range of uncertainty. Also the 
maximum values of the temperature in the final state are affected by the transverse shape of initial dis- 
tributions and since the photon emission has a very strong temperature dependence, these effects are 
seen in the calculated photon spectrum. To estimate the uncertainty in the calculation of thermal photon 
emission we explore the following values of the input parameters for the initial conditions: 

• 2000 < dNch/dy < 4000, 

• O.lfm/c < Tq < 1.0 fm/c and 

• e(To, r) oc [TA{r)]'^ for k = 1 and 2. 

6.63 Initial conditions and the evolution of fugacities 

Gluon dominance in the initial production will lead to small (anti)quark fugacities. This turns out to have 
an important effect on the emission rates. At fixed temperature the decrease of fugacities leads to smaller 
densities which decreases the emission rate. However, at fixed energy density the decrease of fugacities 
increases the temperature which in turn increases the emission rate. Thus the effects from the lack of 
chemical equilibrium are complicated and the question if the emission rate increases or decreases can 
only be answered with a detailed calculation 1 176 1. This means that we have to determine not only the 
initial values and the evolution of the temperature T but also of the fugacities Ag, and Xq. 

We determine the initial values of fugacities in the transverse plane at the initial time tq from the 
information on jet cross sections ag and cTq aq and the cross section for transverse energy production, 
a{Ej.). In doing so we also have to impose limits for the values of fugacities to assure thermodynamic 
consistency. 

In Sec. 16.621 we constructed the initial energy density from a{E^), Eq. (15 It . We can define the 
number densities of gluons and (anti)quarks, Ug and Ug respectively, in the same way by using ag and 
Uq instead of a{E^). If the temperature is known, the fugacities can be determined from the relation 
between the equilibrium thermal density nj(T) and the real density n^: 

Ui = Xihi{T), i = g,q. (52) 

The equilibrium density nj(T) is given as hi = biT^ with bg = 16C(3)/7r^, (^(3) = 1.20206 for gluons 
and bq = bq = 9Nf({3)/27r'^ for (anti)quarks. In fact, to determine the fugacities and the temperature the 




Fig. 54: Evolution of fugacities starting from tire pQCD + 
saturation initial state. Ag(r, r) and Aq(r, r) are shown as 
function of time for r = and r = 6 fm. 



Fig. 55: The evolution of temperature for the same situations 
as for the fugacities in Fie. 1541 



equations both for the number densities and the energy density, Eq. (147 1 must be solved simultaneously. 
An extra complication arises from the fugacity dependence of the bag constant B. 

We will not go to all the details of solving the equations for the temperature and the fugacities 
here; it obviously can be done. However we must discuss the range of values for fugacities for which 
a consistent equation of state can be constructed. In Sec. l6.6ll we stated a lower limit for the fugacities 
to ensure that the hadron gas will be the stable low temperature phase. It turns out that if we solve the 
fugacities and the temperature from the initial conditions using Eqs. d52b and d47l we will obtain, at 
the edges of the transverse nuclear overlap, fugacities which violate this condition. Also, in the central 
region of the transverse plane, gluon fugacities with values > 1 are obtained. If the fugacity is 
written in terms of chemical potential, values exceeding one will lead to a singular unphysical phase 
space distribution function. To avoid this, we limit the gluon fugacity to be < 1. If solving the 
Eqs. d52t and Wfl gives Ag > 1 we take A^ = 1 or equivalently Ug = ng{T) and keep the quark fugacity 
unchanged. New temperature and number densities can then be recalculated using Eqs. d47l and (I52t . 
The lower limits at the edges of the nuclear overlap are imposed in the same way. It turns out that when 
these physically motivated limits are applied at the initial time, they are not violated by evolution of 
fugacities during the expansion on the parton matter. 

The evolution of fugacities is determined by the rates of number changing reactions for gluons 
and quarks. There is a considerable uncertainty in the parton cross sections which are needed in the 
rate calculations. We follow the studies of Biro et al. lOTl and Elliott and Rischke l92l ') and include the 
processes gg —>■ qq and gg —>■ ggg as well as their inverse reactions into the rate equations. The rate 
equations for the gluon and quark fugacities can be written as (see Eqs. (23) and (24) in Ref. 1 9110 



d^{Xini{T)un = Ci{T,Xg,Xg), i = g,q. (53) 



with the collision terms 



Cg = hgRsXgil - X g) " 2hgR2Xg (^1 " , (54) 

Cg = hgR^Xg 1 - ^ ^ , (55) 



where, even though we have taken Xg = Xq, we have written the quark and antiquark fugacities expUcitly 
to show the origin of different terms. For the rate parameters Ri we use the results derived in L91J : 

i?2 w 0.24:NfalXgT\n{l.65/asXg) , (56) 
i?3 = 1.2alT{2Xg - Xiy/\ (57) 

Equations (I53l57t are solved numerically simultaneously with the hydrodynamic equations and the 
change in the EoS due to the change in fugacities is included in the calculation. 

In the beginning of the expansion with large parton densities the evolution towards thermal equilib- 
rium is fast but it slows down as the matter expands and rarefies. In Fig.|54]we show the time dependence 
of the fugacities in the parton matter at r = and r = 6 fm starting from the pQCD + saturation initial 
state at tq = 0.1 fm/c. At r = gluon fugacity Xg = 1, but since the number of quarks and antiquarks is 
small, the rate for gg —>■ qq reduces the number of gluons faster than the the inverse reaction and reaction 
99 ~^ 999 can balance it. 

The evolution of fugacities affects also the evolution of the temperature. In Fig.|55lthe temperature 
is shown as function of time for the same conditions as the fugacities in Fig.|54] To see the significance of 
the lack of chemical equilibrium, we show the evolution of the temperature for the same initial conditions 
also when full equilibrium is assumed. Initially the temperature is higher when the matter is out of 
chemical equilibrium. This happens because the effective number of degrees of freedom is l6Xg + 
2\NfXg/2 instead of 16 + 2\Nf/2 of the equilibrium case. As the number of effective degrees of 
freedom increases with time the temperature drops faster in non-equilibrium matter and approaches that 
in equilibrium matter. However, as will be seen below, from the point of view of photon emission, the 
higher initial temperature turns out to approximately compensate for the lower parton densities. 



6.7 Hadron spectra from a hydrodynamical calculation 

H.Niemi, P.V. Ruuskanen, S.S. Rdsdnen 

The calculations of hadron and photon spectra proceed somewhat differently. Since the mean free 
paths of photons are larger than the size of the fireball, they are emitted throughout the expansion of 
matter from the whole thermal volume. Hadrons escape only from the surface of the fireball or after 
the density in the whole volume is so small that no further interactions take place. We first consider the 
hadrons and then the photons. 

As the matter expands, distances between particles become large, collisions cease and momentum 
distributions freeze out. The condition for the freeze-out to happen is usually expressed locally in terms 
of the energy density or temperature reaching a given value. This determines a three-dimensional freeze- 
out surface a^{x) in space-time. The prescription of Cooper and Frye 11771 to convolute the flow and 
the thermal motion along the freeze-out surface 

oo „ 

= rT[-p,h{n^rv/f)K,{n^r'^)dT 

+mTlo{njrVr^)Ki{njr^) dr] , 

is used to calculate first the spectra of all hadrons and hadron resonances. The second expression above is 
valid for cylindrically symmetric, boost invariant flow with Vr the radial flow velocity, 7^ = 1/y^l — t;^, 
and Kq and /q Bessel functions. We then follow the chains of all possible two and three body decays 
and collect the spectra of final stable hadrons. The stable hadrons can be interpreted as the absolutely 
stable hadrons or also those with weak decays or both weak and electromagnetic decays. The treatment 



(58) 



(59) 



of hadrons with weak decays can be important, e.g. in the case of A's, depending if the measurement can 
separate the feed-down from heavy hyperons or not. Electromagnetic decays, as a source of photons, are 
of special interest but we display later also the spectrum of vr'^'s before their decay. 
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Fig. 56: Transverse momentum spectra of positive 
pions, kaons and protons measured by PHENIX Col- 
laboration at -y/s = 130 GeV shown with results 
from a hydrodynamical calculation with initial con- 
ditions from the pQCD + saturation model. 
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Since the transverse momentum distribution of the decay photons depends strongly on the trans- 
verse spectra of hadrons, we show in Fig.|5^the measured spectra of positive pions, kaons and protons at 
RHIC I178II with results from a hydrodynamical calculation using the initial conditions from the pQCD 
-I- saturation model I145II . The results describe the data quite well at small momenta, but it is clear that a 
hydrodynamic calculation cannot describe the hadron spectra at transverse momenta larger than ~ 4 — 5 
GeV. The fraction of hadrons with p^, > 5 GeV from all hadrons is small and in thermal models the domi- 
nant fraction of low transverse momentum hadrons is assumed to come from the thermalized matter with 
collective motion. From this point of view the calculation of decay photons using the hadron spectra 
from a hydrodynamical calculation gives a reasonable estimate in the range of low momenta, say up to 
~ 3 _ 4 GeV. 

Similarly, the spectra of hadrons from the thermal fireball at LHC energy can be calculated for any 
given choice of initial conditions. As an example of hadron spectra at LHC we show here the spectra 
of positive pions and kaons and protons in Fig. |^ obtained with initial conditions from the pQCD 
-I- saturation model L179J which, as already pointed out above, reproduce well the data at RHIC. To 
show the uncertainty in the shape due to the decoupling procedure, results with decoupling temperatures 
Tdec = 120 and 150 MeV are displayed. All spectra become flatter at later decoupling and the role of the 
flow is seen from the change getting larger with increasing mass. We will show later how this uncertainty 
is transferred to the spectrum of photons from hadron decays. 

Note that even though the range of parametrization of initial conditions have been motivated by 
different model calculations it does not mean that such spectra are predictions by the authors of these 
models in some of which hydrodynamic expansion may not be needed at all. When using a thermal model 
for describing the hadron spectra and the electromagnetic emission we have used these calculations to 
argue what could be the reasonable range of input parameters and we always make the assumption that 
there is a hydrodynamical expansion stage between the primary production and the formation of final 
hadrons. As reviewed in Sec. 16.521 e.g. in 11651 no thermal stage is assumed to occur 
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Fig. 57: Transverse momentum spectra of posi- 
tive pions, kaons and protons calculated at ^/s — 
5.5 TeV from the pQCD + saturation initial con- 
ditions 11791 . The two sets of curves correspond 
to two different choices of decoupling temperature, 
Tdcc = 120 and 150 Me V. 
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6.8 Spectra of thermal photons 

H. Niemi, P.V. Ruuskanen, S.S. Rdsdnen, D.K. Srivastava 

The photon spectra are obtained by integrating the emission rate over the whole space-time vol- 
ume of the fireball. We use the emission rates which are discussed in detail in Secs. l6.2l for the equilibrium 
QGP, in Sec. l6.3l for the QGP in kinetic equilibrium but out-of-chemical equilibrium and in 16.41 for the 
hadron resonance gas. 

The local emission rate is a function of the local temperature and the energy of the photon in the 
rest system of the matter, dN{uj* ,T)/d'^x. In the case of flowing matter this is the co-moving frame 
of the fluid element. If p^^ = {E, p) is the four-momentum of the photon in the fixed frame where the 
photon is observed and the four-velocity of the emitting fluid element in that frame, the energy of the 
photon in the co-moving frame is ui* = p-u = p^u^. Since d?\)/ E is Lorentz invariant, we obtain 



At the transition temperature Tc the matter may enter into a mixed phase which is characterized by the 
volume fractions of each phase. In the mixed phase the emission rate is taken to be the sum of rates in 
QGP and hadron gas at T = Tc, weighted with the volume fraction of the phase. 

We next consider the photon spectra which are obtained when n^(r, r) and T{r, r) are calculated 
using different initial conditions and equations of state. We start with the spectrum, shown in Fig. |58j 
from the minijet calculation with the pQCD -i- saturation model, used as a baseline and plotted in most 
of the figures to facilitate the comparison of different results. In this baseline case we assume full 
equilibrium in the plasma phase, both kinetic and chemical, the charged multiplicity is dN ch/dy = 2900 
and we have taken tq = 1/psat ~ 0.1 fm/c. A similar calculation at RHIC energy, y/s = 200 GeV has 
been presented in Ref. 1 180|. 

The large-py region of the spectrum is completely dominated by the emission from plasma during 
a short period at earliest times. The transverse flow plays no role at that time and the slope is determined 
by temperatures close to the maximum temperature. With these initial conditions plasma gives the main 
contribution also at smaller values of transverse momentum. Photons from the hadron gas phase have a 
much steeper slope than those from QGP. This slope is not given by the temperature of the hadron phase 
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Fig. 58: Our baseline spectrum of thermal phiotons: full 
equilibrium with EoSA, dNch/dy = 2900 and to = 0.1 
fm/c. Contributions from plasma and hadron gas are shown 
separately. 
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alone but by the combination of temperature and the flow. Actually the situation is quite similar with that 
of the hadron spectra. This is shown by the fact that the slope of the (almost massless) pions in Fig. |^ 
is very close to that of photons from hadron phase in Fig. 158) 

In Fig.|59lspectra of thermal photons are shown for charged multiplicities dN ch/ dy = 1900, 2900 
and 3800. To see the effect from the change of the multiplicity alone, same value tq = 0.1 fm/c for 
the initial time has been used in each case. As will be seen later, the transverse momentum window 
of interest for observing thermal photons may extend from ~ 2 GeV up to ~ 10 GeV in the 
favourable case. A change in charged particle multiplicity from 1900 to 3800 increases the photon rate 
by more than a factor of five at the lower end and almost a factor of 10 at the upper end of the interval. 
Since the dependence of decay photons is close to linear in hadron multiplicity, a strong dependence 
on the multiplicity in the relative strength of thermal photons to those from hadronic sources can be 
expected. 

A large uncertainty in the photon emission comes from the uncertainty on the production and 
thermalization times of the initial partons. In all production calculations based on QCD the initial state 
is dominated by gluons due to the SU(3) properties of the couplings among gluons and (anti)quarks. 
Thus the system is far from chemical equilibrium after the production stage. It may also be far from 
kinetic equilibrium and there are arguments that equilibrium is reached only on a timescale of 2. . . 3 fm/c 
I169lll70lll71l . We will consider the chemical equilibration separately below. Here we just consider the 
effect of the initial time r o in the hydrodynamical calculation. 

We already argued in favour of a shorter timescale for the initial time of hydrodynamic expan- 
sion. To, than obtained from the thermalization arguments. In a parton-based approach the produc- 
tion timescale is provided by the momentum scale at which the production is expected to saturate, 
"^prod ~ 1/Psat- The momentum distribution will not be thermal initially but the partons can scatter 
from the time Tprod on and some of these interactions produce photons. A kinetic calculation of photon 
emission starting from different initial distributions would be better justified but we would like to argue 
that by taking tq = Tprod we effectively simulate the photon emission also from those interactions which 
bring the matter towards kinetic equilibrium. In Fig.|60lwe show the variation when the initial time is 
changed. It is seen to be very strong at large p^, as is expected since the magnitude of local relative 
momenta or, in the thermodynamic language, the temperature decreases rapidly due to the longitudinal 
expansion. 

When we compare these spectra with photons from other sources, we will see that the observa- 
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Fig. 59: Thermal photon spectra in Pb+Pb collisions at LHC 
energy for different charged hadron multiplicity assuming 
full thermal equilibrium throughout the expansion. 
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Fig. 60: As in Fie. l59l but as a function of initial time tq for 
fixed charged multiplicity dNch/dy — 2900. 



tion of thermal photons, or more generally photons from secondary collisions among initially produced 
partons, is possible only if the produced partons can begin to scatter at early times immediately after 
production when they still have large relative momenta. 

In calculating the spectra in Figs.|59land|60|we have assumed that the radial dependence of initial 
distributions is proportional to [TA(r)]'^, a = 2. Fig. |^ shows the difference between the cases a = 
1 and 2. With a = 2 the distribution is steeper and more peaked in the center than in the case a = I. 
Consequently, for fixed multiplicity the maximum temperature is higher in the former case and the strong 
temperature dependence of emission rate leads to stronger overall emission and more shallow dependence 
on Pj, of the spectrum. 

We next discuss how the deviation from chemical equilibrium affects the photon emission. The 
initial values and the evolution of fugacities are determined as described in Sec. 16.631 When the values 
of fugacities are less than one, Aj < 1, the numbers of degrees of freedom are effectively less than 
in plasma in full equilibrium. Comparing with the equilibrium case at a given energy density, parton 
densities are smaller but temperatures higher. For the photon emission spectrum, it turns out that the 
increase in temperature approximately compensates the decrease of parton densities, over the range 
considered. 

In Fig.|62lwe compare our baseline case in full equilibrium with the non-equilibrium case specified 
in Sec. 16.631 In both cases the charged particle multiplicity is 2900 and the initial time tq = 0.1 fm/c. 
The spectra approximately coincide, the off equilibrium case giving slightly smaller results at small 
and about 20% higher at Pj,>10 GeV. 

We also show in Fig. |62l the result of a calculation using an equation of state with the high tem- 
perature phase fitted to the lattice data. The calculation starts from the same initial conditions as in the 
baseline case. As is seen in Fig. |52j the lattice parametrization corresponds to higher temperature at 
given energy density than our EoS A. In the photon spectrum this shows up as a slight increase in the 
photon emission. There are indications that if lattice results are interpreted in terms of quasiparticle exci- 
tations, the quasiparticle masses do not exactly correspond to the thermal masses of partons. The photon 
emission is calculated with the standard perturbative rates but it might be more consistent to replace the 
thermal masses in the rates with the (temperature dependent) quasiparticle masses obtained from the 
lattice information. 
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Fig. 61 : Photon spectra when the shape of initial condition is 
is taken to be [TA(r)]'^ with a = 1 or 2. dN^^/dy ~ 2900 
and To — 0.1 fm/c in both cases. 



Fig. 62: Our baseline result (solid line) is shown with the 
the photon spectrum for an out of chemical equilibrium cal- 
culation (dashed-dotted line) and for a calculation using the 
lattice EoS (dotted line). In all cases dNch/dy ~ 2900 and 
To = 0.1 fm/c. 



Recently an interesting new plioton source was proposed I181II . It is an outcome of the interactions 
of energetic jets with the thermalized plasma. The idea is similar to that in jet quenching: the bulk of 
the plasma comes from partons close to the saturation momentum but the high energy partons do not 
thermalize. They will traverse the matter interacting with its quanta. In some of these interactions 
photons will be emitted and if the properties of matter and the cross section for photon production are 
known, the emission rate can be calculated. 

In Fig. the spectrum of photons emitted from jets traversing the plasma is compared with 
our baseline thermal spectrum. The calculations of jets interacting with plasma has been performed in 
the case of no transverse expansion and this complicates the comparison with the results from thermal 
emission. In any case, the emission from jet + plasma interactions is related to the existence of plasma; 
it is not a competing signal with thermal emission but they add together, enhancing the possibility to 
observe a photon excess due to the formation of plasma. 

6.9 Comparing thermal and decay photons in a hydrodynamical model 

H. Niemi, P.V. Ruuskanen, S.S. Rdsdnen 

Hadron decays, in particular those of vr'^'s, form the most copious source of photons. In this section 
we show the photon spectra from the hadron decays when the hadron spectra are obtained from the 
hydrodynamical calculation and compare them with the spectra of thermal photons. The decay photons 
from all hadrons in the EoS are included. In a hydrodynamical model the thermalized part of produced 
matter is assumed to dominate the multiplicity and the transverse momentum distributions up to few 
GeV. A hard parton, say Pt>10 GeV will not thermalize but in traversing the thermal part of matter 
it may loose energy before fragmenting into the jet of final hadrons. The loss of energy of these hard 
partons leads to a shift to smaller transverse momenta of the large-p^, hadrons. The jets still dominate 
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Fig. 63: Photons from jets interacting with plasma. Ther- 
mal photons dominate at small but have a steeper slope 
than those from jets. Since the jet + plasma calculation is 
performed with — 0, we show also the thermal spectrum 
in that case. 



at large transverse momentum part of the hadron spectra. In this picture decay photons with Pj,<2 — 4 
GeV come from the hydrodynamic part of the hadron spectra and the large-p^, photons from hadrons 
produced in jets. There is no obvious way to join the two regions smoothly since there is no well-defined 
way to cut off the jet calculation at low transverse momenta. 

We would like to remind that in the hydrodynamical model thermal photons are emitted through- 
out the expansion of the produced fireball whereas the final hadrons emerge only from rarefied low- 
temperature matter on the decoupling hyper surface. Consequently the large transverse momentum pho- 
tons and hadrons come from different regions of the fireball: high-pj, photons are emitted at the earliest 
times from the initial quark-gluon matter at highest temperature reached in the interior of the fireball. 
On the other hand, the strong flow develops at the surface of expanding matter and the large-p^, hadrons 
originate from the fastest matter at the edge of the fireball. Expressed differently, the slope of the tail of 
thermal photon spectrum is determined by the highest values of the initial temperature in the interior of 
the matter but those of hadron spectra by the fastest flow on the surface of the matter. 

We show the calculation from the minijet initial conditions only. The initial time is tq = l/psat> 
equal to 0.1 fm/c at the LHC energy, ^/s = 5.5 TeV, and 0.18 fm/c at the RHIC energy, ^/s = 200 GeV. 
As is seen from Fig. |64j the magnitude of the calculated vr*^ spectrum agrees with that of the measured 
one for Pj,<3 GeV. Above this region the shallow tail of vr^'s comes from the jet fragmentation. The 
spectra of decay and thermal photons cross at p^, ~ 4.5 GeV at RHIC. 

Results at the LHC energy are shown in Fig.|65] The spectra of vr^'s and the decay photons are 
shown for two decoupling temperatures. The lower decoupling temperature, Tjec = 120 MeV gives the 
upper boundaries and the higher T^cc = 150 MeV the lower ones. The solid line shows the thermal 
photon spectrum for both the full chemical equilibrium and the off equilibrium cases: they cannot be 
distinguished. This curve is, of course, not the lower limit for the thermal photon emission. For example, 
as is seen from Figs. |59l and |60l thermal production depends very strongly on the assumptions made 
of the initial conditions. In the case of Fig. |65] the thermal and decay photons cross at ~ 5 GeV and 
for Pj,< 2 GeV the amount of decay photons is more than 10 times the number of thermal photons. 
Depending on the amount of vr'^'s from jet fragmentation a window with thermal photons of the order of 
10 % of the background could exist in the region of Pj, above 2-3 GeV. All sources of photons, thermal, 
prompt pQCD from primary interactions between incoming partons, and decays of both the thermal and 
jet hadrons, are compared in the next section. 
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Fig. 64: Spectra of (i) tt^'s (dashed), (ii) decay photons 
from these vr^'s (dotted) and (iii) thermal photons (sohd hne: 
full equilibrium; dash-dotted line: off chemical equilibrium) 
from a hydrodynamical calculation with initial conditions 
from a minijet calculation at ^/s = 200 GeV. Also the vr" 
spectrum, measured by the PHENIX collaboration is shown. 
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Fig. 65: As in Fig. l64l but at = 5.5 TeV. The bands for 
TT^'s and decay photons are shown for Tdec = 120 (upper) 
and 150 (lower) MeV. The predictions for thermal photons 
with full chemical equilibrium and off chemical equilibrium 
cannot be distinguished and are shown by the solid line. 



6.10 Thermal small mass lepton pairs 

P. Aurenche, F. Gelis, H. Niemi, P.V. Ruuskanen and S.S. Rdsdnen 

Like real photons, also virtual photons, decaying to lepton pairs, can be emitted during the thermal 
stage of collisions. As mentioned in sec. 16.221 the dynamics of production of a small mass pair at large 
momentum is similar to that of a real photon and experimentally this channel provides a complementary 
signature to real photon production. 

From Eq. (fTSt and Eq. d24l i one can derive the approximate scaling relation which should hold true 
when the virtual photon mass becomes small compared to all other scales in the problem: 

dN ^ dN ^ 
i+i- _ i+i- 

dtd^y. dM^d^p ~ dtd^x 2EdEd^p 
a dN^ 
~ 7rM2 dtcPxdFp 

where P = {E, p) and = N'P for the lepton pair. 

In Fig.|551the dilepton rate for = 2 and 4 GeV at T = 0.5 GeV is shown as a function of M /T 
and compared to the scaled photon rate at the same energies E. It is seen that for the pair mass up to 
0.5 GeV/c^ the scaling works very well. The scaling holds true also for the emission from the plasma 
integrated over its whole expansion time in the nucleus-nucleus collision. This is shown in Fig.|^for the 
transverse distribution at different values of the pair mass. Even for M = 0.5 GeV/c^ the scaling works 
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well except for the small-pj, region where the ratio M/E becomes large and, as expected, the scaling 
formula overestimates the rate. 
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Fig. 67: Transverse spectra for thermally emitted pairs 
from the plasma phase in a Pb+Pb collision at LHC for 
charged multiplicity dNch/dy = 2900. Photon spec- 
tram from the same hydrodynamical calculation with 
emission rate scaled according to Eg. 1621 



10^ 

lo' 

> 

O 10" 



^10-' 

^ -1 
10 



-a 



10" 



10" 



- 1 1 1 1 1 1 1 1 
l_ Lepton pairs 


thermal, M = 50 MeV 

thermal, M= 150 MeV 

thermal, M = 300 MeV 




Dalitz, M = 50 MeV 

Dalitz, M= 150 MeV 




Dalitz, M = 300 MeV 


--^ 

: ^x^. ■ 










- ^XX 












"n. "1 


: dN,h/dy = 2900 


" \ = 





kj [GeV] 

Fig. 68: Dalitz pairs from rj and tt" compared with ther- 
mal pairs in a Pb+Pb collision at LHC with dA'ch /dy = 
2900. The spectra of ?7's and tt'^'s are obtained from the 
same hydrodynamical calculation as the thermally emit- 
ted pairs. 



For smaU-mass pairs between 0.2 and 0.4 GeV/c^ the dominant background comes from the Dalitz 
decays of ?7's. For M ~ 0.4 — 0.5 GeV/c^ pairs from rj' — > ee'y, lo — > eeir^ and direct vector meson 
decays become comparable and for larger masses (up to 1 GeV/c^) dominant. To get an idea of the 
relative size of the thermal and background contributions we calculate the Dalitz pairs using for vr'^'s and 
ry's the spectra obtained from the hydrodynamic calculation. For the Dalitz decay distribution we use 
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which, when integrated over the pair mass, reproduces well the measured branching ratios P1841. 

In Fig. l68l thermally emitted lepton pairs are compared with Dalitz pairs for three different mass 
values. At the smallest mass value, M = 50 MeV/c^, the Dahtz pairs are dominated by those from vr'' 
decays and they are well above the thermal rate. The situation with the relative size of thermal and decay 
contributions is qualitatively the same as in the case of real photons, see Fig. |65] For the mass values 
above the 7r° mass it becomes more favourable. Even at small transverse momentum, the ratio of Dalitz 
to thermal pairs is at most 5 and the contributions cross at ~ 4 — 5 GeV/c. 

It should be noticed that the discussion here concerns the correlated pairs only. Since the mul- 
tiplicity of neutral pions is expected to be more than 1000, the number of Dalitz pairs from vr*^ decays 
with the 0.012 branching ratio will be ~ 15 per unit rapidity. In the experiment it is not possible to 
identify unambiguously if both particles of an observed e+e" pair came from the same decay or from 
a single thermal emission (correlated pair), or if they came from a different interaction, either decay or 
thermal, and would constitute an uncorrected pair. The number of observed uncorrelated pairs goes 
like ~ A'^pair(-/Vpair — 1) and therefore should rapidly decrease with p^,. An estimate of the uncorrelated 
background was already discussed in sec. l5.6l in the context of the DPMJET model. In principle, the sub- 
traction of the uncorrelated pairs can be done on statistical basis (using like-sign pairs) but the available 
overall intensities may limit the accuracy and consequently the minimum signal-to-background ratio for 
resolving the signal. Since the decays from the same hadrons produce the background both for photons 
and small-mass lepton pairs, the simultaneous measurement of photons and leptons, or any independent 
measurement on vr'^ and r] spectra, can be used to improve the background subtraction. 

7. COMPARING vrO's AND PHOTONS FROM DIFFERENT SOURCES 

P. Aurenche, H. Delagrange and P. V. Ruuskanen 

In this Chapter we present a compilation of the results of Chapters |3 and |3 and we compare the 
strength of different sources of vr '^'s and photons. For the reader's convenience, we briefly summarize 
the features of each of the three models we have considered: 

- Dpmjet which is able to deal with the whole spectrum, down to Pj, = 0, as it contains both soft 
physics (string formation) and pQCD physics in the LO approximation. In A-i-A collisions a high density 
of soft strings is produced and final state interactions are taken into account by fusion and percolation of 
strings. The parameter determining this feature has been fitted to RHIC data. No reference to final state 
energy loss is needed in this approach; 

- Standard NLO QCD calculations which are valid at "high" Pj. only, where "high" is arbitrarily chosen 
to he. p^ > 3 GeV/c: this rather low value for LHC is justified by the fact that DGLAP physics seems 
to hold true down to low x values, and more practically, because down to 3 GeV/c Dpmjet and NLO 
predictions are still in good agreement. The NLO predictions describe the "primary" collisions modified, 
in the case of A-i-A scattering, by shadowing and energy loss; 

- Hydrodynamic evolution of the quark-gluon plasma and hot hadronic matter during which hadrons and 
photons are produced essentially at low transverse momentum. This is refered to as "secondary" particle 
production. 

We naively expect to obtain a full description of a nucleus-nucleus collision by adding the NLO 
contribution to the QGP and hot hadronic matter contributions which describe secondary interactions. 
We therefore have at our disposal two alternative models: Dpmjet and "NLO QCD -i- Hydro" (the 
"standard approach") to make predictions which we are now going to discuss. 

Since the photons from the decays of hadrons, in particular from vr "'s, form a background which 
is larger than the total emission from all other photon sources we begin by comparing the spectra of 
vr ^'s calculated assuming different production mechanisms. In the second part of the Chapter we present 
the compilation and comparison of photon spectra from different sources and discuss the implications 
on the possibility to observe different dynamical features, such as the production mechanism and the 



equilibration, which are beheved to be important in the formation and evolution of matter in lead-lead 
collisions at LHC. A third part briefly deals with lepton pair production. 

7.1 Features of neutral pion spectra 

Hadrons at large transverse momenta originate from hard collisions of incoming energetic partons fol- 
lowed by the formation of hadron jets through the fragmentation of scattered partons. In proton-proton 
collisions this is firmly established and it should also take place in nucleus-nucleus collisions. However, 
as discussed in Chapter|3 nuclear effects like shadowing, energy loss in traversing the system of other 
produced quanta and possible collective effects are expected to be important and to modify the spectra 
from those in hadron-hadron collisions. Most extreme of these modifications would be the thermaliza- 
tion of all produced quanta. However, this is unlikely on theoretical grounds from estimates of energy 
loss and the experimental results at RHIC OJ |3| El clearly show the persistence of the high transverse 
momentum tails in hadron spectra. These measurements also show a clear change in the shape of spectra 
indicating the energy loss of jet-generating partons in the final state. 

On the other hand the very fact of energy loss of high energy partons in the final state points to the 
equilibration of lower energy partons. An interesting question is the characterization of the energy region 
where the transition from partons able to traverse the whole final state to those stopping and participating 
to the collective behaviour of the rest of produced matter takes place. It should be remembered that the 
high energy partons, say Pj,>W GeV/c, leading to the jet formation, carry only a small fraction of the 
energy of all final state particles. This facilitates the hydrodynamic modelling since it should be a good 
approximation to assume that all the energy is equilibrated. 

In the following we summarize the results of the previous chapters by showing those from each 
calculation which are estimated as the most likely. The results for vr'^ production at RHIC are shown in 
Fig. 69. This figure is the same as Fig. |32lto which we have added the Pj, spectrum originating from 
thermal production. At around = 3 GeV/c the rapidly falling thermal spectrum crosses the NLO 
QCD one which is somewhat below the experimental points. The thermal spectrum agrees with the data 
at < 3 GeV/c while the NLO QCD predictions agree with the data above p^, ~ 4 GeV/c where 
thermal production becomes irrelevant. It thus appears that the "NLO QCD -i- Hydro" picture gives a 
very reasonable agreement with data at 200 GeV/c as does DPMJET. These two models are therefore 
good starting points for extrapolation to LHC energies. 

In Fig. 70 we show the results at ^/s = 5.5 TeV. We display: the DPMJET model vr" predictions 
(open circles), the spectra from the NLO perturbative QCD calculation (solid line) and from the hydro- 
dynamical calculation with full thermalization (freeze-out at T = 120 MeV, dotted line and 150 MeV, 
dashed line). The basic input parameters for pQCD calculation (including shadowing and energy loss) 
are summarized in the figure (see Chapter |^for more details and in particular Fig.l33l. The effect of the 
choice of renormalization and factorization scales is shown as a shaded band. 

The results from hydrodynamics dominate the spectrum at low Pj, possibly up to intermediate 
momenta 4 ... 5 GeV/c where they cross below the NLO QCD predictions. So the transverse momentum 
regions where pQCD calculations or the hydrodynamic description apply are complementary. The full 
prediction in the "NLO QCD -i- Hydro" is the sum of the two curves but unfortunately, there is no clear 
way to join the two smoothly. It seems that the difference in slopes is so large that the region of crossing 
is quite stable against the uncertainties in the calculations. At large transverse momenta the DPMJET 
results agree quite well with the pQCD calculation up to the highest values where the models have been 
compared, as can be seen in Fig. |33j this agreement occurs despite the fact that the two models have very 
different mechanisms to describe the interaction of jets with hot matter. 

One can note the rather different shape of the p^, distributions, in the range to 5 Gev/c, between 
the two approaches. Interestingly, the DPMJET results show an even stronger change of slope than 
that from the high-pj, pQCD region and the low-pj, hydrodynamic region. The change also occurs at 
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Fig. 69: Spectra of tt^'s in Au+Au collisions at = 200 GeV shown for the following calculations: (i) pQCD jet produc- 
tion and fragmentation (solid line), (ii) DPMJET model (open circles), and (iii) pQCD+saturation for production followed by 
hydrodynamics with decoupling at Tdoc ~ 150 (dashed line). 



lower momenta, around 1. . . 2 GeV/c, and the small-p^ part is much steeper than that obtained from the 
hydrodynamic calculation. If the final multiplicity (the total particle multiplicity from the hydrodynamic 
calculation ~ 4500) is as high as or close to what follows from the pQCD + saturation model initial 
conditions, it seems difficult to obtain such a steep slope from a hydrodynamical calculation. The reason 
for this is that it takes a long time to dilute the dense initial energy of a high-multiplicity collision; a 
strong transverse collective flow has time to build up and boost the particles to high transverse momenta 
even in case of relatively low thermal momenta at freeze-out. These differences show that the shape 
of the transverse momentum spectra of hadrons, even at relatively low transverse momenta, contains 
important information on the dynamics in the nucleus-nucleus collision. These differences in shape are 
much more marked at LHC than at RHIC. 



7.2 Comparing photon spectra from main sources 

Photon spectra from different sources are shown in Figs. 71 and 72 for RHIC and LHC, respectively. 
The figures are more complicated than in the case of pions since photons can be produced either directly 
or through the decay of hadrons. 

At RHIC we see that around ~ 3.5 to 4 GeV/c the NLO QCD production of prompt photons 
(solid line) is of the same order of magnitude as the production of thermal photons (dashed-dotted lines) 
and the decay photons from hadrons as obtained in the hydrodynamic calculation with tq = 0.19 fm/c 
(thick dashed line). Thermally emitted photons decrease rapidly with Pj, and compared to the prompt 
QCD photons and photons from hadronic decays they become negligible already around 6 GeV/c. There 
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Fig. 70: Spectra of tt'^'s in Pb+Pb collisions at ^/s = 5.5 TeV. The following calculations are shown: (i) pQCD jet production 
and fragmentation (solid line with the grey band), (ii) Dpmjet model (open circles), and (iii) pQCD+saturation for production 
followed by hydrodynamics with decoupling at Tdec = 120 (dotted line) and 150 MeV (dashed line). 



is no sensitivity to whether or not the plasma is in chemical equilibrium except at the higher values of p^, 
where the thermal rate is negligibly small. One can also observe an overall agreement between DPMJET 
and the NLO QCD + hydrodynamic estimates over the whole Pj, spectrum. 

We next discuss the LHC predictions shown in Fig. 72. The prompt NLO QCD photon spectrum is 
given by the lower solid line and the spectrum of decay photons from vr^^'s, produced through NLO QCD 
jet fragmentation, by the upper solid line. The dotted and dashed lines show the decay photons from 
all hadrons in the hydrodynamic calculation with Tdec either 120 MeV (dotted) or 150 MeV (dashed). 
Thermal photons are shown both from a calculation when the quark-gluon plasma is assumed to be in 
full thermal and chemical equilibrium (dashed-dotted line with filled triangles) and when the plasma is 
dominated by gluons and is not in chemical equilibrium (dashed-dotted line with open triangles). As 
mentioned before kinetic equilibrium is assumed also in the latter case and in both cases the initial time 
in the hydrodynamic calculation is tq = 0.1 fm/c. The two sets of predictions cannot be distinguished. 
Given the same initial energy density from the primary production, the deficit of quarks and antiquarks, 
which reduces the photon emission, leads to higher values of the temperature which enhance the emis- 
sion. The slope of thermal production is steeper than for the prompt QCD photons from the interactions 
of incoming partons but the crossing of thermal and prompt QCD photon spectra takes place at rather 
high value of p^ around 6. . .7 GeV/c. However, if the effective thermalization time is longer the crossing 
region could be as low as 3 ... 4 GeV/c. 

In order to compare more clearly the different sources, we have plotted in Figs. 73 and 74, for 
RHIC and LHC, respectively, the ratio of the spectrum of "all" photons to the background (from hadronic 
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Fig. 71: Photon spectra in Au+Au collisions at ^/s = 200 GeV. Thermal emission is calculated assuming either full thermal 
equilibrium in the plasma phase (dashed-dotted line with dark triangle) or a gluon-rich plasma (dashed-dotted lines with open 
triangle). Decay photons from thermal hadrons (thick dashed line) and decay photons from NLO QCD pions (thin dashed line) 
are also shown together with the prompt photons (solid line) and decay photons from the DPMJET calculation (open circles). 
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Fig. 72: Photon spectra in Pb+Pb collisions at ^/s = 5.5 TeV. Thermal emission is calculated assuming either full thermal 
equilibrium in the plasma phase (dashed-dotted line with dark triangle) or a gluon-rich plasma (dashed-dotted lines with open 
triangle). Decay photons from thermal hadrons are indicated by dotted or dashed lines depending of the decoupling temperature. 
The spectrum of NLO prompt and decay photons are shown by the lower and upper solid lines respectively. 
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Fig. 73: Ratio of all (NLO+Hydro) photons over all decay (NLO QCD + Hydro) photons as a function of p^, in Au+Au 
collisions at ^/s — 200 GeV. Predictions for primary production of prompt photons (NLO QCD) are shown by the dotted line. 
Preliminary data from PHENIX are also displayed. 




Fig. 74: Ratio of all (NLO+Hydro) photons over all decay (NLO QCD + Hydro) photons as a function of p^, in Pb+Pb 
collisions at ^/s = 5.5 TeV. Predictions for primary production of prompt photons (NLO QCD) are shown by the dotted line. 



decays) photon spectrum. 

In Au+Au collisions at ^/s = 200 GeV, the thermal photons from the hydrodynamic calculation 
(solid line) provide a small increase in the ratio 7Aii/7Bkgd in the range 3 GeV/c to 6 GeV/c compared to 
a pure pQCD estimate (dotted line). However this increase appears too small to be observed experimen- 
tally. The thermal results are insensitive to whether the plasma is in chemical equilibrium or not. With 
increasing , due to prompt QCD photons, the fraction of direct photons grows quite rapidly reaching 
50 % already at 10 GeV/c. Comparison with PHENIX preliminary data ISTI shows a good agreement 
for Pj, < 4- GeV/c (where thermal effects play a major role) while the theoretical predictions tend to un- 
derestimate the ratio at larger p^, values where the NLO QCD production mechanism is dominant. This 
disagreement can be related to an overestimate, at large p^, of vr*^ rates, and hence of decay photons, 
compared to data (see Fig. 69). 

In Fig. 74 the ratio 7Aii/7Bkgd is shown for Pb-i-Pb collisions at LHC. The upper curve is ob- 
tained by using the prediction from the plasma in equilibrium for the thermal photons and 150 MeV 
for Tdoc which, for the calculation of the decay photons, gives the steeper slope for the hadron spectra. 
In the case of the lower curve off-chemical equilibrium is assumed in calculating the thermal photons 
and Tdec = 120 MeV in obtaining the hadron spectra for the calculation of the decay photons. The 
dash-dotted line indicates the prediction obtained under the hypothesis of a gluon rich plasma and a de- 
coupling temperature of 150 MeV consistent with the RHIC data (see Fig. 68). Note also the rather flat 
dependence, at around the 10 to 15% level, of the ratio for p^, between 3 and 10 GeV/c where thermal 
effects play a role. Above 10 GeV/c the ratio increases rapidly with p^ so that direct photon production 
represents 25% at p^ = 20 GeV/c and 60% at p^ = 60 GeV/c. 

In conclusion, the overall behaviour of 7Aii/7Bkgd is similar with that in the RHIC results but 
the ratios differ in details, in particular the enhancement due to thermal production is shifted up to the 
interval from 3 to 10 GeV/c and is larger. 

We should add to the rate of photons produced in A-i-A collisions also the possibility that before a 
hard final state parton fragments into a jet of hadrons it may convert to a photon in the parton medium via 
Compton scattering or annihilation as proposed in [1811. The resulting spectrum should have a power 
behavior and could dominate over the thermal photons already for p^, ~ 6 ... 9 GeV/c (see Fig. l63l . 
This mechanism is hard to quantify, however, and the presently available calculation has been done 
assuming plasma evolution without transverse expansion which is not compatible with the one used 
in our discussion. Therefore we do not attempt a quantitative comparison. Nevertheless, it would be 
important to further develop the modelling of this process. 

There are other uncertainties in these results for 7Aii/7Bkgd> the most serious for the thermal con- 
tribution being the uncertainty in tq, the initial time in the hydrodynamical calculation when the matter 
is assumed to reach an (approximate) kinetic equilibrium. If the time for the beginning of the production 
of thermal photons is more than 1 fm/c, possibilities for observing them are not very encouraging. How- 
ever, one should not be too discouraged by those estimates which predict kinetic thermalization times 
over 1 fm/c. It may take time for the initially produced parton system to get equilibrated, but even then 
we can expect it to be a dense system of interacting partons right after their production. When interact- 
ing (and evolving towards equilibrium), these partons would emit photons. The uncertainty then is if the 
hydrodynamical calculation with the assumption of kinetic equilibration gives a reasonable estimate of 
the emission of photons from a nonequilibrium but dense parton matter. 

If the production of photons from secondary collisions is large enough to be observed, then even 
if the produced parton system never reaches equilibrium, these photons carry information on the typical 
scale of relative momenta in the parton matter and on its density. In an equilibrated system the momen- 
tum scale is given by the temperature, and the density and the momentum scale are related. If early 
equilibration does not take place the observation of photons from secondary collisions would still be a 
direct window into the properties of the densest phase of partonic matter formed in nuclear collisions. 



As the discussion above shows, the phenomenology of photon production is far from being on a 
quantitative basis yet. However, our studies seem to indicate that, at LHC, the ratio of direct photons 
over decay photons, whatever their origin, should be 10% or more for Pj, > 3 to 4 GeV/c. At the lowest 
values of the direct photons would be mostly produced in secondary collisions while around 10 Gev/c 
they would be of mixed origin with the perturbatively produced (/. e. in primary collisions) photons 
becoming dominant with increasing p . Many of the uncertainties concerning thermal production would 
be reduced if the model parameters could be fixed from hadron data leaving then the photon spectrum as 
a crucial consistency check of the whole approach. 

7.3 Comparing lepton pair spectra from main sources 

In this section we briefly compare lepton pairs from NLO QCD interactions, from hadron decays, and 
from thermal interactions in the plasma phase. We present the transverse momentum distributions of pairs 
in the mass interval [0.2,0.6] GeV/c^ and consider correlated pairs only. The results, discussed separately 
for each source in previous sections, are summarized in Fig. 75. In all models and calculations the same 
input is used as for the photon summary in the previous section. 



O 10 



■O10 

a 
■a 

CLIO" 
CM 



pairs - \|s^ = 5.5 TeV Pb + Pb 5 % Most Central Collisions 



10 ^ 



10 



10 ^ 



10 = 



DPMJET : J 
■ 

Hydro : 



1 ' \ '~ 
Ranft et al. 
Correlated pairs only 
H. Niemi et al . 
.. 0.2<Mg.g <0.6GeV/c 
.. 0.2<M^.%0.4GeV/c 
. 0.4<m'/ <0.6GeV/c 



i~ NLO pQCD : EKS98 + ELoss - All scales 1 x p . 

CTEQ5M + BFG set II 



0.2<M .^.<0.6 GeV/c 



5 6 
Pj (GeV/c) 



Fig. 75: Lepton pairs from thermal emission are compared witii those from prompt emission in the NLO QCD approach. 



The main background coming from hadron decays is illustrated with DPMJET results (filled 
squares) since the DMPJET calculation agrees with the NLO QCD results for hadrons at large and with 
the hydrodynamic calculation at low transverse momentum. Results from the NLO QCD calculation on 
pair production are shown as the solid line and thermal pairs from the hydrodynamic calculation as a 
dashed line. Contributions of thermal pairs from each half of the mass interval are shown also separately 
to draw attention to the rather strong mass dependence in the thermal production. 

Statistics of the DPMJET calculation is not large enough for obtaining results beyond p^, ~ 
3 ... 4 GeV/c but the indications of the low-p^ behaviour seem clear: pairs from hadron decays domi- 
nate at low Pj, but they have a steeper p^, dependence and seem to join the distribution of thermal pairs 



around ~ 3 GeV/c. Pairs from hard interactions of incoming partons, given by the NLO QCD result, 
overtake the thermal production around ~ 6 GeV/c. If the dependence of different contributions 
is similar to that in the photon production, we would conclude that for Pj, > 3 GeV/c the background 
from the hadron decay pairs has the same size as the sum of thermal and NLO QCD pairs and might go 
below it as Pj, increases. 

In summary, considering the signal-to-background ratio the e+e^ channel appears at first sight 
quite promising for the study of hot matter in nucleus-nucleus collisions. However, our discussion is 
preliminary and lots of work remains to be done in order to confirm if the spectrum of direct lepton pairs 
can be extracted from the large background of uncorrelated pairs within the statistics which could be 
expected to be available at LHC. 

8. PRODUCTION OF PHOTONS AT "LARGE" TRANSVERSE MOMENTUM AND CORRE- 
LATION STUDIES 

F. Arleo, P. Aurenche, J. -Ph. Guillet,O.L. Kodolova, LP. Lokhtin, A. Nikitenko, I.N. Vardanyan, M. Werlen 

It was shown in previous chapters that the signal for direct, thermal or non-thermal, photon pro- 
duction in A-i-A collisions is not too large compared to various backgrounds, making it a challenging but 
not easy task to use inclusive photon production as a signal for quark-gluon plasma formation in heavy 
ion collisions. In this chapter we consider correlation observables in a kinematical regime where the pho- 
ton is produced "promptly" in a hard QCD process and we study the hadrons recoiling from the photon. 
These hadrons are decay products of a jet, most probably a quark. Since the fragmentation properties of 
jets are supposed to be modified by medium effects, differences in the shape of the correlations should 
be seen when comparing pp and A-i-A collisions. In the next section one looks at the global jet and study, 
from an experimental (CMS) point of view, the effect of smearing and efficiency on an assumed energy 
loss. Then a discussion is given on photon-photon and photon-hadron correlations where the hadron or 
the second photon are decay products of the jet. One takes a more phenomenological point of view and 
discuss various functions in the hope of determining the correlation most sensitive to medium effects. 
All studies are preliminary in the sense that they are conducted at the leading order and therefore will be 
modified by NLO corrections which may be very important in some comers of phase space. 

8.1 Photon-jet correlation at CMS 

O.L. Kodolova, LP. Lokhtin, A. Nikitenko, L. Vardanian 

Among other proposed signals of jet quenching (see Chapter "Jets" of this Report l44l '). the p^- 
imbalance between a produced jet with a gauge boson in 7+jet ll46l 11851 and Z+jet 1X86 1 production 
has been identified as being observable in heavy ion collisions using the CMS detector 1 16|. The dom- 
inant leading order diagrams for high transverse momentum 7 -1- jet production are shown in figure ^61 
Contrary to the gluon-dominated jet pair production where one could investigate jet quenching due to 
mostly gluon energy loss in dense matter, 7+jet channel gives a possibility to study quark energy loss. 
The main background here is hard jet pair production when one of the jet in an event is misidentified as 
a photon. The leading vr*^ in the jet is a main source of the misidentification. Table l8. 1 I presents the event 
rates for signal and background processes in one month of Pb-i-Pb beams (a half of the time is supposed 
to be devoted to data taking), i? = 1.2 x 10^ s, assuming luminosity L = 5 x 10^^ cm^^s^^ to that 

N{events) = Ra'Xy^L, 

where production cross sections in minimum bias nucleus-nucleus collisions is obtained from those in 
pp interactions at the same energy {^/s = 5.5 TeV) using simple parameterization cr^^ = A'^a'i . The 
cross section in pp collisions were evaluated using the PYTHIA_6.1 Monte-Carlo generator l62l with 
the CTEQ5L parton distribution function. Note that the influence of nuclear shadowing is practically 
negligible for the region of sufficiently hard 7+jet production, xi.2 ~ a/?/s > 0.2. However, large 



theoretical uncertainties in absolute rates in pp collisions come from choice of the parton distribution 
functions, next-to-leading corrections, etc. It means that measurements in pp or D-i-D collisions at the 
same or similar energies per nucleon as in the heavy ion runs are strongly desirable to determine the 
basehne rate precisely. 



Table 4: Expected rates for 7+jet and 7ro(— > 27)+jet channels in one month of Pb— Pb beams. 



Channel 


Barrel, r/ < 1.5 


Barrel-i-Endcap < 3 


7-f-jet, E^r^""^'^ > 100 GeV 


1.6x10^ 


3.0 xlO^ 


7ro(^ 2j)+}et, Ei^^'^ > 100 GeV 


8.4x10^ 


2.2 xlO'' 



One can see that even for events with Ej' > 100 GeV the background is still dominant. Signal-to- 
background ratio becames close to 1 only above 200 GeV. The identification of the influence of the dense 
medium formation on signal spectra requires the reduction of the background. One of the possibilities 
is to apply some kind of the photon isolation, so called "zero suppression criteria", which requires no 
energy above a given threshold around photon (see section "Photon detection at CMS" for details). In 
this case the signal-to-background ratio at Et > 100 GeV can be improved by a factor about 2.3 at the 
cost of a 14% reduction of the signal. Another possibility is to apply some kinematical cuts, which do not 
have influence on the -imbalance of the process (i.e. not result in shift of maximal value of Ej, — E^^ 
distribution). 

The possibility to observe the medium-induced energy loss of quark-initiated jet using photon-jet 
correlation in heavy ion collisions with CMS detector has been investigated in 1187 16 1. It has been 
found that initial state gluon radiation and finite jet energy resolution (which is much larger than pho- 
ton energy resolution) result in significant smearing of distribution of differences in transverse energies 
between the photon and jet. But it is still symmetric: E^ = E?^^ only in average (not for each given 
event). The non-symmetiic shape of the distribution appears if a jet loses energy: the maximal value of 
the distribution is equal to the average energy loss of the quark-initiated jet at the given energy detection 
threshold, E^r^^ ~ 100 GeV in CMS case. Note that we are not measuring energy loss of a leading quark 
by such a way, but getting total loss of quark-initiated jet outside the given jet cone. 
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Fig. 76: Leading order diagrams for 7+jet (signal) and 7ro(^ 27)+jet (background) production. 




Fig. 77: The distributions of differences in transverse energy between the 7 and jet with E^' > 120 GeV: a) 
without (7r° + jet) background, b) with (7r° + jet) background. The pseudorapidity coverage is | ij^^ jet \< 1-5. 
Different values of jet in-medium energy loss, initial state gluon radiation and finite jet energy resolution are taken 
into account. 



In order to test the sensitivity of 7+jet production to jet quenching, it was considered three 
scenarios with average collisional energy loss of a jet: < AEg >~ 0, 4 and 8 GeV respectively 
(< AEg >= 9/4- < AEg >). The jet energy resolution at mid-rapidity obtained for Pb+Pb colli- 
sions has also been used to smear the energy of the recoiling parton, as well as the jet rejection factor and 
signal efficiency 1 187 1. Figure^Jshows the distributions of differences in transverse energy between the 
photon and jet with E^' > 120 GeV in one month of Pb-i-Pb beams without and with (vr'^ -1- jet) back- 
ground in the pseudorapidity region | rj"'' |< 1.5 for different values of jet energy loss. In this case 
a luminosity L = 10^^ cm~^s~^ was assumed and PYTHIA_5.7 version with the default CTEQ2L pdf 
choice was used. The jet energy resolution leads to a difference between the input values < AEq > and 
the ones obtained from the spectra. The background of vr*^ -contamination results in non-zero negative 
values of the final distributions (figure ITTb ) already in the case without jet energy loss. However one can 
see that the shape of the distribution is well distinguished for the scenarios considered. For the region of 
{E^-E^r^') > there is a difference for almost every bin greater than 1 standard deviation for the rather 
small jet energy loss 8 GeV and even for the loss 4 GeV. In the real experiment it would be possible to es- 
timate the number of background events using the region without the signal (Ej, — Ei^^) < — 100 GeV/c 
and background shape from Monte-Carlo simulation and/or from pp data. A significant difference in the 
shape of the Ej, — Ej^^ distribution can allow to optimize extraction of the signal from the experimental 
spectra. 

8.2 Photon-hadron and photon-photon correlations in pp and A+A collisions 

E Arleo, P. Aurenche, J. -Ph. Guillet 

In the previous chapters it was suggested that thermal production of photons could be seen for 
rather low p^, values. In the following we consider, on the contrary, the case of a photon with a high 
enough p^ so that it is produced promptly in a hard QCD process and therefore not affected by thermal 
effects and we study the decay products of the jet recoiling from this photon. The analogue of the 
inclusive spectrum Eq. Q, at the leading logarithmic order, can be written II89I 

da^^ CD 1 dz^ T,^\f^ dz^ , \ , 

F ,Axi,M) E^,(x2,M) _^ 

where the p^. (resp. k^,-) are the final state particle (resp. partonic) transverse momenta and the scaling 
variables Zi are defined by Zi = Pj.Jkj.-. The quantity |-^^lafe^ccZ '■'^^ matrix element squared, aver- 
aged over spin and color, of the partonic sub-process ah cd. Higher order corrections to Eq. (l64t 
have been calculated |189J and, based on existing data L190lll9li ri92l. extensive NLO phenomenolog- 
ical studies of correlations have been carried from fixed target energies to the LHC ll83l 11931 . If one 
looks at a photon at large enough Pj,^,it will be directly produced and the corresponding fragmentation 
will reduce to 8 {I — z^). Studying various correlation variables will then allow the mapping of the 
fragmentation function D £) j ^{z^i, M p) since all other ingredients (structure functions, matrix elements) 
are known. One could then observe the effect of partonic energy loss in the medium by comparing 
correlation functions in A-i-A collisions with the equivalent ones in pp collisions. In the following we 
restrict ourselves to the case of promptly produced photon and do not impose isolation cuts. We do 
not consider the "background" contribution generated by photons decaying from hadronic resonances, a 
problem presently under investigation 11941 . 

How to express final state interactions through medium-modified fragmentation functions remains 
however far from being clear. To illustrate the effect of parton energy loss in dense media on correlation 
functions, we shall adopt here the effective model suggested in Ref. LI 85 1 . Within this approach, the 




Fig. 78: Comparison of partonic fragmentation functions in tiie vacuum and in the medium for typical values of the energy 
loss parameter uic- Top two figures: k^^ = 20 GeV/c up quark and gluon into a photon; bottom two figures: k^^ = 20 GeV/c 
up quark and gluon into a charged hadron. 



parton energy shift leads to a rescaling of the momentum fraction Zd in presence of a QCD medium, 



_ Prd * _ Prd _ Zd 

"'Td '^Td ^ *=/ '^Td 

where p^,^ (respectively, A;^^) stands for the transverse momentum of the photon or hadron (respectively, 
the parton), and e the energy lost by the hard parton while going through the medium. Consequently, the 
medium-modified fragmentation functions D^'j^{zd, Q^, k^^) may simply be expressed as a function of 
the standard (vacuum) fragmentation functions Dj^/dizd, Q"^) through 11851 

Zd D'Sfdizd, Q^ fc^J = de Vie, fc^J z*^ Dn/dizlQ^). (66) 

Here, V{e, k^^) denotes the probability for a parton with transverse energy E = k^^ to lose an energy 
e L188J . Note that medium-induced fragmentation functions now depend explicitely on the parton trans- 
verse energy krj, . Assuming the soft gluons radiated by the leading hard parton to be emitted indepen- 
dently, this probability distribution (or quenching weight) can easily be related to the medium-induced 
gluon spectrum determined perturbatively, dl/dto, characterized by the energy scale Ii45i 

u^c = ^qL^ (67) 

The gluon transport coefficient q reflects the medium gluon density while L is the length of matter cov- 
ered by the hard parton. Note that q can be very large in practice in a hot pion gas or quark gluon 
plasma 11951 . Based on the estimate from the pion inclusive spectrum in A-i-A collisions at RHIC 
(y^ = 200 GeV), we shall illustrate our results using the following values w c = 30 and 50 GeV for 
the characteristic scale of the dense medium produced in nuclear collisions at the LHC. The distribution 
V{e, k^^) has been given a simple analytic parameterization in L49J which we shall use in the present 
calculations. 

Since fragmentation functions fall steeply with z, even a small shift Az^ = — « Zde/k^^ in 
(l65l may substantially affect fragmentation processes in presence of a hot and dense QCD medium. This 
can be seen for instance in Figure ^8] where fragmentation functions into a photon {top) or a charged 
hadron {bottom) are computed for k^ = 20 GeV/c up quark {left) and gluon {right) traversing the 
medium. In particular, we expect the effects of parton energy loss to be more pronounced as z gets 
larger. This will be further discussed below when correlation functions in nuclear collisions are com- 
puted perturbatively. 

As an example we consider the following kinematical constraints. Both observed particles (photon 
or charged hadron) are required to be in a rapidity interval [—0.5, +0.5] and furthermore one photon 
should have a transverse momentum greater than 10 GeV/c while the other photon or the charged hadron 
is required to have p^, > 3 GeV/c. Integrated over all kinematical variables, with the above constraints, 
we display in Fig. ^9] the photon spectrum da / dpj. for the case of two photon events {left) and the 
charged hadron spectrum for the photon -i- hadron events {right). The various dynamical components are 
shown. In the two photon case they are the "direct" (both photons produced in the hard sub-process), the 
"one-f (one photon produced directly and the other by bremsstrahlung from a fragmenting parton) and 
the "two-f (both photons produced by bremsstrahlung of the final state partons) The interpretation 
of the photon-photon figure is easy. For < 10 GeV/c, the observed photon is the less energetic of 
the two photons because of the 10 GeV/c cut imposed on the other photon: the observed photon can 
only be produced by bremsstrahlung and the figure shows that the "one-f process dominates by far. 

^'We follow for the labeling of the various terms the conventions used in the DIPHOX code 1 196 1. In the definition of 
Chap. 2 a photon produced "directly in the hard sub-process" or a photon produced by bremsstrahlung of a parton produced in 
a hard QCD sub-process are both "prompt" photons since they emerge from a primary collision and are not decay products of 
a resonance. 




Fig. 79: Details of the various dynamical components in pp collisions at ^/s — 5.5 TeV. Left: the photon spectrum in the two 
photon events; right: the charged hadron spectrum in the photon-hadron events. 



When Pj, > 10 GeV/c, the direct component comes into play but again the "one-f process is the largest. 
This confirms that, with the cuts used, the cross section Eq. d64t for the case of two photon production 
is mainly sensitive to one fragmentation function. For the photon-hadron case the "one-f (the photon 
produced directly and the hadron a fragment of a jet) and the "two-f (both particle fragments of jets) 
processes are comparable. A word of caution is necessary: as usual in perturbative QCD, the distinction 
between direct, one-f and two-f is only indicative as the relative weights of the components depends on 
the choice of the (factorization, renormalization) scales, chosen in this discussion to be all equal to p^, . 

In the following we define various correlation functions in 7 — 7 events II196I (see Fig. l80l and 
study their modifications as a function of the assumed energy loss characterized by the value of the 
parameter lUc- They are : 

- the photon Pj, distribution; 

- the diphoton invariant mass m^^ spectrum; 

- the spectrum in the transverse momentum of the pair defined as = p^^ + p^,^ ; 

- the distribution in the scaled momentum fraction z 34 = — p^g • Pta/Pts- 

Considering the spectrum da/dp^ in Fig. [SOlone sees that the curve below p^ = 10 GeV/c 
decreases when Uc increases (energy loss increases), very similarly to the fragmentation functions in 
Fig. ^8] The relative suppression is larger as p^, increases corresponding to an increase of the effective 
fragmentation variable z. For p^ above the cut, on the contrary, the medium and the vacuum spectra 
come closer together when the transverse momentum increases. This can be understood as follows: the 
detected photon with large momentum is produced directly while the other photon is produced predom- 
inantly near the 3 GeV/c threshold and its effective z ~ z 34 decreases when the p^. of the directly 
produced photon increases: since at small z the energy loss in the medium is small (see Fig. l78l the 
in-medium curve and the "vacuum" curve approach each other. 

The di-photon invariant mass spectrum (upper right plot in Fig. l80l shows, in vacuum as well as 
in the medium, a characteristic bump above m-y-y = 20 GeV/c due to the onset of the "direct" production 
mechanism. Only "one-f and "two-f processes contribute below this threshold. When increasing the 
energy loss, these last two processes are suppressed, leading to a larger suppression of the invariant mass 
distribution at low masses and an increasingly marked shoulder since the direct channel is not affected 
by the medium. 

The peaked shape of the transverse momentum distribution of the pair (lower left plot in Fig.lSOl 
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Fig. 80: Photon-photon correlation functions in Pb+Pb collisions at ^/s = 5.5 TeV. 
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Fig. 81: Ratio of photon-photon correlation functions in Pb+Pb collisions and pp collisions at = 5.5 TeV. 



at Pj, 1 GeV/c is an effect of the imposed asymmetric lower cuts on the photon momenta and the 
decrease of the spectrum at the lower q,^ values when the energy loss increases is easy to understand: to 
obtain a small value both photons should have similar p^. which means that the bremsstrahlung photon 
will be emitted from a parton with a large z value, a configuration highly unfavored in the medium. The 
suppression of the spectrum for small g-j, = will be somewhat modified by higher order corrections 
which will smear the direct process contribution shown, in the figure, by a point at g^, = 0. On the other 
hand at large q^, no dramatic medium effect is expected since the dominant kinematics will be a photon 
(directly produced) at large p^ and a bremsstrahlung photon at very small p^, equivalently small z, both 
of which are not affected by medium effects 

Finally, the distribution in the scaling variable z 34 (bottom right plot) nicely reflects the behavior 
of the fragmentation function when medium effects are increasing: the scaling variable and the fragmen- 
tation variable z are very closely related and it is natural that the z 34 be reminiscent of the behavior of 
the fragmentation functions shown in Fig. ^8] when the energy loss increases. 
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Fig. 82: Photon-charged hadron correlation functions in Pb+Pb collisions at ^/s = 5.5 TeV. 



The energy loss features of the observables we just discussed are emphasized if we plot the ratio 
of the correlation functions in A-i-A collisions over the same in pp collisions. This is done in Fig.[^ The 



effect of the energy loss does not reduce simply to a change of normalisation of the correlations but also 
to a change of shape. For the normalized da/dp^, distribution {top left) one recognizes the increase in 
the suppression as p^ increases below the threshold value of 10 GeV/c, corresponding to an increasing 
z value of the detected photon while above p^ = 10 GeV/c the effect of the energy loss is, as expected, 
becoming less and less visible. Particularly interesting are also the normalized da/dq ^ {bottom left) 
and da / dz^i {bottom right) correlations which display a marked dependence on the variable q ^ and 2:34 
respectively. 

All these studies have been carried out in the LO approximation. Clearly, higher order corrections 
will somewhat smear the results specially near thresholds or infra-red singular points (z = 1 or g ^ = 
for example). Very preliminary studies in the NLO approximation, assuming that higher order corrections 
can be calculated in the same way in pp and in A+A collisions, show that the effects of the energy loss 
can still be clearly seen in the correlation functions 11941 . 

As a further illustration of correlation studies we present in Fig.|82lthe results for photon-charged 
hadron observables, with the same kinematical constraints as before: p^. > 10 GeV/c for the photon and 
p^ >3 GeV/c for the hadron. The interpretation of the curves is very similar to the case of photon-photon 
correlations the main advantage of the photon-hadron channel being a much higher counting rate. 

9. OUT-OF-EQUILIBRIUM PHOTON PRODUCTION 

/. Dadic, F. Gelis, G.D. Moore, J. Serreau 

In the previous sections, l^and^ the rate of thermal photon production in the quark-gluon plasma 
is calculated ignoring the finite life time of the plasma. Recently however it has been argued that finite 
life-time effects could considerably increase the rate of production. We discuss this class of models and 
show that they are based on an unrealistic modeling of the initial and final state conditions which leads to 
an infinite amount of energy being radiated away by the photons. These models are therefore unphysical. 
A possible approach to non-equihbrium aspects is then briefly sketched. 

9.1 Transient effects in the real-time approach 

/. Dadic, F. Gelis, G.D. Moore 

In equilibrium, one usually defines the initial density operator at an initial time t , which is infinitely 
remote in the past. However, one may wonder if any changes are to be expected when the initial state is 
defined at some finite time tj. Some recent papers by Wang, Boyanovsky and Ng 1 197. 198 1 suggested 
that one can see transient effects in the photon radiation of a thermal system that exists only for a finite 
amount of time. A similar model has been considered by Dadic, based on [199 200 201T, but there 
are no published results on photon production at the time of writing this report. The main feature of the 
photon emission spectrum obtained in [ 197 198 1 is that it decreases as a power of the photon energy 
instead of the usual exponential spectrum. This implies that already at moderate photon energies, the 
prediction of II197I I198II dominates by several orders of magnitude over the usual thermal rates. 

We start this section by discussing some general properties of an equilibrated system prepared at 
a finite initial time, and discuss what to expect if the particles in the initial bath are non interacting. We 
proceed with a derivation of the results of [197 . 1981 for the photon emission rate in such a model, and we 
discuss also briefly the approach of 11991 12001 120II . Finally, we close the section by discussing the fact 
that this model leads to infinite "vacuum" contributions, precluding any phenomenological application 
of these rates. 



9.11 Starting the equilibrium evolution at a finite time 

Let us first assume that the density operator at ti is p{ti) = ex.p{—(]H), where H is the full Hamiltonian, 
including all the interactions. In order to calculate perturbatively a correlator like (T0(xi) • • • (/){xn)) 
where {A) = Tr {p{ti) A) /Tr {p{ti)), one has to extract the two sources of dependence on the coupling 
constants. The first one is the usual dependence of the Heisenberg field 4>{x) on the coupling constants, 
and it is readily extracted by going to the interaction picture: 

= U {ti,xo)(t>j {x)U(xo,ti) , (68) 

where U is an evolution operator defined as 

U{t2,h) = P expi / ' d^x£int(</',(x)) (69) 
Jti 

with £int the part of the Lagrangian density that contains the interactions and 0^ a free field that coincides 
with the Heisenberg field at xq = ti. 

However, there is a second source of dependence on the coupling constants, which is present in 
the density operator itself |202 203 1. This can be extracted thanks to the following formula that relates 
the full density operator to the one defined with the non-interacting Hamiltonian Hq [2041: 

g-/3H ^ g-^-^op expi / d^x£int(0,(2;)) . (70) 
Jti 

Eqs. d68l . (l69l and (f/Ol explain the structure of the time integration contour in the real-time formalism 
1205 1 1206i i207 1 1208 1 1209 1 : it starts on the real axis at ti and goes along the real axis up to some final time 
tf (which can be +oo, but needs in fact only to be larger than any physical time in the problem, since 
causality forbids any dependence of physical quantities upon this time), then goes back to ti along the 
real axis, and finally goes down to ti — ip following a segment parallel to the imaginary axis. 

The bare propagators in this formalism obey the so-called Kubo-Martin-Schwinger (KMS) sym- 
metry l2T0ll2TTl . which reads 

Go{ti,yo) = aGo{ti-il3,yo) , (71) 

where q = 1 for a real scalar field, a = exp{/3fiq) for a boson that carries the conserved charge q with 
an associated chemical potential p, and a = — ex.p{f3fiq) for a fermion carrying the conserved charge 
q. The same symmetry holds for the second point of the propagator, with q changed into —q. Because 
the number of fermion lines arriving at a vertex is even and because the total conserved charge arriving 
at a vertex is zero, it is easy to verify that any Feynman diagram made of these bare propagators and 
where the time integrations at the vertices run on the previously defined contour is independent of the 
initial time ti 1202112121 . A corollary of this property is that all the Feynman diagrams depend only on 
differences of their external times, even if they have been calculated on a contour that breaks invariance 
under time translation because of the finite ti. 

Physically, this result means the following: if one prepares a system in a state of thermal equilib- 
rium at some finite time ti, all the Green 's functions of the theory are independent of ti and there is no 
way to tell what the initial time was. In other words, there are no observable transient effects in a system 
initially prepared in a state of thermal equilibrium (provided the particles are interacting in the initial 
statistical ensemble). 

9.12 Non interacting initial statistical ensemble 

In the previous argument, the contribution of the vertical branch of the contour is crucial in order to 
achieve a result that does not depend on ti. Indeed, the KMS symmetry that relates the values of the 
propagators at both ends of the contour is of no help if we terminate the contour at ti. If we perform the 



time integrations at the vertices using only the two horizontal branches of the contour (i.e. the so-called 
Keldysh contour), then we lose the property of independence on ti, as well as the invariance under time 
translations of the Green's functions, and there are now transient effects. 

Given that the vertical branch of the time path arises from Eq. dTOb . not using the vertical part is 
equivalent to using an initial density operator p{ti) = exp{—(3HQ) which does not contain the interaction 
terms of the Hamiltonian. In other words, this is equivalent to having initially a statistical ensemble of 
non interacting particles. An identical situation would be the following: start at t = — oo with the density 
operator p{—oo) = exp(— /?ffo) and assume that all the coupling constants are zero before ti and jump to 
their normal value at tj. Indeed, if we start with p{—oo) = exp(— /3i/o) and if we have a non interacting 
dynamics between — oo and ti, then we will have p{t) = exp(— /?ffo) as long as the interactions remain 
zero, by virtue of p{t) = —i[HQ, p{t)] = 0. 



9. 13 Photon production 

Let us now describe how one could calculate the production of photons when the initial ensemble is a 
non-interacting statistical ensemble. This system is implemented by taking p{ti) = exp(— /3i/o)> with a 
time path that has only the two horizontal branches. 

Derivation of the results of 1 1971 [l98l 

The phase-space density of photons at time t in the system is given by: 

2P03Sr = TTT^T? y. Tr ( pit)al{p)ax{p)) = 777^-^ V Tr ( p{t,)al{t,p)ax{t,p)) ,(72) 



'd^xd^p ~ (27r)3 V ^ ^Kvv-AV/^^-Av/^;; - ^^^^^ y ^ . 

pol A pol A 

where p{t) is the time-dependent density operator, a^(p) is the photon creation operator in the Schrodin- 
ger picture and V is the volume of the system. In the second equality, we use the cyclicity of the trace in 
order to transpose the time dependence from the density operator into the number operator, and we have 
defined: 

ax{t,p) = e^^(*-*>)aA(p)e-^^(*-*') . (73) 
This time-dependent creation operator is related to the Heisenberg field via: 

alit,p) = -ij d'x e-'P-^ d.o ^^i^)4iP) , (74) 

where e^{p) is the appropriate polarization vector. Inserting Eq. d74t in Eq. d72l . taking the time deriva- 
tive, and using the equation of motion of the Heisenberg fields, we obtain^*^: 

2po^^^ = j2^v Epoi A 1^2/0-^0 eTr(^p{t,)[J''{xo,p)idy, - ipo)A^{yQ,p) 

+{dy, + ipo)^"(yo,p) J''(xo,p)]) , (75) 



where J'^{xo,p) is the spatial Fourier transform of the current i/j{x)j'^il){x) (ip is the Heisenberg fermio- 
nic field). One can also write Jf^{xQ,p)A'^ {yo,p) = P J'^{xQ,p)A'^{yQ,p), where P denotes the path 
ordering along the Keldysh contour, and the indices ib indicate on what branch of the contour a field is 
kept. At this point, one can expand the Heisenberg fields in terms of the free fields of the interaction 
picture. Since Tr {p{ti)A^^) = 0, one must at least expand one order further in the electromagnetic 
coupling e. Truncating the expansion at the first non-zero order in e, but keeping all orders in the strong 
interactions, we obtain the following formula 



d^k 
(2^ 



^*The spatial coordinates iiave been Fourier transformed. 



+{dy,+ipo) (P A^^^_{yo,p)Af^{zo,k)) (f Jt^^4xo,p)Jain{zo,k)e'^c^?n^'''j 



(76) 



where is the QCD part of the interactions, and where C is the Keldysh contour that runs from ti 

to +00 and then back to tj. (• • •) denotes the average over the initial ensemble, and we have used the 
fact that we can factorize the photon fields and the quark fields because the initial density operator does 
not couple photons and quarks. The first of these correlators contains a factor {27r)^S{p + k) which is 
used to integrate out the vector k. Then the second correlator brings a factor (27r)'^5(0) which should be 
interpreted as the volume of the system, and cancels the factor l/V. At this point, it is a simple matter 
of algebra to work out the derivatives of the photon-photon correlator, and we get: 

dN^ >^ ^ ^ f^'^ duj sm{{po - uj){xo - ti)) 



PO-UJ 



X [(1 + n?^{po))KULO,p) - nO(po)n^l+(^,p) 



(77) 



where ri°(po) is the photon distribution in the initial ensemble, and where Il'^^{io,p) is the Fourier 
transform of the current-current correlator, possibly including QCD corrections to all orders. The stan- 
dard assumption is to assume that the system is small enough so that photons do not accumulate in the 
system. Therefore, in the "photon production rate" one drops all the terms proportional to nS^(po)> i-^- 
the photon absorption and the blocking effects. Note that there are memory effects in this formula since 
the rate at time xq depends on the initial time ti. Moreover, the self-energy under the integral over lu is 
evaluated off-shell as long as the time difference xq — ti remains finite. The limit of infinite initial time 
ti — > —00 is then obtained by using linia^+oo sin(aA;)/A; = TT6{k), and one recovers trivially the usual 
formula l93ll94l for the photon production rate, with an on-shell self-energy. 

Note also that the result derived here is the same as the one obtained in ri97l 11981 . even if the 
starting point is enunciated differently. This implies that the model of I1 97t d98§ is equivalent to a 
description in which a non interacting statistical ensemble would be prepared at t = —00, followed by 
a non-interacting evolution of the system until the finite time ti, at which point one switches on all the 
interactions. Let us also emphasize that the same calculation can be performed while keeping explicitly 
the vertical branch ([ti,ti — i(3]) of the contour, in which case the residual ti dependence drops out and 
one obtains the standard equilibrium result without having to take a limit. This is consistent with the 
fact that keeping the vertical branch of the time path just amounts to have p{ti) = exp(— (i.e. an 
interacting initial statistical ensemble) instead of p{ti) = exp(— /3ffo) (^ non-interacting initial statistical 
ensemble). 

On the approach of Il l99ll20()ll2(m 

The papers 1 1991 12001 12011 do not discuss photon production, but develop a formalism (called 
"projected functions formalism") that can be used to calculate Green's functions of an out-of-equilibrium 
system. This formalism also assumes that the initial conditions are set at some finite time ti, which is 
taken to he ti = for simplicity. Assuming also the initial ensemble to be a bath of non interacting 
particles, the time path is a Keldysh contour that starts at ti = 0, goes to t = +00 and then back to 
ti. For this reason, all the Green's functions are always evaluated with time arguments in the range 
< xq < +00. For a 2-point function with time variables xq and i/q, one can introduce the usual 
combinations sq = xq — yo and Xq = (xq + yo)/2, and one has the constraint: —2Xq < sq < 2Xq. 
This can be enforced in the Wigner transform (Fourier transform with respect to the relative variable) as 
follows: 

/+00 
dPoPxo{Po,Po)Goo{po,p) , (78) 
'OO 

where Gxq is the Wigner transform restricted to the range sq G [—2Xq, 2Xq], Goo is the unrestricted 



Wigner transform and Pxg is a projection operator given by: 



Px.{po,Po)- • (79) 

Since the time evolution starts at tj = 0, all the diagrams must be evaluated with time integrations at the 
vertices restricted by factors 9{xq). In Fourier space, this brings factors like: 

(80) 



where the po,i are the energies carried by the propagators attached to a particular vertex (all defined to be 
outgoing), instead of the usual delta functions that would ensure energy conservation at the vertex. 

The phase space density of photons at the time xq is defined^^ in 1120 II as: 

where Gs = + G |_ is the third non-zero component of the matrix propagator (the other two 

components being respectively the retarded and advanced propagators) in the Keldysh basis of the real- 
time formalism (we use here the standard notations of 12131 . and we refer the reader to this reference 



for more details). One can then use the formalism of 11991 12001 12011 in order to evaluate the photon 
production rate: 



uj^ — Pq + 2i£UJ 



,(82) 



where I[s{Lii,p) = Ils{io,p) — sign{Lo)[Ilji{uj , p) — n^(u;,p)]. At the time of writing this report, it is 
not clear whether Eqs. dTTl and (l82l are equivalent or not. 



9.14 Infinite vacuum contributions and other problems... 

We now discuss some of the consequences of Eq. dTTl . Because the self-energy Il'^_{u!,p) can be 
evaluated off-shell^° in Eq. dTTl . we can have a non zero photon production rate already from the lowest 
order self-energy, i.e. the diagram with only one quark loop. This contribution was evaluated in L1971 
11981 . where one can read: 

dN^ _ 67re2 r+°° duj sin((po - ^^){xo - k)) 
^^d^xd^p (27r)3 TT Po-^ 

X I ^{2[l-(p-fc)(p-q)]nq((7)(l-nq(A;))J(a; + A;-g) 

+ [1 + (p • • q)]nq{q)n^{k)6iu; - k - q)^ (83) 

where we denote k = p + q, p = p/p, q = q/q and k = k/k, and where nq{k) = l/(exp(/3A;) + 1) is 
the quark distribution in the system. 

The first shortcoming of this formula is that it is incomplete: its authors have dropped a term whose 
integrand is proportional to (1 — nq(g))(l — n(^{k))5{uj + q + k). This term can be non zero only for 
negative values of uj, but since the integration variable uj and the energy po > of the produced photon 

^'This formula is written by analogy with a similar relation known to be valid in thermal equilibrium. However, 12011 does 
not give a first principles justification of this formula in the case of a system out-of-equilibrium. 

"""Note that Eqs. Ml\ and J82> contain the usual on-shell contributions. They start when the photon polarization tensor is 
evaluated at the 2-loop order. 



do not need to have the same sign, the contribution of this term to the photon rate is different from zero. 
Moreover, this contribution is divergent at any photon momentum p. Indeed, the result of the integration 
over q for this term behaves like at large uj, so that the integral over uj is not defined. Physically, it 
corresponds to the decay of the vacuum into a qq pair and a photon, and is due to the mismatch between 
the vacuum of the free theory and the vacuum of the interacting theory. In other words, the initial (non 
interacting) vacuum is an excited state in the interacting theory, and tends to decay spontaneously into 
particles. Naturally, this problem disappears if ti — > — oo since sin((po — ^){xo — ti))/{po — u) — > 
'^^{po — ^) in this limit, and has no support at negative uj. Note also that the same problem arises with 
Eq. (l82b : here also, there are infinite vacuum contributions. 

Even if, following 11971 11981 . we decide to ignore the vacuum contributions, Eq. (l83l still has 
another very serious problem: the asymptotic behavior (large po = \p\) of the power spectrum is too 
hard for being integrable. More precisely, one has: 

dN-y 2 

which implies that the integrated energy yield per unit time and per unit volume / d^ppo{dN^ / d'^xd^p) 
is infinite because the integral does not converge in the ultraviolet. 

It has been suggested that the divergence of the total energy might be due to the fact that the 
coupling constant in this model is switched on instantaneously. Indeed, even if the self-energy II'^I {lu, p) 
vanishes exponentially when \uj\ —>■ +oo, the power law decrease of the power spectrum is due to the 
fact that the function sin((po — (^){xo — ti))/{po — ^) has itself a power law decrease with respect to 
Pq. However, this function arises as a result of the integration over the intermediate time z q: 



sin((po - uj){xo - tj)) _ 



dzQ cos((po - ^){zo - Xq)) 

dzo9{zo - ti) cos(po - ^){zo - xq) , (85) 

oo 

where the cosine comes from factors like limj;o_,a.Q ((9j^Q — ipQ){P Af^^{yo,p)Af^{z o,k)) in Eq. S76i . 
and where the 9{zq — ti) reflects the time-dependence of the coupling constant (zq is the time attached 
to a vertex coupling the photon to a quark line). As a toy model, one can smoothen the behavior of the 
coupling constant by replacing 6{z o — ti) by 1/ (exp((tj — zq)/t) + 1) where r is the typical time during 
which the coupling constant evolves from to its normal value. Doing this substitution in Eq. (l85l gives: 

cos((po -oj){zo - Xq)) _ 27rrsin((po - uj){xo - ti)) ^^-{xo-u)/t 

° e(ti-zo)/T _^ I ~ g7r(po-w)T _ g-7r(po-w)T X ^ _ cu)'^T-2 

— oo 

+O(Te-2(^'0-*«)/^) , (86) 

which is a good approximation for xq — ti » r. In this expression, the first term vanishes exponentially 
when Ipo — 1^1 ^ +C)o. The power law behavior has been restricted to the second term, which has an 
exponentially suppressed prefactor when xq — ti S> r, i.e. when one looks at the rate at a time which 
is far from the region where the coupling constant is changing. In particular, one can see that the power 
law behavior becomes dominant over the exponential behavior only for |po — '^l ^ {xq — tj)/7rr^. In 
other words, there is a considerable range of po — lo in which one does not see effects of the power law 
behavior. However, this does not prevent the ultraviolet divergence when one tries to integrate Eq. (l83l 
over the photon energy. 

This residual ultraviolet divergent power spectrum seems to be due to the fact that Eq. (l83b eval- 
uates the photon number of an interacting system with a non-interacting definition of the photon. As a 



consequence, this definition of the photon number "measures" also the virtual photon clouds that sur- 
round all charged particles after the electromagnetic interaction has been turned on. This interpretation 
is also supported by the fact that the photon self-energy in Eq. (I83t is always evaluated off-shell. 

9.15 What have we learned? 

It appears that the models considered in l UTWlfTOSl and ri99"200"2011 suffer from very serious problems 
which forbid any quantitative prediction at this point. In particular, a common problem in these models is 
that they lead to infinite vacuum contributions, which correspond to the spontaneous decay into particles 
of the vacuum of the non-interacting theory. Note that in a more recent paper |214J, an ad-hoc subtraction 
is proposed in order to dispose of all these divergences^ ^ What remains to be seen is whether this 
subtraction can be justified from first principles, and whether one can get rid of the inherent arbitrariness 
that comes with subtracting infinities. 

From a more formal point of view, the important lesson taught to us by the study of this model 
is that in generic non-equilibrium problems, the initial ensemble should be specified with respect to the 
spectrum of the interacting theory. An approach exploring this issue has been proposed in L215.I . 

9.2 The 2PI effective action approach 

J. Serreau 

To follow the space-time evolution of an out-of-equilibrium system of quantum fields one may 
solve the Schwinger-Dyson equations of the theory - written in the form of an initial value problem - 
for given initial conditions. This cannot be done exactly and, in practice, one has to rely on approxi- 
mation schemes. A powerful way of deriving systematic approximations to the dynamical equations is 
to work at the level of effective actions, from which one obtains the equations of motion by functional 
differentiation with respect to the fields. An example is given by the well-known one-particle-irreducible 
(IPI) effective action, which is the generating functional of IPI Green's functions. In particular, the use 
of such functional methods ensures that the global symmetries of the underlying theory are preserved by 
the approximate equations of motion (for instance, this guarantees energy conservation). 

It has however been observed that standard approximations of the IPI effective action can be sec- 
ular in time [ 216[l217ll : the validity of a given expansion in a small parameter - such as a coupling or a 
1/A^-expansion - can be spoiled by the occurrence of terms which grow indefinitely with time. This has 
been a major obstacle for practical studies of the real time dynamics of far from equilibrium quantum 
fields beyond so-called mean-field approximations (leading order in large-A^, Hartree). The latter, which 
neglect direct scattering between (quasi-)particles, are free of secular terms and have been extensively 
used over the last decades in various physical situations I 218II219I . However, these "collisionless" ap- 
proximations are known to fail to describe important physical effects such as late time thermalization, 
or early time damping of unequal-time correlations (see e.g. |220|). Another important case where one 
can perform explicit calculations of the non-equilibrium dynamics corresponds to the classical statistical 
field theory limit I 2211l222ll223l . This can provide a good description of the quantum dynamics when the 
typical occupation numbers of the field modes are large. A great advantage is that one can solve exactly 
the full nonlinear classical dynamics by means of standard Monte Carlo methods together with numer- 
ical integration techniques. However, there are important physical questions, such as the description of 
quantum thermalization at late time, which cannot be addressed with these methods. 

Important progress have been made in recent years with the use of approximation schemes based 
on the two-particle-irreducible (2PI) effective action |224, 225). The latter is a generating functional 
for the correlation functions of the theory, parameterized in terms of the connected one and two-point 

^'Although this is not said explicitly in 12141 . these new results completely invalidate the results and the phenomenological 
discussion of I197||198I . 



functions, i.e. the average value of the field and the propagator (for comparison, the IPI effective ac- 
tion is parameterized in terms of the one-point function only). This approach allows for practicable 
and systematic calculations of the non-equilibrium dynamics beyond mean-field and classical field ap- 
proximations.^^ In particular, it has been possible to demonstrate for the first time the late time quantum 
thermalization of a scalar field in 1 + 1 dimensions from a first principle calculation, by using a three-loop 
approximation of the 2PI effective action 12271 . 

9.21 Current status 

In recent years, the 2PI effective action approach to out-of-equilibrium phenomena has been extensively 
investigated and applied to various situations of physical interest in the context of scalar and fermionic 
field theories (for a recent review see 12281 ). The most studied case is that of scalar field theories, where 
the question of reliable approximations has received a lot of attention. Various schemes have been worked 
out and used e.g. to study thermalization. These include a coupling-expansion up to three-loop order, 
which is the simplest approximation including direct scattering and describing thermalization |227|, as 
well as a 1/A^-expansion up to next-to-leading order (NLO) 12 20. 2291 . In a slightly different context, 
namely the Schwinger-Dyson approach, motivated by a direct truncation of Schwinger-Dyson equa- 
tions, a similar approximation has been studied: the so-called bare vertex approximation (BVA) |230|. 
In terms of the 2PI 1/A^-expansion, the latter partially include field-dependent^'^ NNLO contributions 
[ 2291 . These approximations have been compared to exact Monte Carlo simulations in the framework of 
classical statistical field theory in 1 + 1 dimensions 12301 123 ll 12321 . It has been shown that the 2PI 
expansion at NLO provides a quantitative description of the dynamics already for moderate values of N. 
For instance, the damping rate of unequal-time two-point functions is accurately described at NLO for 
N > A [231 1 . Similar studies of the time evolution of the one and two-point function in the limit = 1 
123011232)1 show that, although quantitatively different, both the NLO approximation and the BVA are in 
qualitative agreement with the exact result near the symmetric regime. Studies of the quantum dynamics 
demonstrate that the 2PI 1 /A^-expansion at NLO does not exhibit any spontaneous symmetry breaking 
for A^ = 1 in 1 + 1 dimensions |233 1, in agreement with general results. Although it is not a controlled 
limit for a 1/A^-expansion, such a qualitative agreement provides a very sensitive test of this systematic 
expansion scheme. For instance, the inclusion of part but not all NNLO contributions, as in the BVA, 
leads to (spurious) spontaneous symmetry breaking in one spatial dimension |233|. This may, however, 
be particular to the limit A^ = 1 and a direct comparison between these two approximations, both in 
classical and quantum field theory, shows that they agree rather well with each other for A^ > 4 |234|. 
First results in 3 + 1 dimensions exhibiting the phenomenon of spontaneous symmetry breaking have 
been presented in L228J . using the 2PI 1/A^-expansion at NLO. 

The 2PI technique has been demonstrated to be a powerful tool to study realistic particle physics 
applications. In particular, it has been used to perform the first quantitative study of the phenomenon 
of parametric resonance in quantum field theory beyond mean field approximations 12351 . This phe- 
nomenon provides a paradigm for situations with nonperturbatively large particle densities, where nei- 
ther gradient nor coupling expansions are applicable. In contrast, the 2PI 1/A^-expansion remains valid 
in such situations. In particular, the use of the NLO approximation have been shown to solve the problem 
of an analytic description at large densities. A full numerical solution of the corresponding equations of 
motion as well as approximate analytic results concerning the nonlinear dynamics have been obtained 
f2351. Another important point which has been studied in details in the context of scalar field theories 
concerns the renormalization of 2PI approximation schemes |23^, 237, 238 ). It has been shown that ap- 
proximations based on a systematic loop-expansion of the 2PI effective action (the so-called ^-derivable 
approximations) can be renormalized at any truncation order if the theory under consideration is pertur- 

^^This approach, supplemented by a gradient expansion, also provides a very efficient way to derive Boltzmann equations in 
quantum field theory |226|. 

In particular, the BVA and the NLO approximation are equivalent in the symmetric regime (see e.g. l229l ) 



batively renormalizable f236\. Finally, the 2PI approach has been successfully applied to the study of 
fermionic field theories. In particular, the late time thermalization of fermionic quantum fields has been 
demonstrated for the first time from a first principle calculation II239I . This has been done for a realistic 
3+1 dimensional theory of Dirac fermions ("quarks") coupled to scalars ("pions") in a chirally invariant 
way, employing a perturbative expansion of the 2PI effective action at lowest non-trivial (two-loop) or- 
der. Bose-Einstein and Fermi-Dirac distributions have been shown to emerge from the non-equilibrium 
dynamics without further approximation. 

One of the main open questions of the 2PI effective action approach concerns the description of 
the dynamics of gauge fields. Gauge invariance makes it a non trivial issue: the 2PI effective action being 
a functional of the one and two-point functions only, it is difficult to find truncations which are consistent 
with Ward identities. In order to deal with this problem, different directions have been followed. One of 
them is to try to enforce Ward identities at a given truncation order by modifying the 2PI scheme (see e.g. 
12401). A similar strategy has proven useful for the calculation of thermodynamic properties of a QCD 
medium in equilibrium [241], where a hard thermal loop approximation has been applied on the top of a 
2PI two-loop truncation. The price to pay is that one may loose the description in terms of an effective 
action, which, as emphasized above, is very useful for non-equilibrium systems. Another interesting 
possibility is to analyze the magnitude of the gauge-fixing dependence of physical results within a given 
approximation scheme. This has been investigated within the 2PI loop-expansion in Ref. [242. .243.1 and 
further developments are to be expected. 

9.22 Application to photon production? 

From very general considerations, the problem of photon emission can be reduced, at lowest order in 
the electromagnetic coupling constant acm, to the calculation of the following connected current-current 
correlator 



where J^(t, a;) is the current which couples to the photon field and where the brackets denote the average 
with respect to the (initial) density matrix of the emitting system. The 2PI approach provides in principle 
a powerful tool to compute the time dependence of the above two-point function. The current is generi- 
cally a bilinear in the fields describing the emitting system and the correlator Eq. |^can be expressed in 
terms of the corresponding four-point function. The latter can be obtained from the 2PI effective action 
by taking appropriate functional derivatives. More details concerning the possible use of these methods 
in the context of out-of-equilibrium electromagnetic radiation can be found in Ref. |215|. 

As emphasized above, the 2PI effective action approach allows one to compute the two-time de- 
pendence of the current-current correlator without making a priori assumptions about this dependence. 
For example, contrarily to the gradient expansion, it does not require any separation of scales between 
the {t + t') and the {t — t')-dependence of two-point functions. In the context of heavy ion collisions, 
there are various situations where one does not expect any such separation of scales. This is for example 
the case when the system rapidly cools down through a phase transition. Other typical non-equilibrium 
effects which can play an important role are initial time and finite time effects. First, the presence of 
the initial condition alone breaks time-translation invariance and it is known that the early time behav- 
ior cannot be described by a gradient expansion. Of course, because of interactions, the system forgets 
about the initial condition on a time scale governed by the inverse damping rate. Only in the case where 
the total photon production rate (integrated over the whole time history of the system) is dominated by 
times larger than the latter can one neglect the presence of the initial time. Moreover, because invariance 
under time-translation is broken, off-shell processes may contribute to the production rate, as discussed 
in the previous sections. The 2PI approach may be a way to address these issues, taking these effects into 
account in a consistent manner. 

Other effects which can be of importance, are those related to "internal" aspects of the time evolu- 
tion of the emitter. An immediate example is provided by situations where the lifetime of some physical 




(87) 



excitations of the system cannot be neglected on the timescales under consideration. This could be the 
case for a gas of resonances and may also be of importance when chiral symmetry gets restored. Just 
as the off-shell contributions mentioned above, these effects would automatically be taken into account 
in the 2PI effective action approach, applied to the appropriate microscopic theory. On the contrary, 
they are known to be very difficult to include in a gradient expansion. Other interesting topics include 
for example the question of photon emission from a disoriented chiral condensate [244.] . which may be 
formed during the out-of-equilibrium chiral phase transition. 

In conclusion, the 2PI effective action approach has been much developed in recent years in the 
context of non-equilibrium quantum field theory and has been demonstrated to provide a powerful tool 
to study realistic physical situations. These methods might prove useful for the study of photon and 
dilepton production in general non-equilibrium situations. 

10. LATTICE CALCULATION OF THE VECTOR SPECTRAL FUNCTION 

P. Petrezcky, F. Gelis, G.D. Moore 



In this report we are interested mostly in the production of energetic real or quasi-real photons 
and the calculations of rates is done in the HTL framework. The same technics allows the computation 
of the production rate of massive static photons 12451 which have recently also been considered using 
non-perturbative lattice technics. Based on existing litterature there appears a contradiction between the 
two methods: in the lattice approach a strong threshold effect is observed and the production rate tends 
to when the photon mass decreases, in contrast to the HTL approach where the rate is increasing. In 
this section we discuss these matters further. 

Lattice calculations provide information on the imaginary time current-current correlator 

G,..(r,p) = (J^(0,p)jt(T,-p)), (88) 

where J^(r,p) = J2x e*^'^ a;) = J2x e*^^ ^(■'"j 3;)7^g(r, a;). This is related to the spectral 
function a^^, by an integral relation. 
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(89) 



where the kernel K{t, lo) = cosh(a;(T — 1/2T))/ sinh(u;/2r)) is the imaginary time free boson prop- 
agator . The spectral function a^,^{uj,p) can be related to the emission rate of real or virtual photons 
|93l|9l. 



From the lattice data on G^j,^ (r, p) the spectral function can be reconstructed using the Maximum 
Entropy Method (MEM) 12461 12471 . The first calculations of a^^ at finite temperature using MEM were 
done in |248|. There it was found that a^^{uj) = X]/i=i/ ~ ^) slightly enhanced over the 
free correlator in the region AT < lo < 7T and strongly suppressed for to < 3T. These features of the 
spectral function are clearly illustrated by Fig. [83j where results at T = 3Tc are shown. Also shown 
there is the one-loop order HTL resummed result for the spectral function showing strong enhancement 
for small lo. One concern about this result is that the reconstruction only uses A'^^ = 16 points in the 
imaginary time direction. However, as noted in |248|, the correlator G^y{T,p) alone already provides 
stringent constraints on a^i,{u>,p) and thus allows to check the validity of the perturbative calculations 
of this quantity. In particular, the value of the correlator at r = 1/(2T) is not sensitive to the behavior 
of a^^{LO,p) for uj > 16T = 1/a (for = 16 the lattice spacing is 1/(16T)) and thus to the lattice 
artifacts L248J . The lattice data on G^^(a;,p) provide an integral constraint (sum rule) on (T^^(u;,p), 

/ du;a^,{u^,p) ^ = G^^ir = l/{2T),p). (90) 
Jo smn.[-^) 



'■2T 

^'*On the lattice, this should be replaced by its lattice version 

''^The two-loop corrections increase the one-loop result even more 12491 . 




Fig. 83: The vector spectral function ai_i^{uj) at T = STc. The error bars in Fig. l83l are statistical and were estimated using the 
method of t246ii247il . The systematic errors will be discussed latter in the text. 

In table 5. the ratio of G^^{t = l/(2T),p) calculated on a 64^ x 16 lattice to the corresponding 
free (but lattice) value is listed for several p = \p\- Note that Gfj_^{T = l/{2T),p) is always close to 
the corresponding value in the free theory. For example, for p = it is only 9% larger than the free 
value. This implies that a strong enhancement of the vector spectral function relative to its lowest order 
perturbative (free) value is ruled out. 



T/T, 


p = 


p = i.57r 


p = 2.22T 


p = 3.UT 


p = 4.UT 


1.5 


1.095(19) 


1.079(14) 


1.052(20) 


1.001(23) 


0.898(85) 


3.0 


1.083(14) 


1.071(10) 


1.060(14) 


1.039(10) 


1.007(15) 



Table 5: The value of G(r = l/{2T),p)/Gfree,iatt{r = l/{2T),p) at T/Tc = 1.5 and 3.0 

Recently, lattice calculations of Gii{T) using unimproved Wilson fermions^^ on anisotropic lat- 
tices [2501 have helped address the question about the small number of temporal points used to recon- 
struct (Tjj(w) (only the spatial part was considered in this study). The use of anisotropic lattices allowed 
to use more points {N^ = 32 — 96) in the imaginary time direction while keeping the temperature fixed. 
On the other hand the use of the unimproved Wilson action for quarks introduces additional lattice ar- 
tifacts. The vector spectral functions obtained in |250| qualitatively agree with those obtained in 12481 
for uj < 16T. They have a broad peak around uj/T ~ 5 — 6 and suppression^' for u < 3T. For 
T ~ (1.4 — 1.5)Tc {Nt = 54) the position of the peak of the vector spectral function roughly coincides 
with the findings of fMEi . 

The drop of the spectral function below = 3T is in contradiction with HTL perturbation theory, 
which predicts al/u behavior for a^^{io) at small u (see Fig. l83b . due to the contribution of processes 
like bremsstrahlung. However, it is expected that this lowest order behavior is modified into l/y^ 
by multiple scattering effects when u becomes small. At even smaller lo, this behavior will be further 
reduced by strong dissipative effects, so that a^^{uj) vanishes linearly with u. Indeed, the small u limit 
of e^a^^/Quj is the electrical conductivity in the quark-gluon plasma. The small uj behavior of the MEM- 
extracted spectral function would imply a vanishing electrical conductivity, but at the moment we lack 
a natural explanation that would lead to a very small electrical conductivity in a quark-gluon plasma. 

""^In 1 248 1 a non-perturbatively improved fermion action was used, i.e. discretization errors of order 0{a) present for Wilson 
fermions were completely removed. 

At T ~ 2Tc strong enhancement of the spectral function over the free spectral function is observed; however, it is not clear 
how statistically significant this enhancement is. 



Related to this, there are also questions about whether a^^{uj), or crjj(cj) (just the spatial part) is the 
appropriate quantity to reconstruct. 

There are systematic uncertainties in the reconstructed spectral function. The most important 
uncertainties apart from those coming from the limited number of points in the time direction (Nr) are 
those associated with discretization errors and with the choice of the so-called default model m{Lo), the 
prior assumption on the form of a^^{uj). For large lu, a^^ can be calculated using perturbation theory 
(due to asymptotic freedom) and one has (t^^{uj) = rn-ow^ for uj » IGeV with mo ^ 0.076 I251II . It 
is natural to build this assumption into the MEM analysis by choosing the default model m(u;) = ruQio'^ 
with niQ = 0.076. The above form of a^f^{uj) (or m{u)) ) is valid only in the continuum; on the lattice it 
will certainly be modified and this fact should be taken into account when choosing the default model. 
The simplest way out, suggested in 12461 12471 . is to use mow^ for the default model but to vary mp and 
study the dependence of the result on rriQ. We note that the result varies little if one varies the default 
model at small uj. For instance, one can consider m{uj) = rriQU}'^ ian\i{uj / {4T)) or itiquj'^ + niiuj with 
mi ~ T and find only small changes in the corresponding spectral function. The result is sensitive to 
the choice of mo (this is not the case at T = 12461 12471 ). In Fig. |84]we show the vector spectral 
function a^^{uj) using m{u}) = mow^ with mo = 0.152, i.e. with twice the perturbative value. The 
gross structure of the spectral function remains unchanged; however, the statistical significance of the 
suppression for small u is considerably reduced. One expects that dependence on the default model 
should become weaker if larger A'^^ are considered, but it is unlikely that it will disappear completely. 
The most practical way to solve this problem is to fix the large u behavior of o"^^ at T = and use this 
as the default model, as at large uj a^^ is expected to be temperature independent. 

0.4 
0.35 

0.3 
0.25 

3 

~~-=L 0.2 
0.15 
0.1 
0.05 





Fig. 84: The vector spectral function at T = STc. The error bars are statistical and estimated according to 12461 12471 . For 
further details see text. 

Before closing this section, we should mention an alternative approach to the problem of extracting 
the behavior of the vector spectral function at zero momentum and small energy, proposed recently by 
S. Gupta [252]. This method tries to extract a Pade approximation of the energy dependence of the 
spectral function from the lattice data for the vector propagator. The extracted slope of the spectral 
function at l<j = gives the following value for the electrical conductivity: 

C7ei«7.0(^e2)r, (91) 

s 

where is the sum of the electric charges squared of all the species present in the plasma. The 

discrepancy between this value of the electrical conductivity and the apparently very small value one 
would obtain from the results of 1 248 1 has not been elucidated at the time of writing this report. 



Hopefully, once all these systematic effects are under control, the above issues can be sorted out 
and the analysis of the vector spectral function can be extended to finite momenta to give an estimate on 
the emission rate of quasi-real photons. 

11. APPENDIX I - HARD SCATTERING CROSS SECTIONS: FROM pp TO pA AND A+A 

D. G. d'Enterria, J. Ranft 



11.1 Proton-nucleus (pA) collisions 

Glauber formalism 

The inelastic cross-section of a p+A reaction, apA, can be derived in the eikonal limit (straight line 
trajectories of colliding nucleons) from the corresponding inelastic nucleon-nucleon A^A^ cross-section, 
am, and the geometry of the pA collision simply determined by the impact parameter h of the reaction. 
In the Glauber multiple collision model I.253J . such a cross-section reads 



I — Q-'^NN'^A(b) 



(92) 



where TA{b) is the nuclear thickness function (or nuclear profile function) of the nucleus A at impact 
parameter b: 

TA{b) = [ dzpA{b,z). (93) 



Tyi(6) gives the number of nucleons in the nucleus A per unit area along a direction z separated from the 
center of the nucleus by an impact parameter b. The nuclear density, pA{b, z), is usually parametrized 
by a Woods-Saxon distribution with nuclear radius Ra = 1-19 • A^^^ — 1.61 • A^^/^ fm and surface 
thickness a = 0.54 fm as given by the experimental data II254I and normalized so that 

d^bTA{b) = A. (94) 
Hard scattering cross-sections 

Though Eq. ^ is a general expression for the total inelastic cross-section, it can be applied to an 
inclusive p + A —>■ h + X process of production of particle h. When one considers hard scattering 
processes (e.g. direct photon or high-p-r vr'' production), the corresponding cross-section fx^^'^'^ is small 
and one can expand Eq. d92l in orders of amTA (b) and then, to first approximation 



^hard ^ I ^hard 
^pA " 



d^a^Z'^TAib) (95) 



"Minimum bias" hard scattering cross-sections 

Integrating Eq. (l95l over impact parameter, and using (l94l . one gets the minimum bias (MB) cross- 
section for a given hard process in pA collisions relative to the same cross-section in pp (or A^A^) colli- 
sions: 

K^r v= ^-4^'' (96) 

From this expression it is easy to see that the corresponding minimum bias multiplicity (invariant yield 

jhard ^hard /^9^° 'i 

^NN.pA — "NN.pA/ "NN.pAJ 



per nuclear reaction: Nj^^^f^ = '^nnpa/'^nndA^ ^ given hard-process in a pA collision compared to 



that of a pp collision is 



/ ]\Thard\ a ^NN nrhard ^ hard tQl\ 

V^pA 1MB — ^ ■ geo ' NN ~ geo ' '^NN ' ' > 

^pA ^pA 



where a^^^ is the geometrical pK cross-section given, in its most general form, by Eq. ( I92t . The average 
nuclear thickness function for minimum bias reactions [making use of Eq. (l94l l reads: 



. _ JdHTA A A 

{Ta)mB = n ,2, - ^ - -geo- (98) 



Thus, for a pPb {A{Pb) = 208) collision at LHC energies ^/s^ = 8.8 TeV with 

aNN ~ 77 mb |255l, and (99) 
^ 2162 mb 12561, (100) 

one obtains: {N^pI^'^)mb ~ 7.4-A'^^'''^, and the average nuclear thickness function amounts to {Tph)MB = 
0.096 mb-i = 0.96 fm-^. 

11.2 Nucleus-nucleus (A+B) collisions 

Glauber formalism 

As in the proton-nucleus case, the inclusive inelastic cross-section cjab for a collision of nuclei A and B 
is given in the multiple-scattering Glauber approximation by: 



1 _ Q-^NN'^ABi^) 



(101) 



where now TAB{h) is the nuclear overlap function of the nuclei A and B separated by impact parameter 
h. TAB{b) can be written as a convolution of the corresponding thickness functions of A and B over the 
element of overlapping area (Ps {s = (x, y) is a 2-D vector in the transverse plane, and b is the impact 
parameter between the centers of the nuclei): 

TAB{b) = j (fsTA{s)TB{\b-s\). (102) 

TAB{b) is normalized so that integrating over all impact parameters one gets: 

j (fbTAB{b) = AB. (103) 

Hard scattering cross-sections 

As in the pK case, for hard processes of the type A + B ^ h + X, Eq. dlOlt . can be approximated by: 

^A^'^ j dHa^i'-^TABib). (104) 

"Minimum bias " hard scattering cross-sections and yields 

Integrating Eq. (I104t over impact parameter and using (11031) . one gets the minimum bias {MB) cross- 
section for a given hard process in A-i-B collisions relative to the corresponding pp cross-section: 

KDMfi= A-B-afr^^" (105) 

Again the corresponding minimum bias multiplicity (invariant yield per nuclear reaction: N^'^^ab ~ 
^nnab/'^nnab^ ^'^^ ^ given hard-process in a A-i-B collision compared to that of a pp collision is 

I ■\Thard\ /I O '^AW ^jhard A - B hard f^c\f.\ 

V^AB )mB — A - U ■ —g^ ■ iV^^ — —geo- ' ^ NN ) (^^O) 
^AB ^AB 



where o"^^° is the geometrical A+B cross-section given, in its most general form, by Eq. dlOU . The 
average nuclear overlap function for minimum bias reactions [making use of Eq. (Il03t 1 reads now: 



, , _ fd^bTAB A-B AB 

[Tab) MB = f-WT — — —FT, , p NO — ~~ged, (107) 

j d^b tt^Ra + RbY cTab 

Thus, for a Pb+Pb {A'^{Pb) = 43264) collision at LHC energies = 5.5 TeV with 



aNN ~ 72 mb [2571, and (108) 
'^PbPb - 7745 mbl257i, (109) 

one gets: {N^Sl,t)MB ~ 400 • N^1[\ and the average nuclear overlap function amounts to {Tpi,ph)MB = 
5.59 mb-i = 55.9 fm'^. 

For comparison for RHIC, at 200 GeV, one has aj^/f^ 40.83 mb and ct^Iau ~ 7085 ± 33 mb. 



Binary collision scaling 

For a given impact parameter b the average hard scattering yield can be obtained by multiplying each 
nucleon in nucleus A against the density it sees along the z direction in nucleus B, then integrated over 
all of nucleus A, i.e. 

where we have made use of expressions d93b and (I102t . In the same way, one can obtain a useful 
expression for the probability of an inelastic NN collision or, equivalently, for the average number of 
binary inelastic collisions, {Ncou), in a nucleus-nucleus reaction with impact parameter b: 

{Ncoii)ib) = am-TAB{b) (111) 

From this last expression one can see that the nuclear overlap function, TAB{b) = Ncou{b) / (ym [mb^^], 
can be thought as the luminosity (reaction rate per unit of cross-section) per AB collision at a given 
impact parameter. From (II lOt and (llllb . we get so-called "binary (or point-like) scaling" formula for 
the hard scattering yields in heavy-ion reactions: 

{b) « (iV,o«) (6) • N^?,"' (1 12) 



Hard scattering yields and cross-sections in a given centrality class 

Equation (11041) gives the reaction cross-section for a given hard process in A-i-B collisions at a given 
impact parameter 6 as a function of the corresponding reaction cross-section in pp collisions. Very 
usually, however, in nucleus-nucleus collisions we are interested in calculating such a reaction cross- 
section for a given centrality class, (o'^'b'^)ci-C2' where the centrality selection Ci — C2 corresponds 
to integrating Eq. (I104t between impact parameters 61 and 62- It is useful, in this case, to define two 
parameters 1.258. 2591 : 

• The fraction of the total cross-section for hard processes occurring at impact parameters bi < b < 

62 {d'^b = 2TTbdb): 

2-K 

fhard{bl < 6 < 62) 



AB 



b2 



bdbTABib). 



The fraction of the geometric cross-section with impact parameter 61 < 6 < 62: 



2tt 



bi 



bdb 1- c-^nnTabW 



I 9eo 



AS ' 



/geo(&l < & < 62) ■ 

[fgeo simply corresponds to a O.X (e.g. 0.1) factor for the X%(10%) centrality.] 



(113) 



(114) 



Hard scattering production is more enhanced for increasingly central reactions (with larger number on 
Ncou) as compared to the total reaction cross-section (which includes "soft", - scaling with the number 
of participant nucleons Npart as well as "hard" contributions). The growth with b of the geometric 
cross-section is slower than that of the hard component. For this reason, the behaviour of Jab and fgeo as 
a function of b, although similar in shape is not the same (see 125811 ): fhard ~ 1 for 6 = 2Rj^, but /^go ~ 
0J5forb = 2RA. 

Similarly to (I107t . we can obtain now the nuclear overlap function for any given centrality class 

!t''dHTAB_A.B hard 



AB)Ci-C2 



£ d% ^Ib heo 

The number of hard processes per nuclear collision for reactions with impact parameter 6i < 6 < &2 is 
given by 



(115) 



(TV; 



hard\ 

AB /C1-C2 



a' 



hard/ 



AB 



ll <b < 62) 



A-B 



^hard f 

^ NN J hard 



a 



geo 
AB 



f, 



(116) 
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which we could have just obtained directly from dllOt and ( 11151 . From ( 11061 ) and dl 161 ) it is also easy to 
see that: 
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hard\ 
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hard \ 
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hard 
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Finally, the cross-section for hard processes produced in the centrality class Ci 
ing to a fraction fg^o of the reaction cross-section) is: 



fhard 

17~ 



hard 

m 



(117) 

C2 (correspond- 
(118) 



geo 



Figure [S5J extracted from ii258il . plots the (top) fraction of the hard cross-section, fhard{^ < 6 < 62) 
(labeled in the plot as /as), as a function of the top fraction of the total geometrical cross-section, 
/geo(0 < b < 62), for several nucleus-nucleus reactions. As a practical application of Eq. (11181) and 




geo 



Fig. 85: Figure 4 of ref. 12581 . Fraction of the hard cross-section, fhard{0 < b < 62), vs. the fraction of the total geometrical 
cross-section, fgeoiQ <b < 61), for several heavy-ion collisions (from left to right): 197+197, 110+197, 63+193, 27+197, and 
16+197. 



the results of Fig. [85j the hard-scattering cross-sections in Pb+Pb for the top 0-10% (fhard = 0.41 for 



/^go = 0.1 from the practically equivalent Au+Au system of figurel85l and 0-20% (fhard = 0.664 for fg^o 
= 0.2) central collisions relate to the pp cross-section, in the absence of nuclear effects, respectively as: 

«b'\-io% = (208)2 . ^ . « 1.7 • 10^ • 4r (119) 
«b'\-20% = (208)2 • ^ • 4r ^ 1-4 • 10^ • a';,%' (120) 

A straightforward way to compute the invariant yield for a given hard process in a given centrality 
class of a nucleus-nucleus collision from the corresponding yield in pp collisions consists in determining, 
via a Glauber MC calculation, the number of inelastic NN collisions corresponding to that centrality 
class via 

{Ncoll)Ci-C2 = {Tab)ci-C2 ■ ^NN , (121) 
and then use this value in the "binary-scaling" formula 

(iVj^7'^)ci-c. = {N,ou)c,-C2-Nli^'''', or (122) 

{(^AB )Ci-C2 = {Ncoll)ci-C2 ■ -(^NN (123) 

O'NN 

The same two formulae above apply to pA collisions (of course substituting A-i-B by pA and computing 
NcoU from Ta instead of from Tab)- 

Finally, to obtain the experimental rates, (A/^^'^^)ci-C2> actually measured in a given centrality 
bin one needs to take into account the expected integrated luminosity Cint [mb~^] as follows: 

{■^Ab'^)Ci-C2 = ^int ■ icr\B"^)ci-C2 (124) 



11.3 Hard scattering yields and cross-sections for pPb and Pb+Pb collisions 

As a practical application of the Glauber approach described here, in Table |6l the values of {Ncou) are 
quoted for different centrality classes obtained from a Monte Carlo calculation |260| for pPb {^/snn - 
8.8 TeV and am - 77 mb) and Pb-i-Pb i^/sNN - 5.5 TeV for an inelastic pp cross-section of am = 72 mb) 
collisions (Woods-Saxon Pb density parametrization with Ra = 6.78 fm and a = 0.54 fm). Using (I122t . 
(I123t and Table |6l we can now easily get the scaUng factors of the cross-sections and yields from pp to, 
e.g., central (0-10%), minimum bias, and semi-peripheral (60-80%, from the combined average 60-70% 
and 70-80%) pFh (8.8 TeV) and Pb-^Pb (5.5 TeV) colUsions : 

For pPh colUsions {a^pl = 2162 mb): 
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For Pb-nPb colUsions (af^"^ = 7745 mb): 



{N'pA-10% = 1670 • iV^r 
(NptPh) 60-80% = 29.1 • nI^^'''^ 
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hard 
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Table 6: Number of inelastic NN collisions, (Ncoii), and nuclear thickness {Tppt} or overlap (Tptpt) function per centrality 
class, in pPb (^/sj^ = 8.8 TeV, cr^^y = 77 mb) and Pb+Pb collisions at LHC (^fsfff^ = 5.5 TeV, cr;^;^ = 72 mb) obtained with 
the Glauber Monte Carlo code of ref. ^60|. The errors, not shown, are of the same order as the current uncertainty in the value 
of the nucleon-nucleon inelastic cross section, cryvyy, at LHC energies (~ 10%). 



Centrality (Ci - C2) 



p+Pb 







{Tpph) (mb-i) 


{Ncoll) 


{TpbPb) (mb-i) 


0-5% 


15.7 


0.203 


1876.0 


26.0 


0-10% 


15.3 


0.198 


1670.2 


23.2 


10-20% 


13.8 


0.179 


1019.5 


14.2 


20-30% 


12.0 


0.155 


612.4 


8.50 


30-40% 


9.9 


0.128 


351.8 


4.89 


40-50% 


7.8 


0.101 


188.0 


2.61 


50-60% 


5.6 


7.27-10-2 


92.9 


1.29 


60-70% 


3.8 


4.93-10-2 


41.4 


5.75-10-1 


70-80% 


2.6 


3.37-10-2 


16.8 


2.33-10-1 


80-90% 


1.7 


2.20-10-2 


6.7 


9.31-10-2 


90-100% 


1.2 


1.55-10-2 


2.7 


3.75-10-2 


min. bias 


7.4 


9.61-10-2 


400.0 


5.58 



Pb-FPb 



11.4 Nuclear effects in pA and A+B collisions 

Eqs. ( l95l and (I104t for the hard scattering cross-sections in pA and A-i-B collisions have been derived 
within an eikonal framework which only takes into account the geometric aspects of the reactions. Any 
differences of the experimentally measured cr^^'^^^ with respect to these expressions indicate "de facto" 
the existence of "nuclear effects" (such as e.g. "shadowing", "Cronin enhancement", or "parton energy 
loss") not accounted for by the Glauber formalism. Indeed, in the multiple-scattering Glauber model 
each nucleon-nucleon collision is treated incoherently and thus, unaffected by any other scattering tak- 
ing place before (initial-state) or after (final-state effects) it. 

If the Glauber approximation holds, from ( l95t and (I104t one would expect a ex A^, and oc 
growth of the hard processes cross-section with system size respectively. Equivalently, since N^^^'^g = 
^nnab/^nnaa ^NNAB ^i*^ A^^^, onc would expect a growth of the number of hard 

process as oc A^/^ , oc A^/^ for pA^ AA collisions respectively. Experimentally, in minimum bias pA and 
AB collisions, it has been found that the production cross-sections for hard processes actually grow as: 

{4T)mb = A^ - cT^JJ^ , and (ct^T')mb = [ABT - ^^^^ , with a + 1 (131) 

More precisely, in high-pT- processes in pA and heavy-ion collisions at SPS energies one founds a > 1 
(due to initial-state pr broadening or "Cronin enhancement"); whereas a < 1 at RHIC energies ("high- 
er suppression"). Theoretically, one can still make predictions on the hard probe yields in pA^ AB col- 
lisions using the pQCD factorization machinery for the pp cross-section complemented with the Glauber 
formalism while modifiying effectively the nuclear PDFs and parton fragmentation functions to take into 
account any initial- and/or final- state nuclear medium effects. 



11.5 Summary of useful formulae 

Finally, let us summarize a few useful formulae derived here to determine the hard-scattering invariant 
yields, cross-sections, or experimental rates, from pp to pA and A+B collisions for centraUty bin C\ — C2 
(corresponding to a nuclear thickness Ta or nuclear overlap function Tab and to an average number of 
NN inelastic collisions {Ncou})'- 

dprdy A,AB Ci C2 (j^prpdy 

(j2 hard\ £ hard 

1 — 1 = {Ta,ab)c\-C2 ■ ^vAM ■ r , (133) 

dpTdy dpTdy 

id'K:Al)Cr-C2 . ,^ . aeo d'a^pf .,,4. 
= C,r.f{TA,AB)c,-C2-<^.A.As-^^ (134) 

= (N^oiih -c ■ (135) 
dpxdy ^ ^ dpxdy 



dprdy am dprdy 

{d m,a.ab)ci~c2 _ . ,^ > (j;;^.AB ^ V .,^7. 

1 1 — ^int ■ V^coll)Ci-C2 ■ ■ T~ (13') 

dprdy am dpxdy 



12. APPENDIX II - STANDARDS FOR LUMINOSITIES AND ACCEPTANCES 

The list of considered collision systems is: Pb+Pb, D+Pb or pPh, Ar+Ar and pp. 

Luminosities: in (cm~^ s~ 1) 
Pb+Pb : 5 * 10^6 
D+Pb : 5 * 10^^ 

pFh : 10^^ 

Ar+Ar : 10^^ but 3 * lO^'^ for ALICE central part 
pp : 10^'^ but 3 * 10^° for ALICE 

Center of mass energy 

The center of mass energy of a coIUsion of nucleus {Al, Zl) with nucleus {A2, Z2) is given by s/s = 
14 TeV * ^yzi * Z2/{A1 * A2) 

CMS Acceptance 

For electron, muon, photon and jet reconstruction, essentially no hadronic particle identification (PID) is 
available. In heavy ion colhsions, CMS understands tracking for \r]\ < 1.5, px > 4: GeV. The tracking 
studies for charged hadrons at p^" > 1 GeV is in progress. This should be a conservative standard. On 
the other hand, the principal CMS rapidity acceptance is: 

I??! < 2.4 for muons and charged hadrons, 

1 77 1 < 3. for photons and electrons, 

Ir/I < 5 for jets. 

Depending on event multipUcity and further studies, a larger acceptance [\r]\ < 2.5 and px > 1-5 GeV in 
the forward directions] may be available for heavy ion colhsions. These numbers may be used as more 
realistic standard. 

ALICE Acceptance 

For hadrons and electrons: \r]\ < 0.9, 

PID pion pT > 100 MeV 

PID kaon, proton pr > 200 MeV 

HMPID small acceptance RICH allows PID of p, K for pr <5 GeV, \r}\ < 0.5 and 57 degrees 
azimuth 

PHOS photon spectrometer covers \ri\ < 0.13 and 100 degrees azimuth 
For muons: 2.5 < t] < A, pT > 1-0 GeV 

Acceptance for charmonium and bottonium reconstruction from electron pairs: |?7| < 0.9, > 0. 
For more details about this acceptance and trigger conditions, contact experimentalists. 
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